IRE 


Transactions 
ANTENNAS and PROPAGATION 


Volume AP-6 APRIL, 1958 Number 2 


Published Quarterly 


TABLE OF CONTENTS 


Sonik, Bstedh WiRRREB cig gn 5 Sasa oe ao oho Sei Ae aa ee 159 
CONTRIBUTIONS 

Mieke or curbuleat Mixing in Seatter Propagation... ............ 6.8.56. e sess es Ralph Bolgiano, Jr. 161 
The Influence of Moisture in the Ground, Temperature and Terrain on Ground Wave Propagation in the 

VO SIN SIA ASIIGI SAS bo SNA Sica Oe Lea 7 ene Be ee Oe ee B. Josephson and A. Blomquist 169 
Distance Dependence, Fading Characteristics and Pulse Distortion of 3000-MC Trans-Horizon Signals... . 

es ee eed Line ches ae cig et ool or Ol eee nce B. Josephson and G. Carlson 173 
Some Microwave Propagation Experiences from a “Just-Below-Horizon” Path.................0.2220-5- 

EE Cee trys Scryer ooh sw ewes B. Josephson and F. Eklund 176 
Sarees nme OMGGhiN eT VLCC es, ce cs Bre seseeipcO te sl eh eee He esis eee he Renee dye ales Paul I. Richards 178 
Determination of a Current Distribution over a Cone Surface Which Will Produce a Prescribed Radiation 

eee er GO ee. ices gone nis ope ane Wa sw Dilees cade nein wee NS H. Unz 182 
The Effects of the Physical Parameters on the Bandwidth of a Folded Dipole....................... ae 

NR OEP AEE Od Lay Syria ga diy ese vee lv hate ee slg es J. P. German and F. E. Brooks, Jr. 186 
thhe Radiation Characteristics of a Zig-Zag Amtenna..... <2... 1... eee ee eee Dipak L. Sengupta 191 
A Note on the Effective Aperture of Electrically Scanned Arrays. ..........---.5--+555: R. W. Bickmore 194 
The Characteristic Impedance of Two Infinite Cones of Arbitrary Cross Section.......... Robert L. Carrel 197 
Microwave Stepped-Index Luneberg Lenses.........-.....-5.-00-05: G. D. M. Peeler and H. P. Coleman 202 


COMMUNICATIONS 


5 Measurements of £250-MC Scatter Propagation as Function of Meteorology...........--++.++++sseeeee 

reer ere ieee viele wee oe ee D. L. Ringwalt, W. S. Ament, and F. C. MacDonald 208 
Design Data for Small Annular Slot Antennas...........-.--.-+.+++->: W. A. Cumming and M. Cormier 210 
a Gio cease sip sc Rete epee ee et eee eee PiU 


PUBLISHED BY THE 


ofessional Group on Antennas and Propagation 


Administrative Committee 
J. I. Bohnert, Chairman 
R. L. Mattingly, Vice-Chairman 


H. G. Booker F. T. Haddock, Jr. W. H. Radford 
Arthur Dorne J. W. Herbstreit J.B: Smyth 
J. W. Findlay D. D. King O. G. Villard, Jr. 
R. K. Moore 
Ex Officio Members 
P. S. Carter D. C. Ports A. H. Waynick 
Honorary Member 
L. C. Van Atta 


IRE TRANSACTIONS® PGAP IS A QUARTERLY PUBLICATION 
DEVOTED TO EXPERIMENTAL AND THEORETICAL PAPERS ON 
ANTENNAS AND WIRELESS PROPAGATION OF ELECTROMAGNETIC WAVES 


MANUSCRIPTS should be submitted to John B. Smyth, Editor, Smyth Research Associates, 3555 
Aero Court, San Diego 11, Calif. Manuscripts should be original typewritten copy, double spaced, plus 
one carbon copy. References should appear as footnotes and include author’s name, title, journal, vol- 
ume, initial and final page numbers, and date. Each paper must have a summary of not more than 200 
words. News items concerning PGAP members and group activities should be sent to the News Editor, 
Mr. Arthur Dorne, Dorne and Margolin, Inc., 30 Sylvester Street, Westbury, L.I., N.Y. 


ILLUSTRATIONS should be submitted as follows: All line drawings (graphs, charts, block diagrams, 
cutaways, etc.) should be inked uniformly and ready for reproduction. If commercially printed grids 
are used in graph drawings, author should be sure printer’s ink is of a color that will reproduce. All half- 
tone illustrations (photographs, wash, airbrush, or pencil renderings, etc.) should be clean and ready 
to reproduce. Photographs should be glossy prints. Call-outs or labels should be marked on a registered 
tissue overlay, not on the illustration itself. No illustration should be larger than 8 x 10 inches. 


Copies can be purchased from 
THE INSTITUTE OF RADIO ENGINEERS 
1 East 79 St., New York 21, N.Y. 


PRICE PER COPY: members of the Professional Group on Antennas and Propagation, $1.15; 
members of the IRE, $1.75; nonmembers, $3.45. 


ANNUAL SUBSCRIPTION PRICE: PGAP members, included in PGAP fee of $4.00; IRE mem- 


bers, $8.50; Colleges and public libraries, $10.00; non- 
members, $17.00. 


IRE Transactions oN ANTENNAS AND PROPAGATION 
Copyright © 1958, by The Institute of Radio Engineers, Inc. 
Entered as second-class matter, at the post office at Menasha, Wisconsin, 


Acceptance for mailing at a special rate of postage is provided for in the 
in Paragraph 4, Section 412, P, L 


under the act of August 24, 1912. 


act of February 28, 1925, bodied 
. & R., authorized October 26, 1927, ae ia ak 


“1968 


ah a 


SA ae ae ee 


% A REPORT FROM THE CHAIRMAN 


_ I should like to summarize briefly what has occurred 
recently in the PGAP and what is planned for the im- 
mediate future. There has been very little communica- 
tion between administration and members, so that a 
“state of the union” report appears opportune. 
__ The big problem in the PGAP since 1956 has been a 
lack of adequate funds. This state of affairs was created 
when our Group published the papers of the Electro- 
magnetic Theory Symposium held at the University of 
Michigan, which wiped out our financial reserves and 
created a backlog of papers awaiting publication in our 
TRANSACTIONS. Since then, the administrative com- 
mittee and officers of the PGAP have spent an inordi- 
“nate amount of time trying to obtain sufficient funds 
for proper operation. These efforts resulted in the neg- 
ect of some other matters, which I shall mention later 
as needing action. 

_ Therefore, it is with some pride, and relief, that I am 
able to report that this backlog of papers was removed 
by publication in the January issue. The necessary 
“money was raised by the sale of more than $6000 worth 
of advertising in our TRANSACTIONS. Delmer Ports and 
Art Dorne are primarily responsible for this work, ac- 

complished despite the arbitrarily high rates estab- 
‘lished by the IRE Executive Committee. This means 

that papers accepted by our reviewers can now be 
published promptly in the following issue of TRANSAC- 
TIONS. 

I do not wish to convey the idea that our financial 

troubles are over. A glance at the following financial 
‘statement shows that the PGAP was still destitute at 


the close of 1957. 
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news and views 


For the Period from January 1, 1957 to December 31, 1957: 


Balance from January 1, 1957............... $ 3,448.13 
Receipts During Period: 
IRE Matched Funds......... $ 2,404.00 
Assessments). 5.0.0. 060004005 9,950.71 
Advertisingysse- 6 ote. eS ya. 2,192.60 
Sale of Publications.......... 1,329.35 
Surplus from Meetings........ 272.32 
Othew Sourcesiy soma stor —0- 
PotaliReceiptsnac. woes Fos cy ste cl) ek $16,148.98 
Total (Original Balance and Receipts)...... $19,597.11 
Expenses During Period: 
Rublications. 00. 0 kta eae $19,514.15 
Membership Service Charges. . 318.12 
Others ee NaAeitern cen. ners 2 —0- 
otal Expenses tire omranvic ty tito sae $19 , 832.27 
Balance as of December 31, 1957............. $ 235.16 Credit 


Our credit balance will change to debit as money is 
received from advertising already sold. However, more 
advertising needs to be sold in order that a working 
reserve can be built up. Any member wishing to help 
should contact Delmer Ports, as described on the back 
pages of recent issues of the TRANSACTIONS. 

A very interesting statistic, and evidence of the vigor 
of the PGAP, is that we have published 2278 pages in 
these TRANSACTIONS from 1951 to 1957 inclusive. This 
is over 700 pages more than any other professional 
group has published in the same period. This record is 
even more impressive when one remembers our financial 
difficulties and realizes that there are eight other pro- 
fessional groups larger than we, two of which are over 
twice as large. 

Membership in the PGAP is increasing rapidly, 
further proof that our Group is healthy and dynamic. 
Here are the figures for our growth since 1952. 


160 
Paid Membership End of Year 
980 1952 
1219 1953 
1276 1954 
1470 1955 
1997 1956 
2406 1957 


In addition, there are about 400 student members. 
There are eleven PGAP chapters as follows: 


Chapter Chairman 
Akron George H. Welch 
Albuquerque-Los Alamos Dick Moore 
Boston J. T. deBettencourt 
Chicago N. J. Sladek 


Herman V. Cottony 
Robert J. Stegen 
Glenn A. Scharp 


Denver-Boulder 
Los Angeles 
Orange Belt 


Philadelphia E. J. Forbes 
San Diego David Proctor 
Syracuse William T. Whistler 


Washington, D. C. Harry Fine 


A chapter is being formed in San Francisco and a joint 
chapter with the PGMTT is being organized on Long 
Island. Certainly, the PGAP is an organization which 
can be recommended to your friends and coworkers. 

As in the past, the PGAP has sponsored technical 
sessions this year at the National Convention in New 
York, at the joint IRE-URSI Spring Meeting in Wash- 
ington, and at WESCON in Los Angeles. Another 
national symposium is planned in October in Washing- 
ton under the leadership of Harry Fine. The general 
topic is propagation in the troposphere and ionosphere, 
and by meteor trials. 

There are several current items which require action 
on the part of the Administrative Committee of the 
PGAP. At our next meeting I will propose that the 
scope of our Group be expanded to include officially 
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the field of Radio Astronomy. This is recognition of 
what we are already doing. Such action will require the 
approval of the IRE Executive Committee. I will pro-_ 
pose also that an award program be initiated within 
our Group, e.g., awards for best papers in the TRANS- 
ACTIONS or at meetings. I will continue to urge a news- 
letter, a much-needed means of communication among 
members as to chapter activities, PGAP participation 
in technical meetings, current actions of the Adminis- 
trative Committee, and other news items of general 
interest. Such a publication should aid in obtaining new 
members. And, of course, the financial question will 
come in for its share of attention. The present conditions 
under which advertising in the TRANSACTIONS is sold | 
are unsatisfactory. A change in policy by the Institute 
is needed whereby the advertising rates are made 
reasonable and continuity of sales is maintained. Every 
effort to improve this situation will be made. 

In summary, I can say that both the present state — 
and the future of the PGAP appear good; indeed, even 
better than in the past few years as attested to by the 
facts submitted earlier. Also, the future of the pro- 
fessional group system appears bright, with the groups 
assuming more responsibility and leadership in the 
Institute. To that end, I refer the reader to the article 
written by the President of the IRE, Donald Fink, in 
the September, 1957 issue of PROCEEDINGS. 

J. I. BOHNERT 
Chairman, PGAP 


NEWSLETTER 


This is the last News and Views column. It will be 
replaced shortly by a Newsletter, and is expected to be 
much more useful than News and Views has been. The 
lead time required to get it out will be much shorter. | 
Additionally, it will appear more frequently: every other 
month. 

I will continue as PGAP News Editor. Accordingly, 
please send news items to me, Arthur Dorne, 29 New 
York Avenue, Westbury, Long Island, N. Y. 


1958 
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contributions 


The Role of Turbulent Mixing in Scatter Propagation* 
RALPH BOLGIANO, JR.t 


; Summary—tThe Villars-Weisskopf-Wheelon theory describing 
_ turbulent mixing of an established gradient is shown to contain a 
_ contradiction which necessitates its being discarded. To fill the gap 
thus created, the theory of isotropic mixing is extended to account 
_ for the presence of a gradient. The results indicate that “mixing-in- 
_ gradient” cannot be employed to explain the wavelength dependence 
_ characteristic of much of the radio data. On the other hand, it is 
_ shown that experimentally determined spectra of refractive index 
fluctuations lend strong support to the mixing theory herein set forth. 
_ Hence, the conclusion is reached that scatter theory, as it is currently 
4 based on atmospheric turbulence, can provide, at best, an incomplete 
description of transhorizon propagation. 


I. INTRODUCTION 


. EVERAL years ago Batchelor,! on the one hand, 
S and Villars and Weisskopf,? on the other, came to 
the conclusion quite independently that turbulent 
mixing of established mean gradients of refractive index 
must play a major role in the generation of the index 


* Manuscript received by the PGAP, September, 14, 1957. Por- 
tions of this paper were presented at the URSI-IRE meeting in 
Washington, D. C., May 25, 1957, under the title “The spectrum of 
turbulent mixing and its application to scatter propagation.” This 
research was supported by the Wright Air Development Center under 
contract AF 33(616)-3236. : 

+ School of Elec. Eng., Cornell University, Ithaca, N. Y. 

1G. K. Batchelor, “The Scattering of Radio Waves in the Atmos- 
phere by Turbulent Fluctuations in Refractive Index,” School of 
Elec. Eng., Cornell University, Res. Rep. EE 262; September 15, 
1955. 

2F. Villars and V. F. Weisskopf, “On the scattering of radio 
waves by turbulent fluctuations of the atmosphere,” Proc. IRE, vol. 
43, pp. 1232-1239; October, 1955. - 


fluctuations which are believed to account? for a signifi- 
cant, if not the principal, portion of signal propagated 
well beyond the horizon. In their same paper Villars and 
Weisskopf proposed a theoretical description of such a 
turbulent mixing process which leads to a fluctuation- 
size dependence for the scattering coefficient o of 
L,’. (L,=X/[2 sin 6/2], the scattering eddy-size.) For 
the troposphere this yields o proportional to X. 
Wheelon!® subsequently was led to the same result by 
independent reasoning, basically dimensional in charac- 
ter. Since there is a large and growing body of data*~® 
which, when considered within the framework of scatter 


’ See for example: H. G. Booker and W. E. Gordon, “A theory of 
radio scattering in the troposphere,’’ Proc. IRE, vol. 38, pp. 401-412; 
April, 1950. 

4A. D. Wheelon, “Radio frequency and scattering angle depend- 
ence of ionospheric scatter propagation at vhf,” J. Geophys. Res., vol. 
62, pp. 93-112; March, 1957. 

* A. D. Wheelon, “Spectrum of turbulent fluctuations produced 
by convective mixing of gradients,” Phys. Rev., vol. 105, pp. 1706- 
1710; March 15, 1957. 

6 L. G. Trolese, “Characteristics of tropospheric scattered fields,” 
Proc. IRE, vol. 43, pp. 1300-1305; October, 1955. 

7 J. H. Chisholm, et al., “Investigations of angular scattering and 
multipath properties of tropospheric propagation of short radio waves 
beyond the horizon,” Proc. IRE, vol. 43, pp. 1317-1335; October, 
1955. 

8 K. Bullington, et al., “Results of propagation tests at 505 mc and 
4,090 mc on beyond-horizon paths,” Proc. IRE, vol. 43, 1306-1316; 
October, 1955. 

9 J. H. Chisholm, et al., “Experimental investigations of the angu- 
lar scattering and multipath delays for transmissions beyond the 
horizon,” presented before Combined Technical Session, URSI, 
Washington, D. C., May 23-25, 1957. 
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phenomena, has been interpreted’ as being in agreement 
with such a \ dependence, many investigators have re- 
cently been inclined to accept scatter propagation 
theory, based on “mixing-in-gradient,” as rather well 
established. 

Unfortunately, Wheelon’s analysis contains a con- 
tradiction which, as it is fully appreciated, casts grave 
doubt upon the entire L® law. It is the intent of this 
paper, first, to set forth clearly this inconsistency. 
Secondly, by extending Batchelor’s isotropic theory to 
the case of turbulent mixing in the presence of a gradi- 
ent, it will be shown that the Z1* form, previously pre- 
dicted by both Batchelor and Silverman," is still ap- 
propriate for a large range of fluctuation sizes. Further, 
it will be brought out that, if the presence of a gradient 
of mean refractive index has any effect, it must be to 
weaken (make the exponent smaller) this dependence 
at large scales, not strengthen it. Thirdly, by analyzing 
refractometer data directly, support will be introduced 
for the contention that the refractive index fluctuations 
themselves have a mean square spectral density which 
is in consort with the L"/? law. Finally, three methods of 
calculating values of o will be discussed. 


II. INCONSISTENCY OF THE L> Law 


In Section 4 of his paper‘ Wheelon assumes a priori 
that, in the inertial subrange, the spectrum of mean 
square fluctuations is independent of the molecular dif- 
fusivity D, the wave number k(=27/L) being taken to 
be much smaller than k,, the viscosity cutoff. By purely 
dimensional arguments he then arrives at his L® law. 
However, this implies the spectrum of molecular “dissi- 
pation,” 2Dk*E,,(k), is proportional to k—. (E,(k) is the 
spectral density per unit volume of mean square fluctua- 
tions of refractive index and is directly related to the 
comparable spectrum of, say, moisture vapor density 
or electron density fluctuations. Z,(k) is synonymous 
with Wheelon’s B(k).) Thus, the distribution of molec- 
ular effects has its maximum at a wavenumber much 
smaller than the range of values he is considering. There- 
fore, molecular diffusion cannot be neglected in this 
range, in contradiction with is a priori assumption. 

It is immediately evident the difficulty is not confined 
to Wheelon’s analysis. An L> form implies a k-* mean 
square fluctuation spectrum in the corresponding 
wavenumber range, regardless of the method by which it 
was derived. It then follows, as above, that molecular 
diffusion is a prominent agent in this range, which is 
contrary to the premise (Villars and Weisskopf?) that 
this is the inertial subrange. Consequently, the L5 law 
must be abandoned as a theoretical description of mix- 
ing-in-gradient. 


%° See for example: K. A. Norton, ef al., “The use of angular dis- 
tance in estimating transmission loss and fading range for propaga- 
tion through a turbulent atmosphere over irregular terrain,” Proc. 
IRE, vol. 43, pp. 1488-1526; October, 1955. 

4 R. A, Silverman, “Turbulent mixing theory applied to radio 
scattering,” J. Appl. Phys., vol. 27, pp. 699-705; July, 1956. 
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II]. TuRBULENT MIXING OF GRADIENTS 


In order to pave the way for the theoretical treatment 
which follows, first consider the problem of turbulent 
mixing of a mean gradient of refractive index in a 
heuristic manner. The arguments to be presented here 
must not be accepted in lieu of the more precise ones of 
the next section but only as aids in visualizing the entire 
process. 


A. Qualitative Description of Mixing 


At all scales present in the velocity field the eddies do, 
on the average, mix-in fluctuations from the mean gradi- 
ent proportional to their respective sizes. But the local, 
instantaneous gradient on which the smaller eddies 
operate is, on the average, more intense than the mean 
gradient as a result of the stretching of fluid lines and 
surfaces by the turbulence, thereby bringing equivalue 
surfaces closer together. Hence, the smaller eddies 
mix-in progressively more and more than their share 
figured on a mean gradient basis. This is a manifestation 
of the “internal mixing” or inertial transfer mechanism, 
the conversion of fluctuations of one size to others of 
smaller scales. This transfer across the spectrum (see 
Fig. 1) continues until the local gradients become so 


ie 
ae eae 
Yasir BY “INTERNAL MIXING" 
td ean i 
SAT 


DISSIPATION BY MOLECULAR DIFFUSION 


Fig. 1—Conceptual model of the mixing, transfer, and diffusion 
processes in turbulent mixing. 


intense that molecular diffusion acts to dissipate them. 
In the intermediate range, in which inertial transfer is 
the principal process at work, a statistical steady-state 
is reached. Then one has a simple relation between the 
spectral density, E,(k), the net rate of transfer of mean 
square fluctuations across k, S,(k), and the “lifetime” of 
a given size eddy, 71, 


Srlk)rz ; 
rape: LS (kb) tx. (1) 


E,(k) ~ 
The universal equilibrium theory of turbulence! leads, 


in this inertial subrange, to 


rr & Lt, 


% See, for example, G. K. Batchelor, “The Theory of Homogeneo 
ay ee Cambridge University Press, Cambridge, Boer p. 114 
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Thus, in the region where S,(k) is approximately con- 
stant, the mean square fluctuation spectrum is propor- 
tional to L*/*; and at larger scales, where the effects of 
mixing-in-gradient should be more pronounced, the rate 


_ of transfer, S,, increases with k (decreases with L) so 


_ that the spectral density will depend upon L according 


to 


Bath) ery. Vaso soe (2) 
These arguments have much appeal. Nevertheless, 
there are fundamental questions which frequently arise; 


the following two are important examples. What evi- 


_ dence is there that the local gradients are, in fact, more 


intense than the average gradient and that, at the smal- 


- ler scales, the inertial transfer mechanism therefore 


ft x 


cee ee 


ote he 


dominates the mixing-in process? Is there basis for the 
assumption, on which the above reasoning squarely 
stands, that it is mean square fluctuations which are 
conserved in this inertial transfer? 


_B. Analytical Development 


In order to find answers to these questions and to 
place the whole development on a more substantial 
foundation, consider explicitly the interaction of the 


_ velocity field with a property of the atmosphere prin- 


cipally responsible for variations of the index of refrac- 
tion, say, moisture vapor density in the lower tropo- 


sphere, or electron density in the ionosphere. Let the 
-symbol y, a scalar point function of position and time, 
_represent such a density. Then, from the conservation 
_of matter (assuming the fluid is dynamically incompres- 


sible) it is known that y satisfies 


oy 


Uj; 
Ox; 


ary 


O%X,0X> 


oy 
at 


(3) 


Here u;(=u) is the vector velocity, Dy is the molecular 
diffusivity, and standard tensor notation has been em- 
ployed—a repeated index means summation over all 
the vector components. This equation shows that the 


time rate of change of w is determined by the net 
_ effect of convection and molecular diffusion, neglecting 
such factors as evaporation and condensation since 


we are later going to consider only steady-state con- 
ditions. Writing y as the linear sum of its average 
value and a fluctuation component 6y, 


(x, ) = V(x) + OV, 4), 


the overbar indicating time, and/or ensemble, average, 
and similarly 


u;(x, t) Ss 1;(x) os 5u;(x, t), 


and then subtracting the average of (3) from (3) itself 
yields a conservation equation for the fluctuation com- 


ponent, 
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076 
Se, 4 
ymax (4) 


On the reasonable assumption that mean values of the 
fluctuation components are nearly independent of posi- 
tion the term 6u,(06W/dx;) in (4) is of smaller order of 
magnitude than the other terms and may be neglected. 
This may be said generally of spacial derivatives of 
mean values of the varying components and is tanta- 
mount to making the usually satisfactory assumption 
of local homogeneity. 

1) The Relative Intensity of Local Gradients: Multiply- 
ing each of the remaining terms in (4) by 26y and again 
taking averages leads to a statement of the conservation 
of mean square fluctuations in which the significant 
terms are 


ee O6y\? 
— pi/0K: ( “) 
ot Ox; Ox; 


= 2K,(grad J)? — 2Dy(grad Y)?, —_(5) 


where —6y6u; has been expressed in terms of the eddy 
diffusivity, Ky, and the mean gradient. The first term 
on the right represents the average generation of y~- 
fluctuations by interaction of the turbulent velocities 
with the mean gradient, while the last term describes 
the dissipation of these fluctuations by molecular action. 
Under steady-state conditions such that the time de- 
rivative is zero, it becomes apparent that the magnitude 
of the gradient of the fluctuation component far exceeds 
that of the mean value since it is well known empirically 
that the eddy diffusivity is several orders of magnitude 
greater than the molecular diffusivity. This then is the 
evidence sought in answer to the first question in Section 
III-A above. 

2) The Conservation of Mean Square Fluctuations: 
Restricting the consideration to a finite volume V and 
substituting for y and uw; in (4) in terms of their trans- 
forms in three-dimensional Fourier space, 4y(k) +64 y(k, 
t) and A;(k)+64,(k, t), respectively, it may be shown, 
assuming 0Y/dx; to be constant in V, (see Appendix [). 


d| 5Ay(k, ¢) |? ow 
AEC = — 2@{5A,(k, 16 Ay*(k, )} — 
ot Ox; 
+ Oy(k, k’)dk’ — 2Dyk®| 6Ay(K, 4) |? (6) 
k—space 


Here R{ i stands for the real part of the quantity with- 
in the brackets, * indicates the complex conjugate, 
dk=dk,dkodk;, and Q,(k, k’) is the spectral density of 
the rate of transfer of mean square w-fluctuations to the 
Fourier component of wavenumber k from that of 
wavenumber k’. Again in this form the first and last 
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terms on the right describe the excitation and dissipa- 
tion, respectively. The extra, integral term represents 
the internal mixing, between all other Fourier com- 
ponents in the spectrum and that particular component 
being considered, of wavenumber k. On close inspection 
it is found that the function Qy is antisymmetric in its 
two arguments (see Appendix I);i.e., 


Qy(k, k’) = — Qk’, &). 


Thus the exchange term contributes nothing to the over- 
all fluctuations since its integral over all k-space will be 
zero. In turn this implies the conservation of mean 
square fluctuations in the inertial transfer process and 
provides the answer to the second question in Section 
III-A. 

3) The Transfer Function and the Fluctuation Spec- 
trum: Integrating (6) over all directions of the vector 
wavenumber k and suitably normalizing with respect 
to the volume V over which the Fourier transforms 
were defined, yields, in the steady-state, (see Appendix 
J) the relation 


aw re) 
T(t) + [Puce ea = 2D HH, 
Xi 


0 


where 


Oe) a steak eMule 
Blt) = f Tsavte, 9 haa, 


87° 
the spectral density of mean square y-fluctuations per 
unit volume (m.s. 6y/vol.). dQ(k)=elemental solid 
angle in k-space. Once more the first term on the left is 
the excitation, the second is the exchange, and the term 
on the right is the molecular dissipation. Now, integrat- 
ing (7) from a given value of wavenumber to infinity 
one has, taking account of the antisymmetry which Py 
inherits from Qy. 


-) ow 2 
2D, [ kh’? Ey (Rk) dk” — =f Ty (kh) dk" 
k OX; 5, 
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The following properties of Sy() should be noted. 
a) If there is a large k-range over which 


if "WH aR ip “REE (R dR’ (9a) 


and 


f "Tapl!))dk"))) if “Tap AR, (9b) 


S,(R) is sensibly constant in this range. This corresponds 
to the inertial subrange in the universal equilibrium 
theory. 

b) If, for the eddy sizes of interest, molecular dif- 
fusion is of secondary importance; i.e., Léuz,/Dy>1, 
where 6uz, is the rms turbulent velocity characteristic 
of L-size eddies, (9a) above is satisfied for the corre- 
sponding k range. Then, over the range of interest, Sy(R) 
is either constant, or as a result of the gradient of J, an in- 
creasing function of k. 

c) Over that range of & such that, in any small wave- 
number interval, the total of mean square y-fluctuations 
produced by the excitation and removed by the dissipa- 
tion is of a smaller order of magnitude than that ex- 
changed among the various Fourier components by the 
inertial transfer process—the cumulative nature of the 
excitation and dissipation insures this condition is sat- 
isfied well beyond the inertial subrange in both direc- 
tions—E,(k) depends only upon S,(k), S(k), and k. 
Under these circumstances, on dimensional grounds 
alone, Zy(k) and S,(k) are related in general magnitude 
by 


Ey(k) ~ Sy(R)[S(k) H1 3-878, 


Combining (b) and (c) it is apparent that the result (2) 
reached in Section III-A by heuristic reasoning must in 
fact be true. Fig. 2 depicts the relations among the vari- 
ous spectra and the transfer function S,(k). 

4) The Scattering Coefficient: It has been shown! that 
the scattering coefficient 


(10) 


angular density of scattered power per unit volume of scatterer 


(6 at 


power density of the wave front incident upon the scatterer 


= if ak!" it Py (k", Rak’ = Sy(k). - (8) 


Thus S,(k) is the net rate of transfer of m.s. dy/vol. 
across wavenumber k. This transfer function is analo- 
gous to the corresponding kinetic energy transfer, S(k), 
which many authors® have considered in the case of 
turbulent flow, except that there decaying turbulence is 
treated rather than the steady-state situation postu- 
lated here. 


13 Tbid., pp. 126 ff. 


is related to | 54 n(k, t)|? by 
4m? sin? x| 5A,(k. t) [2 


At V 


e= (11) 
6A,(k, t) is the analysis, in 3-dimensional Fourier space, 
of the index of refraction variations within the volume 
V, x is the angle between the incident electric field at V 
and the ray to the receiver, and k, is defined in Figwae 
If |84,,(k, t)| is isotropic at ky, 


on 8r* sin? x En(ke) 
rf Re? 


o 


(12) 
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Fig. 2—Relative shapes of the spectral and transfer functions which 
play a role in turbulent mixing. (a) Excitation, (b) dissipation, 
(c) internal transfer, (d) spectral density. 
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Fig. 3—Geometry of the scatter problem. 


Fluctuations of the refractive index are, in general, 
attributable to fluctuations of one or several properties 


_ of the atmosphere (electron density in the ionosphere, 


moisture vapor concentration and temperature in the 
troposphere) each of which is conserved according to 
(3).14 Consequently, E,,() will take the form of (10) and 
(12) becomes 


824 sin® x By 


Sy (Re) [S(he) >" *ket!8, (13) 


Cm 


where By, is the conversion factor from m.s. dy to m.s. 
dn. Thus, as long as k,<k,, so that molecular diffusion is 
of little importance, o obeys an L,""/*, or weaker, law. 
The presence of a gradient of f, if it has any effect, can only 
cause S,(k) to increase with k (decrease with increasing L) 
and therefore weaken the dependence. 


IV. MEASURED REFRACTIVE INDEX SPECTRA 


There is now some direct evidence regarding the 


spectrum E,(k). Within the last several years Crain, 
and others, have used the airborne refractometer to 


14 Temperature itself is not a conservable property; heat content 


is, however. Under the assumption of small Mach number, tempera- 

ture, being a good measure of heat, obeys (3) to a close approxima- 

tion. For a more complete discussion of this matter see reference 20. 
15 


Crain, “Survey of airborne microwave refractometer 


Cy M. 


measurements,” Proc. IRE, vol. 43, pp. 1405-1411; October, 1955. 
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Fig. 4—Experimental verification of the k-5/3 form for the spectrum 
produced by turbulent mixing. Each curve is the mean of the 
normalized spectra derived from the data of individual runs as 
indicated. 


gather data on the microvariations of refractive index. 
Many of these data have been translated into mean 
square fluctuation spectra. True, they are one-dimen- 
sional spectra, ®,1(k), derived from data along a line, 
whereas E,(k) is representative of the index variations 
throughout a volume. Nevertheless, assuming isotropy, 
which is strongly supported both theoretically and ex- 
perimentally so long as the scale considered is small 
compared to the largest eddies, it may be shown (see 
Appendix II) that these two types of spectra are related 
by 


1 aOnh) _ 
iho 


d®,(k 
, aba) 


(14) 


Thus, if they are power functions of k, they are identical 
except for a multiplicative constant. Consequently, it is 
of some interest to study log-log plots of such spectra 
to determine what k dependence is indicated. Each curve 
in Fig. 4 is the mean of several individual spectra de- 
rived from data collected in a given region at a par- 
ticular altitude and time.'’ This tends to smooth out the 
random character of the individual samples. Note that 
the value of the slope is very close to —5/3, in excellent 
agreement with that predicted by (10) for the inertial 
subrange. 


V. EVALUATION OF ¢ 


There are two routes by which the value of the scat- 
tering coefficient may be calculated. The first, and more 


16 The first half of this relation has been given previously by L. S. 
G. Kovasznay, M. S. Uberoi, and S. Corrsin, “The transformation 
between one- and three-dimensional power spectra for an isotropic 
scalar fluctuation field,” Phys. Rev., vol. 76, pp. 1263-1264; October 
15, 1949. However, they did not discuss the second half, which is 
essential to the present argument. 

17 C, M. Crain, “Index of Refraction Fluctuations and 1046 MC 
Phase Difference Measurements Along an 18,000 Foot Elevated 
Path,” University of Texas, Austin, Texas, Elec. Eng. Res. Lab., 
Rep. No. 6-15; June 1, 1956. 
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direct, involves the determination of En(kz) following 
the methods of Section IV and the substitution of this 
value in (12). If, due to limitations of experimental 
technique, it is not possible to get ®,(k) at ke, it may 
still be possible to determine E£,(k) for some smaller 
wavenumber and extrapolate to the value at k,, em- 
ploying the k dependence found in the manner of Section 
VG 

The second route makes use of equation (13) and 
depends upon knowledge of the transfer function 
S,(k). Assuming , lies in the inertial subrange there are 
in turn, two ways to calculate Sy(z). Since virtually all 
the production of mean square y-fluctuations occurs at 
k<k, and similarly all the dissipation of such fluctua- 
tions occurs at k>>k,, (5) and (8) and the inverse of (8) 
(i.e., the integral of (7) from zero to k) may be combined 
to yield 

2K,y(grad py)? a3 Sy (ke) = 2D,(grad 6y)?. (15) 

The eddy diffusivity K, may be approximated by the 
product of the scale / and the rms turbulent velocity 
6u of the turbulent field 


Ky ~ lou ~ 13[S(1/D }1"8. (16) 


The inertial subrange relations, though not strictly ap- 
propriate, have been employed in the final conversion. 
Assuming further that S(1//) is approximately equal to 
S(k.), substitution in (13) gives 


1674 sin? x By 
Ww 


a~ Pa(erad p)*ke M8; (17) 


In this formulation it must be noted that, because of 
the density stratification of the atmosphere, it is not the 
absolute gradient of y but rather the deviation from an 
adiabatic profile!® which is required. That is, for (grad 
Y) in (17) one must substitute (grad Y+Ty), where Ty, 
the adiabatic lapse rate of y, is given by 


(18) 


g is the acceleration due to gravity, p is the density of 
the composite fluid, p is the pressure, and y = 1.4. 

If (grad W)?, can be determined, since it has been 
shown in Section III-B, 1) that (grad 6y)?)))(grad ~)2, 
and it therefore follows that (grad 6y)? is approximately 
equal to (grad )?, the right side of (15) may be used to 
evaluate S,(k,). The calculation of o is then straight 
forward except that a value must be assumed for S(k,), 
which corresponds to the rate of dissipation of turbulent 


18 That is, a vertical distribution of y such that a sample of the 
fluid which is transported vertically through the gravity-induced 
pressure gradient and which undergoes the corresponding expansion 
or compression adiabatically will always have a value of ¥ equal to 
the y of its surroundings. 
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energy € when &, is in the inertial subrange. Brunt’s 
estimate! that € is of the order of 5X10~4 m? sec~* may 
be employed as a first approximation. One other prac- 
tical difficulty arises with this method. As may be seen 
by comparing (5) and (8), contributions to (grad y)? are 
heavily weighted in favor of the high wavenumber end 
of the spectrum. Hence, it is unlikely that with present 
refractometers, having a resolving distance of the order 
of a meter, true values of (grad y)? will be measured. 

Finally, it must be pointed out that the first two 
methods of calculating « lead consistently to results 
which are several orders of magnitude larger than those 
obtained from radio data.2® On the other hand, the 
third method yields values on the bases of Crain’s meas- 
urement?! of (grad 2)? which are too low. In view of 
the closing comments in the previous paragraph, to- 


gether with the possibility that Brunt’s estimate of € 


may be too large, this is not surprising. It remains to be 
explained, however, why the other techniques, espe- 
cially that involving E,(k) directly, predict values of 
scattering coefficient which are too large. . 


VI. CONCLUSION 


It has been shown that the Z* law which Villars and 
Weisskopf and Wheelon associate with “mixing-in- 
gradient” is untenable. It has also been shown that the 
L11/3 Jaw given by Batchelor and Silverman for isotropic 
mixing is equally appropriate for a large range of scales 
when mixing of an established gradient takes place and 
that, if the presence of a gradient has an effect, it must 
be to reduce the exponent, not increase it. This implies 
that the scattering coefficient should be proportional to 
\-3 and leaves unexplained the \ dependence which 
many of the radio data appear to indicate. That em- 
pirical spectra of refractive index fluctuations them- 
selves agree with the Z1/3 law has been demonstrated. 
The fact that the calculation of ¢ in accordance with the 
L"/8 law leads to vaiues which are in excess of those 
measured by radio means is also unexplained. All of 
these results, taken together, seem to indicate that 
scatter theory, as it is currently based on atmospheric 
turbulence, can provide at best an imcomplete descrip- 
tion of trans horizon propagation. 


APPENDIX | 


DERIVATION OF THE TIME RATE OF CHANGE 
OF THE SPECTRUM OF MEAN SQUARE 
y FLUCTUATIONS 


Defining the three-dimensional Fourier transforms of 


y, v, and 6éy and u;, a;, and éu; within the finite volume 
V by 


™ D. Brunt, “Physical and Dynamical Met @ i 
University Press, p. 286; 1941. : iets ber Sle Peed | 
a R. eg Sanne |” ata Mixing and Its Role in Radio 

ering,” Corne niversity, School of Elec. : 
Bae Reso ont y, School of Elec. Eng. Res. Rep. EE 

21 Private communication. 


1958 


Ay(k, t) = { ve, t) exp | —jk-x}dx 


Ay(k) = J ¥@ exp | —jk-x}dx and = A,(k) = [ a, t) exp Sera tes: 


5A,(k, ¢) = [ a, t) exp { —jk-x} dx 
Vv 


respectively, 

Ay(k, t) = Ay(k) + 64y(k, #) 
and 

A,(k, t) = A,(k) + 64,(k, 2). 
Now, if dy satisfies the conservation equation 


by ow by 0*6y 
rie ory ibs Soret At gore 
Ot Ox; Ox; OX,OX, 


(20) 


and if grad ¥ may be treated as constant throughout V, 
by substituting for y and u; in terms of their Fourier 
transforms one is led to 


dd Ay(k, ¢ oy 
[54y(k, 2)] aay e) ee 
ot Ox; 
J 
Ba, Be, ki Ak — k’,t)5Ay(K’, t)dk’ 
82° k-space ( 4 


— Dyk?5Ay(k, t). (21) 
Then, expanding 
al | 64y(&, 4) |*] 
Ot 
and making use of the assumed dynamic incompressi- 
bility of the air, 


SOT  oaladdesdbAy(k,)} 
al ot 


a 


e 8x3 


f fedettae 4ile = Hi, DdAGK, 9 
k—space 


— §Ay(k, t)A*(k — k’, 5 Ay*(k', )} dk! 


— 2Dyk?| 6 Ay(k, A) |?. (22) 
Taking averages throughout, 
5Ay(k, t) |? ay 
al] 84y(4, |? = —2@a{84,(k, )d4y*(k, 1} 
Ot Ox; 
- Oy(k, k/)dk’ — 2Dyk?| 5Ay(K, t)|*, (23) 


k—space 


where 


ki be y 
On(k, RY = — 9 {SAy*(k, DAK — K’, 18 Ay(K’, D) 


— 5Aj(k DARK = K’, 8 AY*(K,, OD}. 
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A,(k, 1) = J mc, t) exp fees 


(19) 


6A,(k,t) = i du;(x, ¢) exp { —jk-x}dx 
Vv 


If one now defines the spectral density of mean square 
y-fluctuations per unit volume, E,(k), such that 


[o) 1 Sete sy) eee 
f E,(k)dk = — f | dy(x, 2) |2dx, (24) 
0 Vy 
application of Parseval’s theorem yields 
k? SEE 
E,(k) = —— f | 8Ay(k, t) |2dQ(K). (25) 
8r3V Ar 


On integrating (23) over all directions of the vector 
wavenumber k and dividing both sides by 87° V, one has 


OL, (Rk) ow ce 
Se Ht P,(k, k')dk’ — 2Dyk?Ey(k), (26) 
dt Ox; 0 
where 
2 
I(t) = ——— [a SAKK DEAF D} a0) 
8rFV oJ ay 
and 
hp’? 
Pyle, Bt) =—— ff ul, K)aaK) a0). 
8r°V 4n VY 4 
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THE RELATION BETWEEN THE ONE- AND THREE- 
DIMENSIONAL SPECTRA OF SCALAR POINT 
FUNCTIONS 


Let 5n(x) be a scalar point function defined in a finite 
volume V and 5A,(k) be its Fourier transform, such that 


1 
n(x) = air 54,(k) exp {jk-x}dk, 


us k—space 


| 6A,(k) = f anc) exp {— jk-x} dx. (27) 


Define the covariance of 67, @nn (€), as 


donn(—) = sf antayintx + &)dx 


aa f | 84,(k) |? exp [7k-E}dk. (28) 
823 V k—space 
Writing 
2 
[PAF us) 
8r3V 
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one has 
il 
,(k) = —~ un(é) exp | —jk-E}dé, (29) 
81° é—space 
since @nn(~)—0 rapidly as | E| —«, Then 
1 ce) 
Vly 0 
and 
k(t) = Bf &4(4).d2(K). (30) 
An 


Now, along any diameter of V, of length M, define 


1 . 
8n'(1) = 6n(a, 0, 0) = ae 5A,(k) exp {jkixi}dk 


7 k—space 


1 fe) 
= —| 6A,}(ki) exp | jkiar} dha, (31) 
TY 2 


where 


1 ee) co 
5A,(k1) = —| f 5A, (k)dkodks 
Wa? Sco) ed 


M/2 
= f 6n1(%1) exp | —jkixy} day 


—M/2 


Then the covariance of 6,1, @nn!(&), is 


it M/2 
Pnn' (Ex) = —{ én(x1, 0, 0)dn(24 ap ei, 0, O)dxy 
MJ _wv}2 


1 ne} 
ee |S Ae 
are (h,) 


2 exp {jkiéitdk, (32) 


Writing 
6A,1(Ri) |? 
| Ants) |? ©,1(K), 
27M 
one has 
2 2} 
ai(bi) == [ dbunt(&) cos fasti}dér, (33) 
TJ 0 
since @nni(—&:)—0 rapidly as &—2*. Now, either 


dnni(é1) =Gnn(E1, 0, 0), or they converge stochastically. 
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But 
boalti 0,0) = f Bak) exp {rats} dk 
k—space 
= f ©,,(h1) cos | Rigs} dha, (34) 
0 
where 


©O,(R1) = 2f f ®,,(k)dkodk; hs 


2 -) 
=f nn (1, 0, 0) cos | Rigs} dé:. 
T 0 


Thus, ®,1(R:), the spectrum of data along a straight 
line, is equal to ©,(1), the one-dimensional contraction 
of the three-dimensional spectrum ©®,(k), or at least 
they also converge stochastically. 

Finally, therefore, if 6”(x) is isotropic, so that 
E,,(k) =41k?®,,(k) or ®,(k) = En(k)/ (4k? ], 


@,'(k1) = On(hi) = ai if 


E,(k) 
Ark? 


dkodk3 


ce Ce PatR 
=i  amkik = [ ©) ie (35) 
(ky| 27k? 7 area 3 
and 
d®,,1(k3) a dO@,,(k1) Le E,( | ki| ) 
dk, dky | 1 | 
Or 
dO), (k d®,1(k 
dk dk 
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illite Influence of Moisture in the Ground, Temperature and 
Terrain on Ground Wave Propagation in the VHF-Band* 


E Summary—Results are given from a general investigation re- 
_ garding the different factors that govern vhf ground wave propagation 
at low heights and short distances, i.e., mainly in the Sommerfeld 
_ region. In this work the influence on field strength from ground 
_ moisture and temperature, snow, topography, and vegetation, and 
_ also field strength variations along mixed land-water paths have been 
_ studied. The measurements refer to terrain typical for middle 
_ Sweden. 
% Suitable formulas for field strength calculations over smooth 
_ ground are developed. The moisture in the ground is found to play 
a principal role and a method is shown for calculation of changes in 
_ the path attenuation caused by known or predicted changes in the 
water content of ground. The influence of temperature is found to be 
~ small, even when the ground freezes, but a covering of snow may be 
of greater consequence. Characteristics of local variations and ex- 
“ponential attenuation of the field in a forest is demonstrated, the 
fluctuations being very large and the average exponential attenua- 
_tion factor mostly very small. Measurements along a mixed land to 
fresh-water path verifies Millington’s recovery effect. The field over 
- ahill is found to have a maximum value a certain distance below the 
top towards the transmitter. 


EE eee AS 


ie THEORY 


OR a smooth ground, the influence of the ground 
a Preonstnts on field strength is determined by the 
: complex dielectric constant (relative value) 
- & =€—j60oA =€(1—j tg 6). By introducing a “ground 
_ function” ¢={o-e® differently defined for horizontal and 
for vertical polarization, 


(1a) 


e i ? 
: 1 
~ Gos - 


? 
¢; — 1 


(1b) 


e general formulas for both kind of polarizations can be 


_ developed. The reflection coefficient at small grazing 
~ angles will be 

4 K = — (1— v9), (2) 
4 

if 

} yi K1, 

2 

2 


and the space wave E,, (vector sum of the direct and 
_ the ground reflected waves) 


E = 
3 ee = | wVF + jel, (3) 


J * Manuscript received by the PGAP, July 5, 1957. 
i + Res. Inst. Natl. Def., Stockholm, Sweden. 
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where 
47H, 
= ——— <I] 


2) 
AR 


f,and Hz are antennaheights, k = distance, Ey = +/30PD, 
P=transmitted power, D=directivity of transmitting 
antenna. 

The surface wave is 


su = a8 ? pPreaie, (4) 
Rp 
where the numerical distance p is given by 
p = tR/(Xo). (5) 
The height-gain function at small heights will be 
beak eae ly (6) 
AVF 


At vhf, over land and fresh water, tg 6 is negligible 
and the above formulas then give: 
Horizontal polarization: 


EG i 
yea (7a) 


a 
Esp = En 4/14 :G, (7b) 
aes 
where 
4rH,H» 
En =s Eo UR ) (8) 
r r 2 
Gi + ) : (9) 
ee 27H» 


At horizontal polarization the surface wave can be 
neglected, (7a) thus giving the total field strength. 

At vertical polarization it is suitable to calculate the 
total field strength from the surface wave and the 
height-gain functions, 


Etot ae Eyu:¢(H1)¢(H2), (10) 
fiat ee eimig (11) 
jae ares 
Sat ae 
igri Vike (= \, (12) 
2 
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if 


3r\e < R < 12000-A¥? (AX and Rin™m). (13) 


From (7)—(13) it can be seen that the field strength 
at low heights in the Sommerfeld zone at vhf obeys an 
inverse distance squared law. 

Therefore the influence of a change in e€ is equivalent 
in this case to a change of the antenna efficiencies; this 
indicates that only the properties of the ground near the 
antennas affect the received signal. 

With this in mind, and in view of the law of reciproc- 
ity, the dielectric constant at both ends of the trans- 
mission path should have the same influence on the 
field strength. For the surface wave at vertical polari- 
zation over homogenous ground this influence is, ac- 
cording to (11), represented by the factor e//e—1 at 
each end. If the dielectric constant has the value & at 
one end and e& at the other end of the path the surface 
wave therefore will be 


Eo rv €1€2 
Re Tv V/ (ex as 1) (€2 = 1) 


Euu(e1, €2) = (14) 


This has been shown more rigorously by Millington. 


THE INFLUENCE OF MOISTURE IN THE GROUND 


Vhf field strength recordings at vertical polarization 
over smooth ground have shown variations up to 14 db 
over the same paths at different seasons. The measured 
field strengths correspond to e-values from about 3 to 
30, as calculated from (10)—(12). 

Investigations have shown that the value of the ef- 
fective dielectric constant is determined mainly by the 
water content of the ground and is relatively independ- 
ent of the type of ground. 

Very dry earth has effective e-values between 2.5 and 
3, and for wet earth € is approximately proportional to 
the percentage water content, w per cent. Therefore, 
in vhf field strength calculations the empirical relation 


€ = 0.78-w + 2.5 (15) 


may be used. 

In most parts of Sweden the ground humidity has an 
average annual variation with maxima in the fall and in 
the spring and a low minimum in the summer. The range 
of vhf communication sets with low antennas varies 
correspondingly up to a ratio of more than 2:1. In 
military applications it is often important to know the 
range to be expected in a particular case. It is, however, 
difficult to determine the value of w in (15), which 
would enable the calculation of the range. Therefore, an 
apparatus has been developed for measurement of the 
effective e, based on measuring the slope of the E vector 
from an adjacent test transmitter with vertically polar- 
ized antenna. 

In some cases it is only required to know how the 
propagation conditions within a certain area have 
changed from one time to another. In solving this prob- 
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lem the following expression has been found to yield a 


i 
es 


useful measure of the variations in the water content of 


the ground, 


w=C+) (p- 9), (16) 


where 


p is the precipitation (mm), 
e is the evaporation (mm), 
C is the water content at the start of observation. 


The quantities p and e are relatively easy to deter- 
mine. A linear relation between € and w, as expressed 
by (15), is assumed and from measurements at different 
times the e-w-curve can be drawn. The use of this rela- 
tion is explained with reference to Fig. 1, showing a 
sketch of a nomogram between w, ¢, A, distance R be- 
tween transmitting and receiving antennas and the 
field strength E at the receiver relative to the free space 
level at the same distance. At a certain time the field 
strength E; is measured, which gives the values ¢ and 
w,. If thereafter the water content would decrease with 
Aw mm a field strength EZ, may be predicted at the re- 
ceiver. 

This method has proved to give good results in time 
intervals of a few months. 


Fig. 1—Sketch of a nomogram for field strength calculations. 


THE INFLUENCE OF TEMPERATURE AND SNOW 


The temperature coefficient of the dielectric constant 
of ground is very small, of the order of —0.005 per °C, 
so no direct influence of temperature on vhf propaga- 
tion may be expected. To check this, an analysis has 
been made of field strength measurements from differ- 
ent seasons, and no correlation with ground or air tem- 
perature was found. 

However, an indirect influence may exist in so far as 
the ground moisture content and thus e might decrease 
during an extended warm period. 

When the ground becomes frozen or covered with 
snow, the dielectric constant will change with an amount 
which depends on the original e-value. The higher this 
is, the greater will be the change even in relative meas- 
ure. From propagation measurements at 50 mc over 
ground with e=15 originally, the curves in Fig. 2 have 
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i 20 60 | 100 cm 
Depth of ice, snow ond frost. 


Fig. 2—Influence of ice, snow and frost in the ground on the effective 
dielectric constant, measured at 50 mc. Curve a=Ice on a lake, 
Curve b=Frost without snow, Curve c=Snow without frost, 
and Curve d=Snow and frost. 
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Fig. 3—Field strength variations, (a) at an open field, (b) in a forest. 
Frequency 44 mc. 


been calculated. It is seen that a covering of snow 


causes a much greater decrease in € and the field strength 


_ than does the freezing of ground. 


FIELD STRENGTH VARIATIONS AT OPEN FIELDS 
AND IN Woops 


A great number of field strength recordings over dif- 
ferent types of terrain have been made within the vhf 
range both at vertical and horizontal polarization. In 
carrying out these tests the receiving and recording ap- 
paratus were placed at a fixed position while the trans- 
mitter was moved at a constant speed. The marked dif- 
ference between open fields and woods can be seen from 
Figs. 3(a) and 3(b) respectively, which show typical 
recordings at 44 mc and vertical polarization. An analy- 
sis of such curves measured over approximately plain 
ground has shown that the field strength median over 
a 100 m distance obeys an inverse distance squared law, 
as expected. 

The amplitude of the local variations depends on the 
vegetation, the frequency and the polarization. 

The 10 per cent and 90 per cent levels referred to +50 
m distance from a fixed point have been determined. 
Fig. 4 shows the difference in db between these two 
levels in a wood and at a field as a function of frequency. 
It is seen that at higher frequencies even horizontally 
polarized waves are much affected by the trees. 

A relief picture of the field around a tree is shown in 
Fig. 5, which refers to 35 me and vertical polarization. 
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Fig. 4—Frequency dependence of local field variations within a 
distance of 100 m in a forest and at an open field. 
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Fig. 5—The field around a tree, 35 mc; vertical polarization. 


A variation of 11 db occurs within a distance of 1/3 
wavelength from the tree. 

These measurements were made in a pine-forest of a 
height up to about 20 m. 


ATTENUATION IN A FOREST 


In certain cases an exponential distance variation of 
field strength has been observed over short paths in a 
forest, t.e., E proportional to d~?-e~ *. 

For frequencies below 100 mc the exponential factor is 
insignificant in high pine-forest with little underwood. 
In a forest with dense leafy underwood attenuation 
factors at vertical polarization as high as 0.02 db/m at 
30 mc and 0.04 db/m at 100 mc have been measured at 
small heights and over short paths. For longer paths 
(but in the Sommerfeld region) the attenuation factor 
a decreases. It is thus evident, that only the trees within 
a certain distance from the antennas affect the field 
strength appreciably at the receiving antenna. This is 
similar to the influence of the dielectric constant of 
ground. 


FIELD STRENGTH VARIATION NEAR A SHORE 


Fig. 6 shows a measured field strength curve a at a 
frequency of 44 mc along a mixed land to fresh-water 
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Fig. 6—Curve a= Measured field strength curve for a mixed land- 
water path, Curve b=Calculated curve for e=8, and Curve c 
= Calculated curve for mixed path, « =8, ee =80. Antenna heights 
1.5 m (fixed antenna) and 4.5 m (movable antenna). Frequency 
44 mc; vertical polarization. 


path. The receiving antenna was located over land 400 

m from the shore at a height H;=1.5 m. The transmit- 

ting antenna was moved at a constant height H,=4.5 m. 
The dashed curve 6 is calculated from (10)—(12): 


Ey — Eyu(e1)o(A1, €1)¢(H, €1) (17) 


This curve corresponds to € = 8 over land. The dashed 
curve c is calculated from (10), (12), and (14): 


Fy = Eales, €2)6(Ai, €1)¢(H2, €2) 


‘ where €,= 80. 
The transition between the two curves occurs over a 
distance AR, which measured in wavelengths closely 
corresponds to the quantity 


AR vs €1€2 
Nitin (ieee 1) 


1.e., the proportionality factor of the surface wave over 
the mixed path according to (14). 

The measured increase in field strength between the 
minimum and maximum points of the curve (Milling- 
ton’s “recovery effect”) is 3 db which well agrees with 
calculations from the above formulas. 

If the antenna over land were at greater distance from 
the shore (7.e., R>>AR) the field strength gain when pass- 


(18) 


(19) 
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Fig. 7—Curve a= Measured field strength curve over a hill, 35 me; 
vertical polarization, Curve b and Curve c=Calculated curves 
for smooth ground, uniform and mixed respectively, and Curve 
d=Path profile. 


ing from land to water at a low height (7.e., if (2) =1) 
according to (17) and (18) would be 


€2 ér i! 


€1 €éo-+ 1 


which amounts to 9.5 db for e, =8, €2 = 80. 


VARIATION OF FIELD STRENGTH OVER A HILL 


In Fig. 7 measured field strength curve a is shown 
over a hill adjoining a lake. The dashed curves a and b 
are calculated as the corresponding curves in Fig. 6. 

The maximum increase in field strength above the 
inverse distance squared value occurs a short distance 
(about 5 to 6 wavelengths) before the top of the hill is 
reached. In the region behind the hill the field strength 
is determined by a combination of diffraction and the 
recovery effect over the water surface. 
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Distance Dependence, Fading Characteristics and Pulse 
Distortion of 3000-MC Trans-Horizon Signals* 


B. JOSEPHSON} anv G. CARLSON} 


Summary—Statistical analysis of signal recordings at 10 cm 


_ wavelength has given the following results. 


The distance dependence of the half-hour signal median F,, over 
a sea path corresponds to a scattering parameter which is approxi- 
mately inversely proportional to the height. Between F,, and the sur- 


_ face value of M a correlation coefficient of 0.62 is found, and between 


F,, and the dew point 0.55. A unit increase of M seems to be followed 
by a signal increase between the limits 0.3 and 0.8 db. The fading 
range, as defined between the 10 per cent and 90 per cent of time 
levels, is 12-14 db, and the fading frequency mostly about 2 c /sfora 
300 km path. 

Over a 260 km land path the pulse broadening was normally 
0.1-0.2 4s which corresponds to a scattering parameter varying in- 


_ versely as the first to second power of the height. A severe pulse 


distortion caused by selective fading is found to be very common. 
Some other propagation characteristics are also discussed. 


I. GENERAL 
(Cex waveene recordings of field strengths at 10 


cm wavelength have been made during periods 

of several weeks at different seasons of the year. 

Paths from 200 to 300 km lengths, both over land and 

Over sea, have been used. The measurements were 

carried out with pulsed signals, and apart from the study 

of field strength variations, the shape of received pulses 

has been observed. Vertical polarization was used 
throughout. 

Geographically the investigation refers to middle and 

south Sweden including the east and west coast waters. 


II. Test EQUIPMENT 


For the sending end, two radar transmitters have 
been used, one with 200 kw and the other with 120 kw 


_ pulse power. At the earlier tests the pulse width was 


2 us, later 0.8 us. The bigger transmitter had a rec- 
tangular shaped antenna with a directivity of 37 db. 
The other transmitter was equipped with a centerfed 
paraboloid reflector of 300 cm diameter having a direc- 
tivity of 37 db. 

At the receiving end the last mentioned antenna type 
and also a 130 cm paraboloid (29 db directivity) were 
used. 

In the first experiments the receiver had a noise 
factor of 12 db and in the later 8 db. This was achieved 
by using cascodes according to Wallman and balanced 
mixers. The amplitude range was about 50 db. The re- 
ceived signals have been recorded on an Esterline 
Angus recorder for median values and during short in- 
tervals on a Brush recorder for rapid fading analysis. 
Pulse distortion has been studied on a Tektronix type 
514 AD and oscilloscope. 

* Manuscript received by the PGAP, July 5, 1957; revised manu- 


script received, November 19, 1957. 
+ Res. Inst. Natl. Def., Stockholm, Sweden. 


Statistical analysis of the recordings was carried out 
by the aid of a semiautomatic curve-analyzer, which 
has been specially developed for propagation studies. 


III. REsuLtTs 


In the analysis of the measured data, the propagation 
factor F is-introduced, by which the received power 
P, is referred to the calculated free space power level 
Py, given by 


(1) 


where 


P,=transmitted power, 
G,=directivity of transmitting antenna, 
A,=absorption cross section of receiving antenna 
(m*), 
d= path distance (m)._ 


The propagation factor F is then defined by 


P, = F2- Pyy. (2) 


A. Medians of the Propagation Factor 


From the recordings median values of F during half- 
hour intervals have been determined. In Fig. 1 the dis- 
tribution of such medians, called F,,, are shown for four 
different cases. Curves A and B refer to measurements 
over sea (south Sweden) at distances of 190 km and 290 
km respectively along same line from the transmitter, 
carried out continuously during 5 weeks in the months 
of October and November. Curves C and D refer to a 
260 km path over land (middle Sweden) and represent 
continuous recordings during December and January 
respectively. Curve B shows signs of two different dis- 
tributions. 

It is not possible to decide from these measurements 
to which extent the difference between the curve B and 
the curves C and D is due to a characteristic difference 
between sea and land paths. The test paths in both cases 
are almost of the same length, but obstacles at the two 
terminals on the land path will give about 4 db extra 
propagation loss. Furthermore the surface reflection 
on the sea path may give some interference gain. 
Weather conditions were quite different, such that the 
air mass in the scattering region was colder and had 
lower relative humidity in the C- and D-cases than in 
the B-case. 

In Fig. 2 the distance dependence of the F,,-factor 
over the two sea paths is elucidated. A distance power 
law F,2~d-" is assumed, and the curve shows the time 
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Curve A=190 km sea path October-November, 1954 
Curve B =290 km sea path October-November, 1954 


Middle Sweden 


Curve C =260 km land path December, 1956 
Curve D =260 km land path January, 1957 


Fig. 1—Distribution of 3000 mc half-hour median propagation factor 
(log-normal-distribution). 


Percent! of holl-hour intervals ordinole (3 exceeded. 


Fig. 2—Distribution of 3000 mc half-hour median distance 
dependence (log-normal-distribution). 


distribution of m. It is seen from this curve, that for 
shorter periods (approximately 10 per cent of total time) 
the F,,-factor for the remote receiver has been larger 
than that for the nearer receiver. An inversion layer 
at a height and a location suitable for the longer path is 
a possible explanation. 


B. Fading Characteristics 


Statistical analysis of short time intervals of the 
scattered signal shows a Rayleigh distribution of the 
rapid fading (Fig. 3, lower curve). For longer intervals 
the distribution changes to normal (upper curve). The 
fading range, as defined between the 10 per cent and 90 
per cent levels, amounts to 12—14 db, and is relatively 
independent of the momentary signal median, provided 
this does not exceed the long time median too much. 

The fading frequency (number of times per second 
the median is passed in one direction) is normally about 
2 c/s at a 300 km path and is approximately propor- 
tional to the path length. 

It is interesting to know how often the signal has 
stayed below a given level for a given time. This has 
been evaluated for a number of cases, an example of 
which is shown in Fig. 4. Such analysis is conveniently 
carried out on the above-mentioned curve analyzer. 
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Fig. 3—Distribution of instantaneous 3000 mc signal level 
(log-Rayleigh-distribution). 
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Fig. 4—Number vs duration of fast fades at 3000 mc on a 
290 km path. 


C. Pulse Distortion 


Pulse distortion has been studied over the 260 km 
land path. The shape of the transmitted pulse is shown 
in Fig. 5(a), pulse length being 0.8 us at midheight. The 
received pulses were photographed with 1/25 sec ex- 
posure time, thus giving the average shape of 20 pulses. 
Normally a pulse broadening of 0.1 to 0.2 us was ob- 
served, as in Fig. 5(b). Also Figs. 5(c) and 5(d) show 
about 0.2 us broadening, however in combination with a 
severe distortion caused by selective fading. This type 
of distortion is very frequent, and in general the shape 
of the received pulse is continuously changing within 
about 0.2 us increased total width. 

Not so often a much more pronounced pulse broaden- 
ing, up to 0.8 ws, has been observed. Fig. 5(e) shows a 
pulse of 1.2 us width, 7.e., 0.4 us broadening. 


D. Signal Dependence on Meteorological Factors 


From meteorological observations near ground at 
both ends of the 290 km sea path the dew point tp °C 
and the modified refractive index M were calculated. 
Mean values of tp and AM (M=320+AM) for both 
ends of the path at the times 7 Am and 7 pm each day 
are shown for a 15 days period in Fig. 6. In the same 
diagram half-hour medians of the propagation factor 
at times 7 AM and 7 pM are plotted. 

From these curves it is evident that a correlation 
exists between the propagation factor median and both 
the M value and the dew point (or the partial water 
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Fig. 5—Pulse distortion over 260 km land path. (a) Transmitted 
path, (b)—(e) Received pulses. 


y 


_ vapor pressure e = 6.1-exp (0.072¢p)mb). The correlation 


coefficient has been calculated to 0.62 between F,, and 
AM, and 0.55 between F,, and tp. The correlation anal- 


_ ysis for F,, and M gives the two regression lines a slope 


of 0.3 db and 0.8 db per unit change of M. 

It is recognized that these figures are very uncertain 
in view of the short period analyzed. 

The measurements over the 260 km land path which 
now have been extended over a much longer period will 
be statistically analyzed in the same way. This will 
indicate to which extent the deviations between the 
curves C and D for the land path and curves A and B 
for the sea paths in Fig. 1 may be correlated with aver- 
age deviations in refractive index or dew point near 
ground. 


E. Various Observations and Calculations 


As mentioned in Section III-A relatively high signals 
over the longer sea path are sometimes associated with 
low signals over the shorter path. At such instances the 
fading frequency over the longer path is much lower 
than normal, often by a factor of 10, while at the same 
time a severe pulse distortion occurs of a type indicat- 
ing multipath propagation. 

An opposite situation is also frequently observed, 
where a high signal level, up to the free space value over 
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Fig. 6—Propagation factor and simultaneous surface values of 
modified refractive index and dew point. (290 km sea path, 
3000 mc.) 


the shorter sea path is concurring with a very low signal 
at the remote receiver. This is associated with duct 
propagation over the shorter path. A narrow peninsula 
of about 70 m height on which the receiver for the 190 
km path is situated, prevents the formation of a con- 
tinuous surface duct along the longer path. The high 
attenuation over this path may be explained by the 
change of shape of the lowest lobe caused by the duct. 
It is also possible that the intensity of the turbulence is 
smaller at these weather conditions. 

The signal received by the scattering mechanism is 
calculated from the quantity (scattering coefficient) 


a = S(h)-f(9), (3) 


_where @ is the scattering angle and S(h) = (Ae/e)?/s is 


the scattering parameter (s=scale of turbulence). The 
scattering parameter as a function of the height h may 
be assumed to obey a power law, 


(4) 


Different p values give different distance dependence of 
the signal and different pulse widening. 

Fig. 2 shows that the half-hour median of the propa- 
gation factor squared during 50 per cent of time has 
been inversely proportional to approximately the third 
power of distance. This corresponds to a p value in (4) 
of slightly more than 1. 

As mentioned in Section III-C, a pulse broadening of 
0.1 to 0.2 us was normally observed over the 260 km 
land path. This corresponds to p values from 2 to 1. 

For the type of distortion shown in Fig. 5(e), the pulse 
widening may be estimated to 0.4 us, indicating a height 
difference between the propagation paths of about 3000 
m. An inversion or a turbulent layer at a height of about 
4500 m would therefore be a possible explanation of this 
distortion. 
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Some Microwave Propagation Experiences from a 
“Just-Below-Horizon” Path* 


B. JOSEPHSON} anp F. EKLUNDt 


Summary—Based on comparisons between field strength records 
at 2300 mc and simultaneously measured M curves it is concluded 
that the most important propagation disturbances on the test path are 
caused by reflections from discontinuities in the refractive index 
profile. The occurrence probability for such discontinuities is shown 
to have a marked seasonal and diurnal period with a maximum in 
the summer and around sunrise. 


I. Test PatH AND MEASUREMENTS 


HE path profile is shown in Fig. 1 for a modified 
eee radius a,=4a/3. The line between trans- 

mitter and receiver is broken mainly by two ridges 
covered with wood. 

Field strength recordings over this path have been 
made almost continuously since the spring of 1955 ata 
frequency of 2300 mc. During some periods, meteoro- 
logical soundings with the aid of thermistor sonds 
carried by a captive balloon were made at a point on 
the path up to a height of about 300 m. These measure- 
ments were made mostly when the prevailing weather 
situation was favorable for air stratification at low 
height. 
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Fig. 1—Path profile. Earth radius 4/3-6370 km. 


II. THEORETICAL CONSIDERATIONS 
A, Diffraction 


The propagation factor F (the ratio between actual 
and calculated free space field strengths at receiver) 
corresponding to diffraction signals has been calculated 
by two methods: 1) Diffraction over an equivalent 
sharp edge (calculation carried out by aid of Cornu’s 
spiral). This gives F equal to —18.5 db. As seen later 
this value is far too high. 2) Diffraction over four 
spherical surfaces determined by the antenna heights 
over surrounding terrain and the height and position 
of the two ridges.! The result is F equal to — 35 db. 


* Manuscript received by the PGAP, July 5, 1957. 

+ Res. Inst. Natl. Def., Stockholm, Sweden. 

1K. A. Norton, P. L. Rice, and L. E. Vogler, “The use of angular 
distance in estimating transmission loss and fading range for propaga- 
tion through a turbulent atmosphere over irregular terrain,” Proc. 
IRE, vol. 43, pp. 1488-1526; October, 1955. 


B. Scattering 


The propagation factor corresponding to the scattered 


signal median has been calculated according to the 


Booker and Gordon theory. The scattering parameter 
was assumed to be independent of height, which is justi- 
fied over this short path. An F value of about —50 db 
was found. 


Thus the scattered signal is negligible compared with — 


the diffraction signal. 


C. Tropospheric Reflections 


A plane horizontal discontinuity in the atmosphere : 


is assumed where the refractive index changes from 


n to n(1—An). At nearly grazing angles y the reflection 
coefficients for vertically and horizontally polarized — 


waves will be 
1 — 2An — V/1 — 2An/yp? 
1 — 2An + V/1 — 2An/yp? 

1— /1 — 2An/¥? 
EAT = Aa ee 
This gives for both kinds of polarizations the relation 

2K 
fuk Kes 
where W is measured in milliradians and AM =An- 108, 

In view of the fact that no sharp discontinuities exist 
in the atmosphere the above formulas may yield too 
high values of the reflection coefficient. From a rigorous 
treatment of this problem by Millington? it may be con- 


cluded, that (1) be fairly accurate if the change AM 
take place in a height Af fulfilling the inequality 


Ah < d/y. 


vest Paar 


(1) 


AM = y- (2) 


(3) 


When the field at the receiver is composed of dif- 
fracted and reflected waves of about the same magni- 
tude this will be observed as interference fading and an 
about 3 db increase of signal median over the normal 
diffraction level. This case would occur if the reflection 
coefficient is equal to the diffraction propagation factor 
as defined in Section II-A. Thus a diffracted signal of 
case 2 in Section II-A corresponds to K = 107): =0.018. 

The necessary values of AM, calculated from (2) to 
give this value of K over the test path (Fig. 1) are: 


150 200 400 900m 
1 Sutera Oe 2G 


Height of discontinuity 
AM 


2 G. Millington, “The reflection coefficient of a linearly graded 
layer,” Marconi Rev., pp. 140-151; October-December, 1949., 


~ index have been analyzed with a view of distinguishing 
the influence of the different propagation mechanisms 
just discussed and to show the probability of their oc- 
currence during different parts of the year and of the 


day. 
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D. Refraction 


From the path profile it may be calculated that a 


1 direct ray can reach the receiver by refraction if the 
_ curve has a gradient 


dM/dh < — 0.13 m- 
up to a height of about 150 m. With M gradients in be- 


_ tween the normal atmosphere value +0.118 m-! and 


the above figure, the diffracted signal will be between 


4 its normal value and approximately the free space level. 


III. Comparison BETWEEN MEASURED AND 
CALCULATED DaTA 


The measurements of field strength and refractive 


A. Normal Atmosphere 


When the meteorological soundings have shown a 
normal M curve, the received signal has always been 
very stable. The propagation factor was then —37 db 
in good agreement with the calculated diffraction level 
over spherical surfaces according to case 2 in Section 


II-A. 


B. Low Level Stratification 


Many soundings have shown fairly sharp low level 
discontinuities in the M curve as demonstrated in Fig. 
2 together with the simultaneously taken field strength 
records. 

Fig. 2(a) shows a rather small discontinuity. The sig- 
nal median is somewhat higher than the normal dif- 
fraction level. The fading is rapid and has an amplitude 


of about 10 db. If the assumption is made that the 


fluctuations depend on varying reflected signal inter- 
fering with a steady diffraction signal, the corresponding 
reflection coefficient of the layer at 120 m height will be 


- approximately 0.02 as calculated from the record. 


When calculated from (1) and the M curve, K is ap- 


3 proximately 0.05. 


De ted «ks pl 


In Fig. 2(b) a record is shown from an occasion, at 
which the signal during several hours has exceeded the 
free space level (i.e., F>1). Extremely deep fading is 
also noticed. Calculations from the M curve shown in 
Fig. 2 and other M curves measured during the same 
period given reflection coefficients varying between 0.2 


and about 1, mainly because of changes in the height 
_ of the layer. The condition (3) is not quite met in this 


case, but the error will be small. The concave shape of 
the reflecting surface may possibly further enhance the 
reflected wave. 

Thus it seems probable that the high signals observed 
in this case and the fading in the first case have been 
caused by reflections. In fact, at every occasion when 
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Fig. 3—Distribution curves for half-hour medians. Curve a 
=month of July, Curve b=month of December. 


simultaneous meteorological measurements and signal 
recordings are available over this path, and the signal 
has exceeded the normal diffraction level more than 
about 5 db, this could be explained as caused by low 
level tropospheric reflections on the same indications 
as is exemplified above. On the other hand, a refractive 
index gradient such that a direct wave would be re- 
fracted from transmitter to receiver has never been 
observed. 


IV. STATISTICAL ANALYSIS 


A. Distribution of Signal Medians 


Characteristic distribution curves for half-hour me- 
dians of the propagation factor, F,, are shown in Fig. 3 
for one summer and one winter period. One finds that 
the winter offers much more stable propagation condi- 
tions than does the summer. 
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The curve for the summer period seems to be com- 
posed of two different distributions, one connected to 
diffraction propagation with a median of —37 db, and 
the other influenced by low level tropospheric reflec- 
tions. The transition between the two parts takes place 
at a level around 5 db over the median. The left part 
of the curve would by extrapolation give a median be- 
tween —50 and —60 db, and it is interesting to compare 
this value with the value —50 db calculated for scatter- 
ing mechanism. 


B. Occurrence of M-curve Discontinuities 


A minimum signal increase of 5 db over the median 
value is chosen as a criterion on existence of a reflecting 
tropospheric layer over the test path. Based hereupon 
the signal recordings during 18 months have been an- 
alyzed to give seasonal variations of the occurrence 
probability for such discontinuities and the result is 
shown in Fig. 4. During this summer the layers exist 
about 10 per cent of the time, during winter only 1 or 2 
per cent of time. 

Fig. 5 demonstrates the diurnal dependence of the 


probability during a summer period (months of July and - 


August, 1955), which is seen to be very marked. Reflect- 
ing discontinuities have never been found around noon. 
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Fig. 4—Occurrence probability of low level reflecting discontinuities. 
5—Measurements 1955, 6—Measurements 1956, 7—Measure- 
ments 1957. 
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Transients in Conducting Media* 
PAUL I. RICHARDS+ 


Summary—The exact fields generated by transient electric and 
magnetic dipoles in an infinite conducting medium are derived. 
Specialization to brief, approximately triangular pulses then brings 
out the salient practical features of such signals. The peak signals 
are attenuated as r~ or r“‘ rather than exponentially, and, in a prac- 
tical sense, they travel very slowly (mean effective velocity about 
equal to sonic speed for a range of 1 km in sea water). 


INTRODUCTION 


HE electric field excited by a pulsed magnetic 

dipole, or current loop, in an infinite conducting 

medium has been obtained in explicit form by 
Wait.1 The purpose of the present paper is twofold: 1) 
to extend Wait’s results by deriving explicit formulas 
for the magnetic field and for the fields excited by an 
electric dipole; 2) to discuss these solutions and point 
out some practical aspects of transient electromagnetic 
signals in sea water. 


* Manuscript received by the PGAP, July 6, 1957; revised manu- 
script received November 8, 1957. This study was supported by the 
Office of Naval Research. 

+ Technical Operations, Inc., Burlington, Mass. 

J. R. Wait, “A transient magnetic dipole source in a dissipative 
medium.” J. Appl. Phys., vol. 24, pp. 341-343; March, 1953. 


These practical features are, first, that even the peak 
values of wide-band signals are not exponentially at- 
tenuated but are attenuated according to an inverse 
power of distance. This, of course, is in sharp contrast 
to the attenuation of sinewave signals and the distinc- 
tion can be of great importance when considering the 
weak fields at “large” distances (100 meters or more in 
sea water). Secondly, the effective signal velocity is very 
slow and decreases with range to values as low as sonic 
velocities at distances of the order of 1 km in sea water. 


METHODS AND NOTATION 


For our purposes, it is expedient to write Maxwell’s 
equations for a homogeneous medium in the special 
form (mks units) 
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a=r/20T, m= 1/27 (2) 


then 
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Va Ms E = v0B/dr. 


Assume that all fields vanish (“almost everywhere”) 
when ¢=0 and take Laplace transforms with respect to 
the variable 7. This merely replaces 0/07 with » and 
the resulting equations may be solved by taking any 
solution of 


WAN Sas Serle eid (3) 
End setting either 

E = — poZ, B=VvV.XZ (4) 
or 

pie ta. Piet Vie OZ: (5) 
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MAGNETIC DIPOLE SOURCE 


If the source is a distribution of currents confined to 
“a very small region of space, its magnetic dipole moment 


is defined by 
fff x JdV 
5 iP ok TONE. 


To determine the fields generated by such a source, it is 
simplest to observe that the vector potential for a static 
magnetic dipole? is A=(u/47)MXQV(1/r) and the mag- 
netic field (in the midplane) is thus proportional to 
pM /4rr*. The solution of (3) and (4) which will give 
the same result at very small distances for arbitrary 
ptt) is 


M(t) = (6) 


Zn(P) “s 4 (20T)-*m(p) X Vala! exp (—av/p(p+2))] (7) 


where m() is the Laplace transform of M(t). Substitut- 
ing (7) into (4), considerable manipulation gives 


E= (u/8arTr?)ro x MF (t) (8) 
B = (u/4rr*)ro X (to X Mo)Gi(t) 
+ (u/4rr*) [3(ro-Mo)ro mo]|G2(2) (9) 


where subscript 0 denotes a unit vector and where, in 
terms of the normalized variables defined in (2), the 


2See, for example, J. A. Stratton, “Electromagnetic Theory,” 
McGraw-Hill Book Co., Inc., New York, N. Y., 1st ed., p. 237; 1941. 
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where Laplace transforms of Gi, Go, and F are 
T=¢/o, v=1/Ven (1) g(p) = a°p(p + 2) exp | —alp(p + 2)]!/2} m(p) 

are the characteristic time constant and characteristic g2(p) = {1 + allo + 2)]'?} 
velocity of the medium. Thus if we introduce normalized ex 5 

’ 4 2) ]1I2 
space and time variables defined by ona se 

Sp) = pgx(p). (10) 


The Laplace inversion of g:(p) is readily obtained from 
standard references.* The other two functions, go and 9 
were essentially inverted by Wait,! whose results can 
be stated as the following Laplace Transform pair 

£{ b%(2 — b2)e~7,[(2 — 6212] U4 — b)} 

= {1+ d[p(p + 2)]*/*} exp {—d[p(p + 2)]¥?} 

— {1+ b(p + 1) + 38°} exp [—d(p + 1)] (11) 
which may also be checked directly by evaluating the 
double integral of the pair (14) (with »=2).4 In (11), 
U(x) is the unit jump function defined by 

OQpit.ae <0) 
nar , Nie 
ile ak, area ae) 


and the now standard notation, [,,(Z) =i-"J,(iZ), de- 
notes a Bessel function of imaginary argument. We shall 
have frequent use for the following “kernal” functions. 


K,(@, 7) = aer7(r? — a?)—" 27, | (7? — a?)1/2], (12) 


With the relation (11) all of the functions in (10) may 
be inverted: 


Gi(r) = U(r — a)a®e*[M'"(r — a) + 2M"(r — a)] 
+ U(r -— of K,(a, u)[M’(r — u) 
+ 2M'(r — u)|du (13) 
G2(r) = U(r — a)[ae*M"(r — a) 
+ (1+ a+ 4a2)e*M(r — a)] 
+ U(r -— df Kea u)M(r — u)du (14) 
F(r) = U(r — a)[ae*M""(r — a) 
+ (1 + a+ $0°)M'(r — a)] 
+ U(r - a) J Kalo u)M'(r — u)du (15) 


where primes denote differentiation with respect to the 
argument. 

Eqs. (8), (9) and (13)-(15) constitute a complete 
solution (in our normalized units) for the magnetic 
dipole source. If we set M(t) =6(t), (8) and (15) become 
Wait’s solution! for the electric field (as simplified by 
recurrence relations for the modified Bessel functions), 
while (9), (13), and (14) constitute an explicit expression 


Magnus, Oberhettinger, and Tricomi, 


3 For example, Erdélyi, 
McGraw-Hill Book Co., Inc., New 


“Tables of Integral Transforms,” 
York, N. Y.; 1954. 
4 Tbhid., p. 201. 
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for the magnetic field. As will be seen below, Gi(r) 
represents most of the transient magnetic field while G2 
contains the dc component if there is any. The latter 
fact is evident from the character of the vector factors in 
(9) and may also be proved (using limp+o bf (p) =lim, +2 
F(r), if the latter exists) directly from (10), which shows 
that 


Gi(r) ~ 2a?M’(r) 
G(r) ~ M(r) 
F(r) ~ M'(r) 


aS eo 


The practical aspects of these various results will be dis- 
cussed below. 


ELEcTRIC DIPOLE 


With a small charge distribution for a source, the 
dipole moment may be defined as 


po= fff ear 


For a static electric dipole, the scalar potential® is 
(—1/47e)P-V(1/r) and the electric field (in the mid- 
plane) is proportional to P/4mer*; except for the factor 
(1/eu) this field is the same as the magnetic field from a 
magnetic dipole. Comparing (4) with (5), we see that 
the correct fields will be obtained by using (5) instead of 
(4) and the following modification of (7) 


(16) 


1 —SSSSSSS_— 
2. = =~ (20T)*a(0) X V.[a-! exp (—av/p(p + 2))] (17) 


where q(p) is the Laplace transform of P(t). Substitut- 
ing (17) into (5), considerable manipulation gives the 
results, 


= (u/8nTr?) Po X roR(t) (18) 
E = (1/4mer®)ro X (ro X Po)Si(é) 
+ (1/42rer®)[3ro(ro-Po) — PolS2(t) —(19) 


where as before subscript, 0, denotes a unit vector and 
the time functions have the Laplace Transforms, 


si(p) = a°p?q(p) exp { —alp(p + 2)]!?} 
so(p) = p(p + 2)*g(p) {1 + alp(p + 2)]1} 
-exp {—a[p(p + 2)]/2} 

(p + 2)50(p). 


ene if e’"P'(u)du 
0 


P(r) — dew f e“P(u)du 
0 


r(p) = 
| If we define 


II(r) 


(20) 


then the functions 5;, s2, and 7 may be inverted as before, 
and the following results are obtained. 


5 J. A. Stratton, op. cit., p. 175. 
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Si(r) = U(r — a)a®e* P(r — a) 


LO (aoa) he K,(a, u)P’(r — u)du (21) 


So(r) = U(r = a)[(1 + a + $0”)e“Il(r — a) 
+ ae~Il'(r — a)] 
+ U(r — a) f Kxl6, u)II(r — u)du 22) 
R(r) = U(r — a [C1 + a + 30%) e P(r — a) 
+ ae*P’ (7 — a)| 
+ U(r — a) f Kilo u)P'(r — u)du. (23) 


‘ 


t 
; 


With regard to the function II(7), note that it is essen-— 


tially a “smeared out” version of P(r). The main effects” 


of this smearing may be estimated by differentiating 


(20) which gives II’ = P’ — 211. Thus II’ =0 when2 Il=P’, 
and it follows that (in normalized units), 


Meee = FP axe 
If P(¢) does not change appreciably in times of the order 


of 27, the equality in (24) will be a good approximation. | 
In this case, of course, there are no asymptotic dc | 
fields unless P(r) increases like 7 indefinitely. As before, — 


we can derive the asymptotic relations for r>~, 
Si(r) ~ a?P’(r) 
So(r) ~ 2P"(7) 
R(r) = P’(r). 


ATTENUATION AND SIGNAL VELOCITY 


(24) 


In considering the practical implications of the above | 


results, it must be recalled that distances and times 
have been expressed in units of 2vT and 27, respec- 
tively. For sea water 2T=3X10- sec and 2vT £1 cm. 


Thus the numerical values of a in the above equations © 


are essentially equal to the range in centimeters, and at 
1 km range, for example, a=10°. These large values of 
the variables have an essential effect on the general be- 


havior of the functions K,,(a, 7) defined in (12), and it is 


necessary to examine their properties in some detail. 


We shall show that the integral terms in (13)—(19) and : 


(21)—(23) determine the peak values of the fields. 

When 7 =a, the fields first start to rise, and the values 
of the K functions are K, = (a?/2)e~* and Ky =(a?/8)e-, 
but because the modified Bessel functions grow expo- 
nentially, the K functions gradually rise far above these 
very low initial values. When 7 is very large, we may 
use the asymptotic expression,® 


I,(Z) ~ (awZ)~'2@2 (25) 
which'is Valid (@orreallay when : 
Z> (an? — 17/8. (26) 


° G. N. Watson, “Bessel Functions,” Cambridge University Press, | 


New York, N. Y., p. 203; 1952. 
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Substituting in (12), the asymptotic form of K, (for 
fixed a and large 7 —a) is expressible in the form 


are—@ /2r 


K,(a, 7) ~ emia : 


(27) 
(In fact, this is even a fair approximation when T=6G; 
it then contains a factor e-*/?2 and is, therefore, exceed- 
ingly small at ranges of interest in most practical prob- 
lems.) Note that the final decay of K, is not exponential 


_ but rather of the form 7-*, 


To locate the intermediate maximum in K,, it is 
most convenient to assume at first that it coincides with 
the maximum of (27), which occurs at 


Tm = a3 /(20:+ 1). (28) 


Fora >30 (range > 30cm approximately), Z = (Tm?—a?)1/2 


satisfies (26) more than adequately, thus justifying the 


_procedure. By (27) and (28), the maximum values of 


the K,, are then 


CB eneve = 0.462, 6 Qe = 1.83/40? (29) 


occurring at 7 =a?/3 and a?/5, respectively. 

To illustrate the peculiar behavior of these functions, 
consider a range of about 1 km in sea water, a=10°. 
Since a?/3 =10!°/3 and 2T is about 3X10-" sec for sea 
water, we see that, while K; and Kz start to rise (the 
signal “precursor” arrives) at about 30 usec, the peak 
does not occur until about ¢=1 second! The initial, 


“precursor” jump in K, is virtually zero, namely 


(a3/2)e-* = 104-999 as compared to the maximum value, 
~0.462; note also that the peak of Ky» is only about 


~ 10-!° times that of K;. Even for a range of 1 meter 


(a=10?) in sea water, (Ki)max = 2500 (Ke)max and these 
peaks occur at 1 usec even though the “precursor” 


_ (magnitude about 10~-*%) arrives at 1/30 psec. Detailed 


values for a=10! (100 m range in sea water) are given 


Bin Fig. 1 and Table I. 


It is clear from these examples that for ranges of prac- 


r tical interest, at least in sea water, K2 is usually much 


smaller than K, and that, in any case, the integral terms 
in (13)—(15) and (21)—(23) will be much larger than the 


_ other terms, which involve factors of e~*. In this sense, 
the K functions are essentially the “delta function re- 
sponse” of the conducting medium. 


Thus consider a brief pulse (either electric or mag- 


netic dipole) of width W and with a more or less tri- 


angular shape. By “more or less triangular shape” we 
mean that the maxima of the derivatives of M or P are 


approximately given by 


M,,' = M/W, (30) 


M,,/' = M,,/W? 


where the subscript m denotes maximum value. If the 
pulse width W satisfies 


W/T K @? (31) 


then the kernal functions K,, are essentially constant in 
regions of width W/T near their maxima, and the inte- 
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a= 10% 


(Y= 100m in sea water) 


: 
~ 


7 «107 (FOR SEA WATER, t/3, msec) 


Fig. 1 


TABLE 1 


VALUES OP THE FUNCTIONS Ki, Ke FOR a=10! (100 METER 
RANGE IN SEA WATER) 


t, msec for 
5 us Ks Sea Water 
10! 5 X 1074332 1.25 X 10-4382 0.003 
108 7.6X10-20 7.6X10-%6 0.3 
3X 10° 4.510% io Ome 0.9 
6X 108 0.00652 1.09107 1.8 
107 0.0850 8.51079 & 
1.5107 0.245 1.631058 4.5 
2X107 0.368 1.841078 6 
2.5107 0.432 Tais lOs* 7.5 
3<107 0.459 153X108 9 
3.33107 0.462 1.391078 10 
4X107 0.452 Tels 1085 12 
5x<107 0.415 8.03 X10~° 5 
6xX107 0.374 6.23 X10~° 18 
8X<107 0.298 See LO 24 
108 0.242 2.421079 30 
2108 0.110 ed Opa One 60 
5X<108 0.0323 6.4610! 150 
10° 0.0120 Le ZO KL Ome 300 
3109 0.00239 7.98 X10-8 900 
101° 0.000398 3.98104 3000 


gral terms in (13)—(15) and (21)—(23) can be approxi- 
mated by taking the K functions outside the integral. 
In this way, for example, it follows from (13) and (29) 
that the peak value of G(¢). is approximately 
4 (Mm! +2Mm) = Mn(it1/2W) in normalized units or 
M,,(1+T7/W) in arbitrary units. In this way the re- 
maining functions can be approximately evaluated and 
the peak fields found. The results are (in ordinary mks 
units) : 

Magnetic dipole: 


) ee = (M,,/ror')ro Dx Mo (32) 
ie. 
Bmax = (uM,,/4rr°) (1 a -) ro X (ro x M>) 
+ (uMnW/mor®)[3r0(r0-Mo) — Mo]. — (33) 


Electric dipole: 


Fax = (P»/4roWr’*) ro x (ro x P») 
+ (Py,/2moWr3)(20T/r)?|3r0(r0- Po) — Po| (34) 


Baw = (P3/ror)( Po X80). (35) 
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These approximate values are valid for a pulse width, 
W<T(r/2vT)?=por?/2 and for a range, ry >30(2vT) 
=(60/c) (e/u)"?. The second terms in (33) and (34) 
are usually much smaller than the first terms. 

The results of (32)-(35) show that even the peak 
fields from a short pulse are attenuated algebraically 
rather than exponentially, in contrast to cw signals. 
Moreover, as already anticipated, the arrival times are 
surprisingly late. According to (28), the peak values 
(32)-(35) occur approximately when 


t = tm = (T/2)(r/20T)? = por?/4. 


This arrival time is perhaps more strikingly represented 
in terms of an average “practical” signal velocity, 
(vert = 1/tm = 8v2(T/r) = 8/rop. (36) 


For r=1 km in sea water, this is 1.3 X 10? m/sec or about 
sonic velocity. 
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April 


CONCLUSION 


Exact expressions have been derived for the fields — 


excited in an infinite conducting medium by magnetic 
or electric dipoles with arbitrary time variations. These 
rather unwieldy results are then studied for the special 
case of brief pulses in media such as sea water. Simple 
estimates for the peak transient fields at moderate to 
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very large ranges are presented, which show that the | 


attenuation with distance behaves like r~* or r~™ instead 
of the very strong exponential cutoff characteristic of cw 
signals. The time of arrival of these peak fields is ex- 
pressed as a mean effective “practical” signal velocity 
which varies with range and can be as low as sonic 
velocities. 

By regarding a more general time—and space—vari- 
able source as a superposition of pulsed dipoles, these 
simple results can be used to obtain rapid approximate 
estimates of the important practical features of the re- 
sulting signals. 


Determination of a Current Distribution over a Cone Surface 
Which Will Produce a Prescribed Radiation Pattern* 
H. UNZt 


Summary—The problem of determination of the current dis- 
tribution over a cone surface, which will produce a prescribed radia- 
tion pattern, is formulated. F pattern and G pattern are defined and 
expressed in terms of a “potential” function #, in case of currents 
polarized in the direction of the generating lines of the cone. 

A matrix relationship is found between coefficients of a series ex- 
pansion of and coefficients of a series expansion of the currents. 
The relationship becomes a direct one to one correspondence in case 
of an infinite biconical system. 


INTRODUCTION 


HE problem of finding a current distribution along 
“ae surface of a cone in order to produce a pre- 
scribed far zone radiation pattern has become very 
important lately. While the mathematical formulation 
of the problem is not difficult and is based on already 
existing ideas, the method of solution in this paper is 
rather new. 
The method of excitation of the required current dis- 
tributions over the surface of the cone is not discussed 
herein. It is assumed that once the required current 


* Manuscript received by the PGAP, September 30, 1957; re- 
vised manuscript received January 27, 1958. This work was supported 
by the U. S. Navy at the University of California, Berkeley, Calif., 
under Contract No. N7onr-29529. 

+ Elec. Eng. Dept., University of Kansas, Lawrence, Kan. 


distribution has been found, it can be excited over the 


surface. The boundary value problem, which is quite 


difficult, also is not considered in this paper. We are in- 
terested here in the relationship between the far zone 
radiation field and the current distribution. 

It is known that a problem of this kind is not unique, 
namely there are infinite numbers of current distribu- 
tions which will produce, within a certain approxima- 
tion, the same far zone pattern. By expanding the cur- 
rent and the far zone field in two double series, matrix 
relationships between the coefficients of the expansions 
can be found. Any two sets of coefficients, which obey 
those relationships, will be a solution of the problem 
of finding the far zone field which corresponds to the 
current dustribution or vice versa. In the case of an in- 
finite biconical system, we have a one to one corre- 
spondence between the coefficients of the expansion of 
the far zone pattern and the coefficients of the expansion 
of the current distributions. It is similar to the result 
derived in the much simpler case of continuous line cur- 
rent distribution and its far zone pattern.! 


+H. Ung, “Linear Arrays with Arbitrarily Distributed Elements, ” 
University of California, Berkeley, E. R. L. Rep. Ser., No. 60, 


Issue No. 168, Navy Contract No. N7onr-29529, pp. 54-56; Novem- 
ber 2, 1956. 
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P(R;6;) 


K=1(2,¥) i 


FORMULATION OF THE PROBLEM 


Spherical coordinates, with the origin at the apex of 
the cone (Fig. 1), are used. The surface of the cone is 
given by 0=6. The coordinates of the observation point 
P are P(R; 6; @). The coordinates of the source point Q 
are O(p; 40; y). 

The far zone field is found by using the general form- 
ulas given by Silver? 


E i can [ |G i) 
gt ae e° 46 
? 4rR Vv 
te /~ Onis) eika-ingV (1a) 
b 
10 ca 
E. = —_— cn f | Gem) 
' 4irR v e 
_ /~ TInt) eita-indV (1b) 
MK 
where 


ad=radius vector from the origin to source point Q, 
J,=electric current density vector, 
Tm =magnetic current density vector, 
ir; ie; i»=spherical coordinates unit vectors at 
RRO Gy 
k=21r/N=w/c t=VJV-1, 


with harmonic time variation e+??'. 
From (1) it is seen that there is a complete similarity 
between the far zone fields received from electric cur- 


2S, Silver, “Microwave Antenna, Theory and Design,” McGraw- 
Hill Bake Con Inc., New York, N. Y., M.I.T., Rad. Lab. Ser. No. 
12, p. 89, (134); 1949. 
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rents (J.) and the far zone fields received from magnetic 
currents (J,,). In order to discuss both of them in gen- 
eral, and since we discuss currents distributed over the 
surface of a cone; 7.e., surface currents K, rather than 
volume currents J, let us define the radiation pattern 
functions as 


F(@;) = if JK nerras (2a) 


G(9;¢) = i (K -ig) ete ind S, (2b) 
S 
F(@; ) is called the F pattern and G(@; ¢) is called the 
G pattern. In case of electric surface currents, the F pat- 
tern will be the radiation pattern of Es and the G pattern 
will be the radiation of Ey. In case of magnetic surface 
currents (or slots), it will be the other way around. 
The radius vector @ from the origin to the source 
point on the surface of the cone Q(p; 00; W) is a straight 
line along a generating line of the cone and may be 
written 


ad = pip. (3a) 
The surface current density K is polarized in the di- 
rection of the generating lines 


K = I(p; W)t, 


where J(p; W) is the current distribution over the surface 
of the cone 6 =6. 

Substituting (3) in (2) and integrating over the cone 
surface we get 


(3b) 


I 


F(@;¢) = sin af 1(p;¥)av |” poli, ine (4a) 


2a pl Latest 
G(0, @) = sin [10s Wav f  pdolin-tedero-e, (4b) 
9 Po 
RELATION BETWEEN F PATTERN AND G PATTERN 


From the vector analysis identities it may be easily 
shown’ that 


(6) 
te = 19 (5a) 
00 
0 . 
— tr = i, sin 0. (5b) 
dp 
Since p and 7, do not depend on @ and ¢, we can write 
or 
oe etkp (ptr) = etke (tp *R) i Rp (=) (6a) 
06 06 
04 
| etke (ip 22), — erke (ip . iR)ikp (: A =) ; (6b) 
dg dg 


3H. B. Phillips, “Vector Analysis,” John Wiley and Sons, Inc., 
New York, N. Y., p. 98, problem 14; 1933. 
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Substituting (5) in (6) and comparing with (4), we 
can write 
0, sin’ by G2" CUS mel ee 
FO; 6) =~ 2 f “rep; wan | dpetesem (a) 
00 tk 0 Po 


as be 0 sin 4 
(6; $) sin as 


Qa Py 
J I(p; wav J dpete(e -t).(7b) 
0 Po 


From (7) we see that we can write the far zone field in 
terms of a “potential” function ®(6;¢) 


0 
F(; 8) =~ (0; 9) (8a) 
G(6; ¢) sin @ = es (0; d) (8b) 
dg 
where 
SimuOon fee YEA 
&(4; ¢) a ik if I(p;~)dy ip dpe**e (% tk), (9) 


From (8) we get the identity 


* 70,4) =~ (G0; ) sin 6] (10) 
== -¢) = — ; p) sin 6]. 

Od : 00 

Eq. (10) gives a direct relationship between the F pat- 
tern and the G pattern. 


DETERMINATION OF THE CURRENTS 


From (8) we see that by giving the “potential” func- 
tion ®(8; @) we prescribe both the F pattern and the 
G pattern. In order to determine the current which will 
produce a prescribed radiation pattern, the integral 
(9) must be solved [®(6; ¢) is a given function and 
I(p; W) is the unknown current distribution function 
over the surface of the cone]. We solve this problem by 
transforming the integral (9) into infinite matrix rela- 
tionships between the coefficients of the expansions of 
@(6; @) and J(p; y). Since we are interested only in a 
radiation pattern which will approximate the required 
one, (9) could be solved by cutting off the infinite 
matrixes arbitrarily and inverting them. 

It may be found that 


ir = sin 6(7, cos ¢ + i, sin d) + 7, cos 0 (11a) 
i, = sin 4(%, cos + 7, sin) + i. cos.  (11b) 
From (11) we get 
7,-tk = sin % sin 8 cos (W — ¢) + cos cos 6. (12) 
Substituting (12) in (9) we get 
smvGy, .(* 2" Py 
#0;¢) == [ros wav f "ae 
tk 0 Py 
-exp ikp|sin 4 sin 8 cos (YW — ) + cos 6p cos ot. (13) 


Expanding the current function and potential func- 
tion in a complex Fourier series as follows 
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+o 
I(e;¥) = Do talele™ (14a) 
(0; ¢) = Bmn(0)e™?. (14b) 


m=—co 


Substituting (14) in (13) and using the identity* 
Pig 
f giz cos ngimndn = Qri!™T\m(Z) | m| = integer (15) 
0 


we get, after comparing coefficients of the F ourier series 
on both sides,* 

sin 0 py 
Bm(@) = 2m sin A aim f bn(p)dpe*? cos 69 cos 8 

ik x 


J\mj(Rp sin 6) sin @). (16) 


Expanding B,,(@) in terms of associated Legendre 
polynomials 


Bn(6) = >, Da™P,!™\(cos 6). (17) 
n=0 
Using the orthogonality properties® 
f P,,'™'(cos 6) P,'™!(cos 6) sin 6d6 
0 
2 n+ |m|)! 
Te | | nt» (18) 
2n+1 (n— | m | )! 
the coefficients in (17) become 
pags ! 
nen es (n | m|)! 
2 (n+ |m|)! 
. il Bn(0)Pn'™!(cos 6) sin 6d6. (19) 
0 
Substituting (16) in (19) and using the identity® 
if etkR cos a cos 07, (RR sin a sin 6)P,,(cos 6) sin 6d0 
0 
= 2777, (RR)Pnr™(cos a) (20) 


where j,(kR) =spherical Bessel function; we get 

2rsin®  In+1 (n — | m|)! 
1 m 

ik 2 (w+ |m|)! 


D,™ = 2i"P,,'™! (cos Oo) 


if "bm phn bo)dp. (21) 


Let us assume that b,,(p) may be expanded as 


_ 4E, Jahnke and F. Emde, “Tables of Functions,” Dover Publica- 

tions, New York, N. Y., pp. 149, 131; 1945. The introduction of 
|m| instead of m is a simple extension of the original identity. 

* Similar relations have been developed by: K. M. Siegel, e¢ al., 
“Studies in Radar Cross Sections X XII—Elementary Slot Radiators,” 
University of Michigan, Ann Arbor, Mich., E. R. I. Rep. 2472-13-T, 
p. 16; November, 1956; and G. Held, University of Washington, 
Seattle, Wash., private communication. F 

J. A. Stratton, “Electromagnetic Theory,” McGraw-Hill Book 
Co., Inc., New York, N. Y., p. 403, (16) and (17), p. 411, (69); 1951. 


. 
> 


bm(o) = Y~ di™7 (kp) (22) 
l=0 
and the series (22) converge for both large p and J. Let 


us denote the definite integral of spherical Bessel func- 
tions: 


. 


f “n@ineas = Tae: #003) 


Tv 
' ~ [Q+ Qn + Hp? 
where I7,, (20; 21) is a definite number. 
Substituting (22) in (21) and using (23) we get 


2? (n—|m| )! 


D,™ = intimi—1[ 2+ 1] 1/2 
(n+ | m|)! 


Be sin OP,,|"!(cos 00). 


=. Ti,m(p0} p1) 
2 ——_§_—$<— q,™), 
2 [27+1]1/2 : (24) 


Since P,,'"!(cos 8) =0 when n< | m|, we see from (24) 
that D,™ =0 when x<|m|. Taking this into account 
we have to rewrite (17) as follows 


Br(6) = >> Dr™P,!™I(cos 6). 


n=|m| 


(25) 


MATRIX RELATIONS ® 


The F pattern or the G pattern need not be exact. 
Therefore the “potential” function ®(6; ¢) may be ap- 
proximated, due to (14b) and (25), as follows 


M N 
6(0;¢) = >> > D,™P,|™\(cos Oem 


m=—M n=|m| 


(26) 


and the current distribution may be approximated, due 
to (14a) and (22), as 


Tio3¥) = De De difilkp)e™. (27) 


m=—M I=0 
The relationship between the coefficients D,“ and 
d,™ is given by (24). 
In Appendix I we prove that 

Tim ©; + 0) a 57,m 
In case of infinite biconical system, (13) will have to be 
written twice, for 6) and t—. If we assume that the 
current distribution in one cone is symmetric with re- 
spect to the apex of the biconical system to the other 
cone or 


(28) 


I2(—p;¥+ rn) =I%(0; Y), (29) 


we will have for the infinite biconical system the same 
type of relation as (13) with po= — © and pp=+o. 
Substituting (28) in (24) we get 


ee, | m|)! 


re | | 1 sin 0.P,.'"!(cos 00)dn™. (30) 
k? n m|\)! 


"2% (m) — 


Eq. (30) gives a one to one correspondence for the in- » 


finite biconical system between the pattern coefficients 
given by (26) and the current coefficients given by (27). 
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The result in (30) is rather interesting. In case of 
n<|m| we see that since P,'"\(cos 6)=0, we get 
D,™=0 and d,™ for n<|m| may be completely 
arbitrary. d;™ in (27) is determined for 1=|m| but is 
completely arbitrary for 1<|m| . This is a proof that 
the solution is not unique and there are different current 
distributions along the surface of the cone which will 
produce the same radiation pattern. 

Eq. (24) is true for every n and every m, so it may be 
rewritten in a matrix form 


Do ™ | = [Min] dim (31) 


where (25) represents infinite matrices for every m. The 
terms of the matrix [Mz,,°™] for every m become small 
very fast from a certain value / or 7 on, as may be seen 
from (24) and (23) and from the behavior of the spheri- 
cal Bessel function j,(z). Therefore, the infinite matrix 
[Mi.™] may be cut off to become a finite matrix for 
approximation purposes. ; 

If we take in (31) | m| =constant, we get a matrix re- 
lationship for each | m| . According to (24), D,™ =0 for 
n< | m| regardless of d;™. In other words, if we use the 
relations in (26) and (27), we want to find (P+1) un- 
known d,™ coefficients for each | m| , but we have only 
(N— | m| +1) independent equations with given D,™. 
If we want to take the number of terms to be equal 
N=P, obviously we have arbitrariness, since we have 
| m| more unknowns than equations. We might take 
N=P+ | m| and then we will have the same number of 
equations as the number of unknowns; the matrix in 
(31) will then become a square matrix and the formal 
solution of P equations with P unknowns may be writ- 
ten as 


ad) = [,,)-D, >]. (32) 

The prescribed far-zone “potential” function ®(6;@) 
should be such that it is possible to approximate it uni- 
formly by the series given in (26), since the series is the 
most general expression in spherical coordinates. Other- 
wise the required current distribution over the surface 
of the cone cannot be found. In case of a finite cone sys- 
tem, we have assumed that there are no currents along 
the base of the cone. 

The results found for the continuous current distribu- 
tion may be transformed to the case of discrete array 
elements over the surface of a cone. We will have to ap- 
proximate the results for continuous currents by sets of 
discrete elements. 


APPENDIX | 


Kaptéyn’ derived the following relation in Bessel 
functions 


8 dz 2 sin (4 = v)r/2 
if PUA AEC agate 2 2 
0 Zz Tv 


[i 


(33) 


7G. N. Watson, “A Treatise on the Theory of Bessel Functions,” 
Cambridge University Press, Cambridge, Eng., p. 404(ff); 1952. 
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where pu; v are real numbers. Substituting 
yv=nt+4 


where m:n are positive integers or zero, in (33) and re- 
arranging we get 


=m +> (34) 


© : sin (m — n)w/2 
— Jmii(2)Inpi1/2(2)ds = (35) 
if etal (z) +1/2(2) (m+ 2)? — 4p 
We have by definition 
W 1/2 
jn(2) > [=| I nxaj2(2) (36) 
22 
where 7,,(z) is the spherical Bessel function. 
From (36) we can rewrite (35) as 
4 sin (m — n)r/2 
im(Z)In(Z)dz = . (37) 
if jul) u(2) (m — n)(m+n-+ 1) 
From the power series for j,(g) we can see that 
G2) = (—) 7,2). (38) 
Let us consider the integral 
Tv 
In.n(— eo) 5 se 2 ) 
(2m + 1)(2n + 1)]!/2 3 
=f inlsiinlsids 9) 


—o 
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If both m, n are even or both m, n are odd integers, 
Inn(— ©, +) will vanish according to (37) because 
(m—n) are even integers. If m=odd and n=even Or 
vice versa, Imn(— ©, +) will vanish according to 
(38), because the integrand in (39) is an odd function. 
Finding the limit ’ 


Fost sin (m — n)1/2 ae (40) 
mn (m—n)\(m+n+1) wnm+i1 
we see that Ba 
is DW (coe + 2) 
[(2m + 1)(2n + 1)]!/? 
ie - (41) 
== ie lis i= Orns 
[ ntojntoae = 5 
where m; n>0. Eq. (41) may be rewritten as 
Lipghs 5 ste 00 ) i On (42) 
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The Effects of the Physical Parameters on the Bandwidth 
of a Folded Dipole* 


J. P. GERMAN} anv F. E. BROOKS, JR-t 


Summary—An expression is derived which yields the feedpoint 
admittance of a folded unipole or dipole at any desired frequency. 
This expression shows that the input admittance is composed of two 
parts—a radiation part and a transmission line part. The radiation 
part is analyzed as a vertical radiator while the transmission line part 
is analyzed by transmission line theory. The terms of this expression 
allow the input admittance of the antenna to be studied in terms of 
the physical parameters. By plotting the feedpoint admittance on a 
Smith chart it is possible to observe at what frequency the currents 
in the feeding transmission line exceed some desired standing wave 
ratio and thus obtain a measure of the bandwidth of the antenna. 

The findings show that the broad-band characteristics of the 
folded diple are obtained primarily through the increase in equivalent 
cross section caused by the folding of the two elements. The band- 
width is further aided by the reactive component of the transmission 
line part being opposite to the reactive component of the radiation 
part, even though the reactive component of the transmission line 


* Manuscript received by the PGAP, May 14, 1956; revised 
manuscript received, May 13, 1957. 

+ School of Elec. Eng., Purdue University, Lafayette, Ind. 

t Collins Radio Co., Dallas, Tex. 


part is considerably smaller than desired. Experimental and cal- 
culated results show that for a given size conductor the bandwidth is 
improved as the spacing between conductors is increased. The spac- 
ing should be as wide as possible but still provide the desired field 
pattern and input admittance. The bandwidth of a folded dipole is 
improved if the conductor size is increased regardless of the spacing 
between elements. Maximum bandwidth is achieved by using large 
conductors with wide spacing. 


INTRODUCTION 


ARIATIONS in the conductor size and the spac- 

\ ing of a folded dipole will obviously influence the 

feedpoint impedance or admittance and, as a re- 

sult, will influence the bandwidth of the antenna. It is 

the purpose of this paper to determine the best arrange- 

ment of the physical parameters in order to achieve a 
maximum bandwidth. 

A measure of the bandwidth of any antenna may be 

had if the feedpoint impedance or admittance is plotted 
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GROUND PLANE 
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Fig. 1—A folded unipole. 


ELEMENT TWO ELEMENT ONE 


HO) 
OE 


Fig. 2—A folded unipole with coaxial feed replaced by equiva- 
lent generators. 


on a Smith chart! so that it may be observed immedi- 


_ately at what frequency the currents in the feeding 


transmission line exceed some desired standing wave 


ratio. In order that the Smith chart plot may be made, 
it is necessary to have an expression which will provide 


the feedpoint impedance or admittance not only at 
resonance but at any frequency off resonance. The feed- 
point impedance at resonance has been undertaken by 
Roberts? and Guertler,’? but their findings do not pro- 
vide sufficient information for plotting the feedpoint 
impedance at frequencies other than resonance. The 
first undertaking of this paper will be to derive an ex- 
pression for the feedpoint admittance of the folded di- 


pole as a function of frequency. Once this expression has 
~ been derived it will then be possible to examine what 
-_ effect the conductor sizes and spacings will have on the 
_feedpoint admittance and the bandwidth of the an- 


tenna. 


THE FEEDPOINT ADMITTANCE OF A 
FOLDED UNIPOLE 


The treatment of the folded dipole will be analyzed 
from the standpoint of a folded unipole, 7.e., the antenna 
resulting when an ordinary folded dipole is cut in half 
and fed against ground as shown in Fig. 1. The opera- 
tion and admittance characteristics are exactly the same 
as for the common folded dipole except that the admit- 
tance is doubled. Any admittance multiplying factor 
due to the folding will be the same for the folded uni- 
pole as it would be for the folded dipole. 

In order to provide a suitable method of analysis it 
will-be convenient to replace the coaxial feed on this 
folded unipole by three radio-frequency generators as 
suggested by Roberts.” Since these generators are purely 
imaginary, it can be stipulated that they have zero in- 
ternal impedance. All generators are operated in phase, 


1p. H. Smith, “An improved transmission line calculator,” 
Electronics, vol. 17, pp. 130-133; January, 1944. ; i 

2W.V.B. Roberts, “Input impedance of a folded dipole,” RCA 
Rev., vol. 8, pp. 289-300; June, 1947. navies ; 

3R. Guertler, “Impedance transformation: in folded dipoles, 
Proc. IRE, vol. 38, pp. 1042-1047; September, 1950. 
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ELEMENT TWO ELEMENT ONE 


Fig. 3—The lower generator acting alone. 


ELEMENT TWO ELEMENT ONE 


ir] Ne 


Sonal 


Fig. 4—The upper two generators acting alone. 


with polarities indicated by the arrows in Fig. 2, two 
with voltages e and one with voltage e. From Fig. 2 it 
can be seen that the generator on the left and the lower 
generator have opposite polarities, placing the lower 
end of element two, the grounded element, at ground 
potential as in Fig. 1. The generator on the right and the 
lower generator, operating in phase and with polarities 
adding, apply a voltage (e+e,) =er to the lower end of 
element one, the driven element. Thus, Fig. 2 is equiva- 
lent to Fig. 1. 

The extra generators have been placed in the prob- 
lem so that if the individual generator currents can be 
determined, they may be added, by means of the princi- 
ple of superposition, to obtain the current entering the 
lower end of the driven element. 

Assume for the moment there is only voltage on the 
lower generator (see Fig. 3). This radio-frequency gen- 
erator is then feeding a radiator composed of two ele- 
ments in parallel. If Ye is the input radiation admit- 
tance of this composite radiator, then the lower genera- 
tor will supply a total current of eYz to this antenna. 
The current that is produced by the lower generator 
acting alone will divide inversely between element one 
and element two as the impedances it sees in the ele- 
ments. The current flowing in element one, the fed ele- 
ment, will be designated as Jiz and the current flowing 
in element two, the grounded element, will be desig- 
nated as Jor. Note that both currents are instantane- 
ously flowing in the same direction. From this it is possi- 
ble to write 


(1) 


Now assume there is no voltage in the lower generator 
(see Fig. 4). Under this condition, the voltage (e+e) is 
impressed upon a transmission line shorted at the far 
end. If the two conductors are unequal in size they will 
present unequal impedances to the two generators; 
therefore, the generators must be unbalanced in order to 
place the point between the generators at ground po- 
tential, and to make the currents in the two conductors 
equal and opposite. The current in this transmission 
line will be designated as Jr and this current has an 


Tin + Lor = GV gt 
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Fig. 5—Ratio of currents in the two conductors of a folded dipole. 


instantaneous direction outward on one conductor and 
inward on the other. From transmission-line theory it is 
possible to write: 


Ir = (e: + e)Y, coth yl, (2) 


where Y, is the characteristic admittance of the trans- 
mission line formed by elements one and two, / is the 
distance from the shorted end and y is the propagation 
constant, generally denoted as y=a+ j8 where a is the 
attenuation constant and 8 the phase-shift constant of 
the line. 

The coaxial feedpoint admittance may be found by 
taking the ratio of the total current entering the lower 
end of the driven element to the total voltage applied 
to the lower end of this same element. 

The feedpoint admittance may then be written as 


Tir + Ir 
me inertredy 
and defining the current ratio as 
Tor 


woe! 


Thre 


Y ) (3) 


(1) becomes 
Vre 


= prea 


It will be assumed that since both elements are in 
parallel, in close proximity to each other, and generally 
under similar conditions, ir and Jr are in phase, and 
hence 7 is a real number. The ratio of the currents in 
two parallel conductors of unequal size can be solved 
using logarithmic potential theory. Fig. 5 shows the 
results of this analysis with dimensions of the conductors 
as shown. 

Now substituting (2) and (4) in (3) the feedpoint ad- 
mittance becomes: 


1 € 
ip 4 - Y.coth yl. 
Hera) eae oh 


Tir (4) 


(3) 
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Fig. 6—Radius of equivalent cylindrical antenna. 


It will be noticed that the first portion of this admit- 
tance expression, the radiation portion, has three parts. 
The term Yp is simply the admittance that a radio- 
frequency generator would see when connected to a 


‘vertical radiator which is fed against ground. The equiv- 


alent cross section of the two elements in parallel may 
be solved for by assuming the conductors are parts of 
one-wire transmission lines far from the earth, and hav- 


; 


ing equal characteristic impedances. Fig. 6 shows the. 


results of this analysis. The input impedance of a simple 
vertical radiator fed against ground will not be reviewed 


in this paper, but the results of previous findings‘ will _ 


be used later. 

Following the Yr part of the feedpoint admittance 
expression there are two factors which, when multiplied 
into the term, modify this admittance. This modifica- 
tion is, of course, introduced by the two elements of the 
folded unipole, the conductor size and the spacing be- 
tween them. The voltage ratio factor [e/(e+e1)], is 
dependent upon finding the neutral plane between the 
elements of the antenna, as shown in Fig. 7. The solu- 
tion for the voltage ratio [(e/e+e:)], may be shown to 
be proportional to the characteristic impedances of the 
individual conductors against the neutral plane. The 
results of this analysis are shown in Fig. 7. 

As a final portion of the feedpoint admittance there is 
an added term due to the transmission line component. 
If the transmission line is considered lossless the term 


_ *E. Hallen, “Admittance Diagrams for Antennas and the Rela- 
tion Between Antenna Theories,” Cruft Laboratory, Harvard Uni- 
versity, Cambridge, Mass., Tech. Rep. No. 46; June 1, 1948 (Office 
of Naval Res., NSORI-76 Contract, Task Order No. 1, NR-078-011). 
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Fig. 7—Cross section of the transmission line. 


Y,. coth yl may be replaced by —jY, cot 27/1/X where d 
is the wavelength on the line. The reactive component 
of this term has a sign opposite to the reactive com- 
ponent of the first part and this is one factor which aids 
in making the folded dipole more broadband than a 
simple dipole. 


Tue FINAL FEEDPOINT ADMITTANCE EXPRESSION 


It is now possible to write the feedpoint admittance 
of a folded unipole at any desired frequency. 


v : ( 1 )( Loe ) ee eae 2nl (6) 
= ee | = Ve Col) 
og a 1 + n Lo2 ar Zu A d 
_ where: 
; 18! 
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Yru =feedpoint admittance of the folded unipole. 

Yr=feedpoint admittance of a composite radiator 
made up of the two elements in parallel. 

n =ratio of the currents in the two elements. 
D=distance between centerlines of the two ele- 
ments. 

R,=radius of element one, the driven element. 

R2=radius of element two, the grounded element. 

Zu =characteristic impedance of element one and a 
zero potential plane between elements one and 
two. 

Zo2=characteristic impedance of the element two 
and a zero potential plane between elements 
one and two. 

Yc=characteristic admittance of the transmission 
line formed by the two elements. 

l=distance from the shorted end of the trans- 
mission line. 
\=wavelength. 


To determine Yr from published literature, it is 
necessary to calculate the equivalent radius of the two 
elements in parallel. This equivalent radius is given by 


R= (Ri D*)¥ (a+) 


where: R’=equivalent radius of the two elements in 
parallel. Other symbols are as defined above. 


ACCURACY OF THE ADMITTANCE EXPRESSION 


The feedpoint admittance expression now allows the 
normalized input admittance of a folded unipole to be 
plotted on a Smith Chart so that a measure of the band- 
width can be obtained. This allows an immediate ob- 
servation to be made of the frequency at which the 
feeding transmission line currents exceed some desired 
standing-wave ratio. 

The results of measured admittances vs calculated 
admittances will be presented so that some idea of the 
accuracy of the feedpoint admittance expression may 
be obtained. Many examples have been calculated and 
measured, however only two curves will be shown since 
comparable results were obtained from all examples 
checked. Fig. 8 shows the calculated and measured 
(normalized to 6.66 millimho characteristic admittance) 
feedpoint admittance of a folded unipole made of two 
equal size conductors, and also the dimensions. 

In order to apply the final feedpoint admittance ex- 
pression it was necessary to calculate the equivalent 
cross section of the composite antenna composed of the 
two conductors in parallel. With this new cross section 
it was possible to go to available literature* and obtain 
Yr, the input admittance of a simple cylindrical ver- 
tical radiator having the calculated cross section. The 
other parts of (6) were then calculated which in turn 
modified Yp to fit the physical dimensions of the exam- 
ple under consideration. The results of these calcula- 
tions were then checked by actual input admittance 
measurements. The measured and calculated data agree 
reasonably well. 
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Fig. 8—Calculated (X) and measured (0) feed point admittance of 
a folded unipole, normalized to 6.66 millimhos. 


Fig. 9 shows the calculated and measured normalized 
feedpoint admittance of a folded unipole made from 
two unequal size conductors. As may be seen from the 
dimensions on Fig. 9 one conductor has a diameter 
twice that of the other. Again the calculated and meas- 
ured data agree very well. 


THE BANDWIDTH CONSIDERATION 


A sample calculation of the input admittance of a 
folded unipole will show that the susceptance of the 
radiation portion [the first term of (6)] is generally 
much larger than the susceptance of the transmission 
line component [the second term of (6) ]. It would ap- 
pear that maximum effectiveness would be obtained by 
making the characteristic admittance of the transmis- 
sion line, Y., as high as possible. 

If the characteristic admittance is made larger by 
decreasing the spacing D of the transmission line, the 
equivalent cross section of the composite radiator is re- 
duced thus resulting in a reduced bandwidth. This un- 
fortunate situation leaves only one method of increasing 
the characteristic admittance—that of increasing the 
size of the conductors. Calculation of a few examples 
soon shows that once the conductors have been in- 
creased in size, the bandwidth can still be increased by 
increasing the space between the conductors even 
though this might appear to be an ineffective technique. 

It should be pointed out that neither conductor size 
vs bandwidth nor spacing vs bandwidth follow a linear 
or simple mathematical relation. The conclusion which 
must be drawn is that the increased bandwidth of a 
folded unipole comes about primarily through the 
increased equivalent cross section caused by the two 
elements rather than the cancelling of the reactive com- 
ponent by means of the transmission line term. The 
cancelling of the reactive component certainly con- 
tributes to the increase in bandwidth, but this term is 
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Fig. 9—Calculated (X) and measured (0) feed point admittance of 
a folded unipole with unequal size conductors normalized to 6.66 
millimhos. 


not sufficiently large to become the controlling factor. 
The equivalent cross section is increased most easily 
by increasing the spacing between elements. The limit 
on the spacing between elements will in general be de- 
termined by other factors. The field pattern is altered 
as the spacing is increased and the input admittance 
becomes larger as the spacing increases and may reach 
a value too large for the transmission line feeding the 
antenna. 


CONCLUSION 


The analysis of the broad-band characteristics of the 


of 
April 


4 


folded unipole and thus the corresponding folded dipole 


have led to the following conclusions: 


1) The broad-band characteristics of the folded uni- 
pole come about primarily through the increase in 
equivalent cross section caused by the folding of 
the elements. 

The broad-band characteristics are aided by the 
reactive component of the transmission line part 
being opposite to the reactive component of the 
radiation part; however, the reactive component 
of the transmission line part is considerably 
smaller than required for complete compensation. 
3) Experimental and calculated results show that for 
a given size conductor the bandwidth is improved 
as the spacing between conductors is increased. 
The spacing should be as wide as possible but still 
provide the desired field pattern and input admit- 
tance. 

The bandwidth of a folded unipole is improved if 
the conductor size is increased regardless of the 
spacing between elements. Maximum bandwidth 


is achieved by using large conductors with wide 
spacing. 
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The Radiation Characteristics of a Zig-Zag Antenna* 
DIPAK L. SENGUPTAT 


Summary—The radiation properties of a single zig-zag antenna 
are reported. This is a special type of broad-band traveling-wave 
antenna which, when properly designed, produces a strong axial 
beam of radiation. The radiation pattern has a half-power beamwidth 
of 22° in the E plane and 28° in the H plane; the corresponding side- 
lobe ratios are 10.5 db and 17.5 db down, respectively, in the two 
planes. Approximate expressions for the radiation fields are also given 
and compared with the experimental values. The broad-band prop- 
erty of the antenna is studied by measuring the radiation patterns 
over a range of frequencies. The performance of the antenna is com- 
pared with the Yagi antenna. The high directivity of the antenna can 
be utilized advantageously in the vhf and uhf ranges. The results of 
impedance measurement also are reported. 


INTRODUCTION 


HE purpose of this paper is to report some theo- 
retical and experimental investigations about the 


radiation properties of a single zig-zag antenna, 
first reported by Cumming.! This is a special type of 
traveling-wave antenna which, when properly designed, 
produces a strong axial beam of radiation with a very 
low sidelobe. 

In his paper Cumming studied the radiation proper- 
ties of such an antenna at a single frequency; he made 
more elaborate experimental investigations about the 
properties of a balanced type of double zig-zag antenna. 
A more thorough investigation of the single zig-zag 
antenna has been necessary in connection with the 
building up of a radio telescope whose antenna assembly 
consists of a rectangular array of six zig-zag antennas. 
For the satisfactory design of the telescope, a better 
understanding of the radiation properties, especially 
the half-power beamwidth, sidelobe ratio, and the band- 
width of operation, has been necessary. The following 
investigations were made with this end in view. 

This antenna may be considered as a special case of 
the well-known helical antenna. Theoretical expressions 
have been given for the radiation field by assuming a 
current distribution similar to that used in the helical 
case. The design condition for the axial mode of opera- 
tion has been derived directly from this expression. Al- 
though, from the nature of the assumptions made, the 
expressions seem to be very approximate, the agreement 
between the calculated pattern and the experimentally 
determined patterns, especially in the H plane, is quite 
reasonable and can be worked with for practical pur- 


poses. 


* Manuscript received by the PGAP, October 29, 1957; revised 
manuscript received, January 16, 1958. The work reported here was 
made possible by the financial assistance offered by the National 
Research Council of Canada under Radioastronomy Project, Grant 
No. G 587. Paper presented at the 1957 Canadian IRE Convention, 

Can. 
Bee co. of Elec. Eng., University of Toronto, Toronto, Can. 2 

1 W. A. Cumming, “A nonresonant endfire array for vhf and uhf, 
IRE TRANS. ON ANTENNAS AND PropaGamion, vol. AP-3, pp. 52-58; 


April, 1955. 


In the last part of this paper a comparative study is 
made between the properties of the zig-zag and the Yagi 
antennas. The better directivity and reasonably broad- 
band property of the former make this antenna a better 
choice than the Yagi in directional reception. Finally, 
the results of the impedance measurements are given. 


APPROXIMATE EXPRESSIONS FOR THE 
RADIATION PATTERN 


The single zig-zag is shown diagrammatically in Fig. 
1. The coordinate system is chosen in such a way that 
the antenna lies entirely in the Y—Z plane, the axis of 
the antenna being along the Z axis. The antenna is fed 
by a coaxial line whose inner conductor is connected to 
the antenna and the outer conductor is directly con- 
nected to the ground plane. Although the ground plane 
effects the performance of the antenna to some extent, 
especially regarding the back radiation, its main func- 
tion is to offer a proper ground for the feed and, as is 
shown later, the size of the ground plane is not critical 
for the required performance of the antenna. The theo- 
retical analysis becomes very complicated if the effect 
of the ground plane is taken into account. The expres- 
sions for the radiation field, given later, have been de- 
rived with the assumption that there is no ground plane. 

To get the field expressions, a single current wave of 
constant amplitude is assumed to travel along the an- 
tenna as is done in the case of the helix.? Theoretically 
such a current wave fails to satisfy the boundary condi- 
tions at the open end and also fails to take care of any 
possible reflection at the corners. However, if the an- 
tenna is long, compared to wavelength, the amplitude 
of the reflected wave will be very small and hence may 
be neglected for practical purposes. The current dis- 
tribution is assumed to be of the following form 


To Jet? (1) 


where, 8 =27/) and s is the distance measured along the 
antenna. The time dependence is of the form e~#°. 

Physically we can think of the antenna as a number 
of V’s connected in series. The length of each arm of the 
V and its angie should be chosen such that the radia- 
tion fron all the V’s add in phase along 6 =0° direction. 
The relation to be satisfied by the armlength (2L) and 
the pitch angle (2a) (Fig. 1) is given later. 

Once the field due to single V is known, the field due 
to the zig-zag can be found by the principle of multipli- 
cation of pattern.! The field also can be found by di- 
rect integration along the length of the antenna using 
the fundamental principles of antenna theory. The re- 
sults obtained by both methods are same. 


2 J. D. Kraus, “Antennas,” McGraw-Hill Book Co., Inc., New 
York, N. Y., p. 184; 1950. 
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Fig. 1—The schematic diagram showing the zig-zag antenna 
with the feeding arrangement. 


With these ideas in mind, if we carry out the analysis, 
it can be shown that the electric field in the Y—Z plane 
(E-plane) due to the zig-zag is given by 


(6) ; 12071 e%8" 


| —128(L—b cos jek BL{1 + sin (6 — a)} 
‘ {1 + sin (9 — a)} 

sin BL{1+ sin (a + 6)} 
{1+ sin (a + 8) 


cos (6 na a) 


cos (a + 0) | a (2) 


where L, 6, and aware as defined in Fig. 1 and b=Z sin a. 
The factor F is given by, 
1 sin mp/2 


aa ‘nm sin y/2 (3) 


where, 2 is the number of V’s, 


and 
J (= 8rb a) ; 
= | —— — — cos @]. 
" . (4) 


From (2) it is found that in order to have maximum 
radiation along @=0° direction, the following equation 
L—b= a ; 
5 
ie. L(A — sina) = d/4) S) 


must be satisfied. This is the fundamental design rela- 
tion for the zig-zag antenna. 

Proceeding in exactly the same way it can be shown 
that the electric field in the X—Z plane (H plane) is 
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given by 
E 8) 2 1207/9 elBr 
Oe BX r 
sin BL(1 + sin a cos 8) | 
— -cosa|F 
(1 + sin @ cos @) 


sin (86 cos @) 
—_—_—_——: tana 


cos 8 


Eee cos@), 


(6) 


the factor F is the same as in (3). 


EXPERIMENTAL RESULTS AND DISCUSSIONS 


A number of zig-zags (for pitch angles 2a@ equal to 
15°, 20°, 25°, 30°) each with six V’s are designed accord- 
ing to (5). For the 20° zig-zag, the length of each arm 
(2L) is chosen to be equal to 1 63/64” so that (5) is satis- 
fied at 3600 mc; 40 is then found to be 2”. For the other 
zig-zags the armlength 2L is taken to be equal to the 
proper value required by the corresponding angle. The 
ground plane in each case is of the dimension 2 by JX 
at the frequency 3600 mc. 

The radio telescope made up with the zig-zags ope- 
rates at 300 mc. But at 300 mc, the physical dimensions 
of the antenna become unwieldy for elaborate pattern 
studies. From the principle of model measurement? it is 
known that the radiation pattern of the antenna at 300 
mc will be the same as that at 3600 mc, provided that 
the physical dimensions of the antenna are scaled up by 
a factor of 12. The choice of 3600 mc as the frequency of 
measurement is found to be advantageous both electri- 
cally and mechanically. All the patterns are measured 
using an automatic antenna pattern plotter, where the 
required field is produced by a standard optimum horn 
and the antenna under test is used to receive that field. 
The receiving antenna is then rotated about the desired 
axis (choice of axis depends on the polarization of the 
field) and the received signal is then used to plot the 
directional pattern with the help of an automatic de- 
vice. 

The E- and H-plane radiation patterns at 3600 mc for 


the 20° antenna are shown in Fig 2 and Fig. 3. The theo-. 


retical patterns are superposed on these for comparison. 
The agreement between the two in the H plane is reason- 
able. The observed E-plane pattern shows a big sidelobe 
at an angle of 30°. The possible reason for this asym- 
metry is the feeding arrangement which causes an in- 
herent asymmetry of the antenna in the £ plane. 

It is found experimentally that the length d of the 
feed (see Fig. 1) has some noticeable effect on the pat- 
tern if it is greater than \/4. With the increase of d 
beyond )/4 the sidelobe at 30° (in the E plane) increases 
and the direction of maximum shifts away from @=0° 
direction towards the feed side. This is due to the fact 
that when d is large, it starts to radiate appreciably 
thereby effecting the end fire property of the antenna. 
The effect of d when it is less than \/4 is not observable 
experimentally. 


*G, Sinclair, “Theory of models of electroma ic sys ? 
, gnetic systems, 
PRoc. IRE, vol. 36, pp. 1364-1370; November, 1948. 
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Fig. 2—The radiation pattern of the 20° zig-zag antenna in 
the E plane. 
measured, — — — — — calculated 
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Fig. 3—The radiation pattern of the 20° zig-zag antenna in 
the H plane. 
‘ ——————— measured, — - —- — — calculated. 


From Fig. 2 and Fig. 3 it is found that considering the 
electrical length of the antenna (1.37A) it has a very 
good directivity. The directive gain of the antenna is 
found to be equal to 13.5 db. The cross-polarization 
component of the antenna in the H plane is found to be 
25 db down. From this it appears that the antenna may 
be used very advantageously for directional purposes; 
it seems to be better than the commonly used antennas 
for this purpose at the vhf and uhf ranges. 

The measured half-power beamwidths of the pat- 
terns and the corresponding sidelobe ratios for the dif- 
ferent antennas at 3600 mc are shown in Table I. 


TABLE I 


Pitch Angle Half-Power Beamwidth Sidelobe Ratio (db down) 


2a E plane H plane E plane H plane 
152 31° 36° 9.5 15 
20° 22° 28° 1055 i Ba 
ZS: 26° 34° 7.0 8.0 
30° 26° 34° 4.4 6.0 


EFFECT OF FREQUENCY ON PATTERN 


The radiation patterns of the 20° zig-zag antenna are 
measured throughout a range of frequencies approxi- 
mately 20 per cent above and below the center fre- 
quency 3600 mc, in both the £ and H planes. These are 
reproduced in Fig. 4 and Fig. 5. It is found that about 5 
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Fig. 4—Measured E-plane patterns of the 20° zig-zag antenna 
‘at different frequencies. 
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Fig. 5—Measured H-plane patterns of the 20° zig-zag antenna 
at different frequencies. 


per cent above and below the center frequency the radia- 
tion pattern stays stable maintaining the maximum in 
the end-fire direction. As the frequency is changed 
further, the main beam gradually broadens and the 
sidelobe increases. When the frequency is about 10 per 
cent off from 3600 mc, the main beam splits up creating 
a minimum at the 6=0° direction and two maxima on 
either side of it. This phenomenon is similar to the case 
of the helical antenna.4 From the terminology used in 


4 Kraus, op. cit., p. 211. 
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the helical case, we can say that the bandwidth of axial 
mode operation of the zig-zag, designed for 3600 mc, is 
400 mc or about +5 per cent of the center frequency. 


COMPARISON BETWEEN THE ZIG-ZAG 
AND THE YAGI 


It is worthwhile to make a comparative study of the 
Yagi and the zig-zag, because the former is an antenna 
very commonly used in the vhf range where the zig-zag 
can be used with great advantage. One conspicuous ad- 
vantage of the zig-zag is its less frequency sensitiveness 
which makes it a very promising antenna for television 
reception. 

The directivity of the zig-zag antenna is better than 
the Yagi. Alfred® reported that a Yagi at 300 mc with 8 
elements and a cylindrical reflector has a half-power 
beamwidth of 39 in the £ plane and 40 in the H plane. 
From the results reported above the properly designed 
20° zig-zag has much better directivity in both the prin- 
cipal planes. 

Alfred also found that the Yagi has a back to front 
ratio of 30 per cent, whereas the zig-zag has a 6 per cent 
back radiation even with a very small ground plane. 
However, the back radiation from the Yagi may be re- 
duced considerably by slightly detuning the elements, 
but this causes a reduction in gain. 


IMPEDANCE CHARACTERISTICS OF THE 
ZIG-ZAG ANTENNA 
The impedance of the full scale zig-zag (20°) is meas 


ured in the 300-mc band using the standard technique. 


5 R. V. Alfred, “Experiments with Yagi aerials at 600 mc/s,” 
J. IEE, vol. 93, pt. Illa, p. 1490; 1946. 
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The results of the measurements are shown in Table II. 


TABLE I] 


Frequency (mc) Input Impedance (ohms) 


290 150—j40 
295 135—j80 
300 160—j80 
305 110—j785 
310 125—j785 
315 152—j0 

320 170—j80° 
325 135—7120 
330 70—795 


CONCLUSION 


From the above report it is evident that the zig-zag 
is a good directive antenna with low sidelobe and can 
be used in the vhf and uhf ranges. We have taken the 
length of the antenna as 1.37) so that the antenna may 
not be large physically. By taking more V’s the direc- 
tivity of the antenna may be increased further. Although 
this is a special case of the helical antenna, it differs 
from the latter in being sensitive to linearly polarized 
waves. 
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A Note on the Effective Aperture of Electrically 
Scanned Arrays* 
R. W. BICKMOREt 


Summary—Since the size of current radar antennas must be in- 
creased in order to detect aircraft and missiles flying at supersonic 
speeds, it is no longer possible to use rapid mechanical scanning. 
Therefore this type of scanning is being gradually replaced by elec- 
trical scanning. 

The popular approximation which states that the effective aper- 
ture length is proportional to cos 0, where @ is the scan angle from 
broadside, is examined and shown to be adequate in many cases 
where the end-fire effect is negligible. An exact effective aperture 
based on half-power beamwidth is derived and compared with the 
approximate solution for various aperture lengths as a function of 
scan angle. 


* Original manuscript received by the PGAP, July 31, 1957; re- 
vised manuscript received, December 19, 1957. 
+ Res. Labs., Hughes Aircraft Co., Culver City, Calif. 


HE popular approximation which states that the 

effective aperture length is proportional to cos 6, 

where 49 is the scan angle from broadside, has 
come to be regarded by many as an exact relation. For 
many arrays, this approximation can be shown to be 
completely adequate. However, the increased use of 
electrical scanning as a replacement for mechanical 
scanning makes a re-evaluation of the effective aperture 
particularly important. 

For the purposes of this paper, consideration of a 
plane square aperture carrying a uniform, linearly polar- 
ized, continuous current will be most convenient. As 
shown in Fig. 1, the aperture lies in the xy-plane with 
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Fig. 1—Far zone geometry for a square array. 


current J, flowing from —a/2 to a/2. The phase dis- 
tribution is assumed to be separable and linear. 

Upon evaluation of the far-field components in the 
usual manner,! the following expressions are obtained in 
the three principal planes: 


F,(0) sin 1 
; um (1) 
@= 7/2 


_ where u,=ka/2 (sin 8—sin 001), 80.=direction of main 
_ beam in the yz-plane. 


sin Us 


E,4(6) | ~ cos 6 (2) 


us a 
¢=0 
where u.=ka/2 (sin @—sin 692), Qo. =direction of main 
beam in the xz-plane. 


sin 21 Sin v2 


E3(¢) | ~ sin ¢ (3) 


V1 V9 
6= 7/2 


where v;=ka/2 (cos ¢—cos $0), d0=direction of main 
beam in the xy-plane, v.=ka/2 (sin @—sin go). 

The element factor is immediately seen to be unity in 
(1), cos @ in (2), and sin ¢ in (3). 

' The array factor is seen to be the same function of u 
' regardless of scan angle in the two elevation planes. 
However, neither ~ nor sin~! w are linearly related to 
physical space except when 6)=0 (1.e., broadside radia- 


a’ 0.443 
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sin 6 — sin 6 
= sin (6 — 4) cos 6) — [1 — cos (@ — )] sin 4. (4) 


The quantity cos 4) is thus factorable when the second 
term of (4) is negligible. This would appear to be true 
only when 4 is small, which corresponds to radiation 
near broadside. 

It can be noted that the approximation completely 
fails to predict the end-fire phenomenon of the array 
factor. That is, the approximation indicates an omni- 
directional pattern when 0)=7/2, whereas the more 
exact expression gives 


ka 
sin E (1 — sin 0) | 
Ew = (S) 


ka a in 6) 
— (1 — sir 
5 nN 


which is far from omnidirectional. 

The approximate solution also fails to predict the lack 
of symmetry of the radiation pattern. The angular dis- 
tances of each half-power point from the peak of the 
main beam are unequal. 

If these angles are designated A and 6, then the total 
half-power width is 


IN 
© = A+6=rcsin [0.443 — + sin | 
a 


a 


nN 
+ arc sin [0.443 — — sin 0] (6) 


where A is the semi-half-power width for 6>6 . 6 is the 
semi half-power width for 0 <4. 

This expression for 0 may now be used to determine 
an “effective aperture” by equating it to the half-power 
width of a broadside array and solving for the required 
broadside aperture length. This procedure circumvents 
the mathematical difficulty of equating the expressions 
for gain, and still allows for obtaining an “effective 
gain” by using the expression G,/G2= [02/0,]?, which 
is a good approximation for arrays larger than about 
10A in length. 

If a’ is the effective aperture length, 


(7) 


nN rN 
sin }1/2 E sin (0.48 — + sin on) + arc sin (0.143 — — sin a} 
a / 


or 


a 


a’ a f r\2 1/2 a 1/2 
—- = —— (1 + 0.443 “) — sin? | -b I¢ — 0.443 *) — sin? | \ : (8) 
ey in : : 


tion). In order that the pattern remain unchanged with 
scan angle 60, the quantity (sin 6—sin 49) would have 
to be identical to the quantity sin (@—6)). However, it 
is easily shown that 


1S. Silver, “Microwave Antenna Theory and Design,” McGraw- 
Hill Co., Inc., New York, N. Y.; 1949. 


The approximation, of course, states that 


a a a 
— ~ — = — cos I. (9) 


Figs. 2(a) and 2(b) show the ratio of these two values 
of effective aperture plotted against the scan angle 0) 
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Fig. 2—Ratio of exact to approximate effective aperture. 


for a/\=10 and a/A=100. A unidirectional array is as- 
sumed, in these computations, and when the A semi- 
half-power point passes into the region 6>7/2, it be- 
comes a point in imaginary space. The curves in Fig. 
2(a) and Fig. 2(b) then give an exact solution to the 
problem of an array above an infinite ground plane. The 
portion of the beam in imaginary space gives rise to an 
imaginary component of effective aperture, and the 
associated energy appears as stored energy in the case 
of an infinite ground plane or becomes available to be 
diffracted or absorbed in the case of a finite ground 
plane. 

Returning to the array factor (3), it can be seen to 
differ from (1) and (2) in that it is the product of two 
complementary array factors. The general effect is that 
of sharpening the pattern in the plane perpendicular to 
the plane of scan. However, for fairly large arrays, the 
sharpening is quite small. This can be seen by observing 
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the pattern in the vicinity of ¢.=0 or go=7/2, where 
one factor will be found to almost completely dominate 
the other. Thus, except for short arrays, the pattern on 
a surface orthogonal to the scan plane remains essen- 
tially unaffected. 

However, the broadening in the plane of scan is not 
negligible and becomes more pronounced as the array 
length increases. This can be illustrated by comparing 
the half-power beamwidths of the array factors for 
6)=0 and 0)=7/2. The quantity sin «/u has amaximum 
value of unity for ~=0 and becomes 0.707 for u=1.39. 
Therefore, 


N nN 7 
@| = 2sin | 0.443 ~| ~ 0.88— fora>d, (10) 
a a vei F 
6. = 0 
and 
nN N 
@| = 2cos-t E ~ 0.443 “] ss 1.88 4/> 
a a 


fora>> . (11) 


00 = a/2 


Thus the ratio of the full end-fire beamwidth to broad- 
side beamwidth becomes 2.14/a/r in the plane of scan. 
As Jordan? points out, however, a comparison of first- 
null widths reveals that the end-fire pattern is only 
+/2a/x times as broad. 

It is readily seen that in order to avoid this beam de- 
terioration, the individual elements comprising a uni- 
directional square aperture that is @ wavelengths on a 
side must be progressively rearranged, as the scan angle 
is increased, to form a rectangle of increasingly greater 
aspect ratio until, at end-fire, the elements form a line 
source 4.57 (a/\)? wavelengths long.* The physical un- 
desirability of such an arrangement illustrates the 
major disadvantage suffered by an electrically scanned 
planar array compared with a cylindrical surface array 
for one-dimensional scanning or a spherical surface ar- 
ray for two-dimensional scanning. 
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? E. C. Jordan, “Electromagnetic Waves and Radiating Systems,” 
Prentice-Hall, Inc., New York, N. Y.; 1950. 
’ Assuming no supergaining in the purely end-fire direction. 
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The Characteristic Impedance of Two Infinite Cones of 
Arbitrary Cross Section* 
ROBERT L. CARREL} 


Summary—A theoretical method of finding the characteristic 
impedance of two infinite cones of arbitrary cross section is pre- 
sented. In this problem, the spherical geometry can be mapped on a 
plane, thus reducing the problem to finding the solution of Laplace’s 
equation by using conformal mapping techniques. This method is ap- 
plied to the inclined biconical antenna and several types of fin an- 


_tennas. The solution for the biconical antenna is Z)=/27 cosh say ob 


where F is a function of the geometry of the antenna. For the fin, the 
solution takes the form Z)=K(k)/K’(k) where K and K’ are com- 
plete elliptic integrals of the first kind and k depends on the ge- 
ometry. 


INTRODUCTION 


HE problem is to find the characteristic imped- 

ance, Zo of an antenna consisting of two infinitely 

long conical conductors which have a common 
apex. The following method is applicable to all cone- 
shaped surfaces of arbitrary cross section, 7.e., surfaces 
which can be defined by the two spherical angles 6 and 
@. Basically, the problem consists of solving Maxwell’s 
equations subject to the condition that tangential E 
vanishes on the surface of the conductors. For the TEM 
spherical wave it can be shown that Maxwell’s equa- 
tions reduce to Laplace’s equation in two dimensions. 
Conjugate function theory is then applied in solving 
this two-dimensional potential problem. 


FORMULATION OF THE PROBLEM 


A TEM spherical wave is defined as having no radial 
components of electric and magnetic field F and H. 
With the usual e*' time convention, it follows from 
Maxwell’s equations VX H=iweE and VX E= —iwuH 
that a TEM spherical wave is represented by H =V X?g 
and E=VXH/iwe where ?=unit vector in the radial 
direction and g=e~**"T(6, @) is a scalar function of the 
spherical coordinates 7, 6, and ¢ such that 8? =w7ye and 
V?T =0. Such a wave can be propagated along two in- 
clined conical surfaces if the function T is constant on 
the surfaces of the conductors and if V?7’‘=0 in the sur- 
rounding space. 

We seek a change of variables that will transform the 
spherical coordinates @ and ¢ into a plane in such a way 
that the identity of the conical boundaries is preserved 
and that Laplace’s equation remains unchanged. 

Bet 


p = f(9) 
o- = > 
V(o, Pc) F T(6, ¢) e 


(1) 


* Manuscript received by the PGAP, September 3, 1957. This 
work was performed under the sponsorship of the U.S. Air Force, 
Wright Air Dev. Center, Wright-Patterson Air Force Base, Ohio. 

+ Antenna Lab., University of Illinois, Urbana, Ill. 


(In this and the following development p, ¢,, and z are 
cylindrical coordinates; 7, 6, and ¢ are spherical coor- 
dinates). If 0V/dz=0T/dr =0 (the initial premise of the 
theory), Laplace’s equation in the two coordinate 
systems is 


a/ ov 1 eV ; 
a (s ae Stag dee @) 
and 
“(sing ~-)+ eee 0 (3) 
00 00 sin 0 d¢? 


After transforming the variables as dictated by (1), (2) 


becomes 
a gi a LOE i 
aol f’(@) 90 f(0) ag? 
where 
vo . 4] 
We ETE 
(3) and (4) will be identical if 
f’(9)/f(8) = csc 8. (5) 


After integration we find that f(@)=c tan 0/2, where c 
is an arbitrary constant of integration which we may set 
equal to unity. The required mapping relation is then 


p = tan ee 
b= ¢ ] 


The transformation maps the surface of any sphere 
r=constant into the (p, ¢.) plane. Further, the upper 
hemispherical surface maps onto the interior of the 
unit circle; the lower hemispherical surface maps onto 
the exterior of the unit circle. Fig. 1 shows the mapping 
of a cone of arbitrary cross section. It can be seen that 
the radial lines ¢.=constant describe the spherical co- 
ordinate ¢, and the concentric circles p=tan 0/2 corre- 
spond to the spherical coordinate @. Since the trans- 
formation (6) also preserves Laplace’s equation, all 
that remains is to find a solution for the cylindrical 
problem in the (p, ¢-) plane. Thus the dynamic field 
problem has been reduced to a static field problem, the 
solution of which yields the capacity per unit length. 
However, we are primarily interested in the character- 
istic impedance, a dynamic field quantity, which is re- 
lated to the capacity for the case of a two conductor 


(6) 


1 For a similar result see W. R. Smythe, “Static and Dynamic 
Electricity,” McGraw-Hill Book Co., Inc., New York, N. Y., pp. 
243-245; 1950. 
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p= ton 0/2 
Sad 


(a) (b) 


Fig. 1—The transformation of a general conical section onto a 
plane. (a) Sperical coordinates, (b) Cylindrical coordinates. 


system operating in the principal mode. Specifically, 
Zo=en/C, where e is the capacitivity, 7 is the intrinsic 
impedance of the medium, and C is the capacity per unit 
length. 

Let us consider some examples involving different 
geometrical shapes and orientations. 


THE COAXIAL BICONICAL ANTENNA 


The coaxial biconical antenna affords a straightfor- 
ward illustration of the technique. The structure is posi- 
tioned in the spherical geometry as shown in Fig. 2(a), 
and yields two concentric circles in the (p, ¢-) plane. 
This is the familiar coaxial line problem for which the 
solution is known. 


Zo = n/2n log b/a, (7) 


where 7 is the intrinsic impedance of the medium be- 
tween the conductors. Substituting for 6 and a gives 


Zo = n/ 2m log cot ¥i/2 cot 2/2. (8) 


This is the exact solution given by Schelkunoff.’ 


THE NONCOAXIAL BICONICAL ANTENNA? 


This structure is positioned in the spherical coordi- 
nate system so that the singular point 0 = 7 is exterior to 
both cones. The plane of the axes of the cones is the 
o¢=7/2, 37/2 plane, as shown in Fig. 3(a). 

The equation of a cone in the spherical coordinate 
system is 


sin 0’ cos 9’ sin 8 cos ¢ + sin @’ sin ¢’ sin 6 sin } 
+ cos # cos@= cosy, (9) 


where y is the cone half angle, and 6’ and @’ describe 


the axis of the cone. Transforming variables is dictated 
by (6), we find that 


2p sin 6’ cos (¢- — ¢’) + (1 — p?) cos 6’ 
= (1+ p”) cosy. (10) 
If we make another change of variables from cylindrical 


coordinates (p, ¢-) to cartesian coordinates (u, v) such 


2S. A. Schelkunoff, “Electromagnetic Waves,” D. Van Nostrand 
Co., Inc., New York, N. Y., p. 287; 1943. 

8S. A. Schelkunoff, “Advanced Antenna Theory,” John Wiley 
and Sons, Inc., New York, N. Y., p. 27; 1952. 


a = tony,/2 
b = ton (w—p/ 


(a) (b) 


Fig. 2—The transformation of two coaxial cones onto the w plane. 


(a) Spherical coordinates, (b) w=-+ iv planes. 


Vv 


(b) 


Fig. 3—The transformation of two noncoaxial cones onto the 
w plane. (a) Spherical coordinates, (b) w plane. 


that p?=u?+v? and ¢.= tan —1v/u, we find that (10) is 
the equation of a circle. After some algebraic manipula- 
tion we obtain 


( sin @’ cos ¢’ ) ( "sin 6’ sin ¢’ ) 
“ ———_———_) + (» —- —_—_—_ 
cos 6’ + cosy cos 6’ + cosy 


sin y id 
= (—=*_) ~ ee 
cos 6’ + cosy 


Thus the two cones in Fig. 3(a) map into the two circles 
in the w=u-+iv plane, as shown in Fig. 3(b). 


This problem corresponds to finding the character- 


istic impedance of a uniform two-wire transmission line. 
Its well-known solution is* 


Zo = n/ 2m cosh“! E =f ej ac et ae =] 
or n/ : 


2rire 


(12) 


Upon substituting 


sin 6,’ sin 6.’ 
Ui vo = = 
cos 6;/ + cos Wy cos 65’ + cos P2 
sin Wy sin 2 
i BOY ea 
cos 6;' + cos Wj cos de’ + cos pz 


and noting that @;’ +62’ =@o where @) is the angle between 
the axes of the cones, we find that 


cos ¥1 Cc = 0 
Zo = n/ 2a cosh | id taeda | (13) 


IL sin y; sin W. 


*W. R. Smythe, “Static and Dynamic Electricity,” McGraw- 


Hill Book Co., Inc., New York, N. Y., pp. 80-82; 1950. 


sf 


_ Fig. 4—The characteristic impedance of a biconical antenna with 


2 


7] 


oriented to exclude the singular point 6 =7, and lie in the 


steps and will be presented later. It is seen that under 


: 


+ 


Nn 


Nn 


Characteristic Impedance, Z,. —ohms 
uw 


equal cone angles. 4 is the angle between the axes of the cones, 
y is the cone half angle. 


It can be shown that (13) reduces to (8) when 0)=7. 
Fig. 4 is a graph of the characteristic impedance of a 
biconical antenna with equal cones showing the de- 
pendence on 4». 


THE COPLANAR FIN ANTENNA 


The coplanar fin antenna is shown in Fig. 5(a). It 
consists of two infinitesimally thin conductors which 
have a common apex and lie in a plane. The fins are 


¢=0, z plane. 


It is instructive to consider first the special case of 


two equal fins,° z.e., Yi=Yo=y. Let 6 =02 =0)/2 where 
0. is the angle between the axes of the fins. The extension 


to the more general case involves several intermediate 


the transformation of variables dictated by (6) the two 
fins map onto two equal slits in the w=u-+7 plane. The 


mapping of the w plane onto the a plane by the function 


Rep. iby. R:-E. Carrel; 
Antenna of Infinite Length,” Tech. Rep. 
3220, January 15, 1957, and is included here for completeness. 


o = cot E (= = ¥)| w (14) 
Je ee 
is shown in Fig. 6(a). 
The Schwarz-Christoffel transformation, 
pe iF of (15) 
0 Vl =9?)(1 — ko”) 


maps the o plane onto the interior of the rectangle 
P,P; P2P2 in the ¢ plane as shown in Fig. 6(b). With the 
help of (15), both K and K’, as defined in Fig. 6(b) can 


5 This topic is the subject of a University of Illinois Antenna Lab. 
“The Characteristic Impedance of the Fin 
No. 16, Contract AF33(616)- 
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970 
w=—ton 1/2 (6,' + 2) 
Wa=—ton V2 (Be'— ye) 
Ww," ton/2(6,'—y,) 
Wa* ton 1/2(6,' +y,) 


(a) 


(b) 


Fig. 5—The transformation of two coplanar fins onto two slits 
in the w plane. (a) Spherical coordinates, (b) w plane. 


(a) (b) 


Fig. 6—Mapping the slit « plane onto a period rectangle in the 
¢ plane. (a) Complex o plane, (b) Complex ¢ plane. 


be expressed in terms of k. We have 


K= fs a (16) 
Jo JA — AG — #P) 
and 
1k ds 
ae 
. ih V(1 = s)(1 — Bs?) 


The latter integral can be brought into a more elegant 
form if we make the substitution s = (1 —k’t?)—/? where 

ki’? + hk? = 1. (17) 
Thus 


fo (18) 


iF di A 
awl = 3) ee RP) 


Notice that (16) and (18) are complete elliptic integrals 
of the first kind of modulus k and Rk’, respectively. The 
configuration in the ¢ plane can be considered as a 
parallel plate transmission line, which has a character- 
istic impedance given by 


Zo = nd/b, (19) 


where 7 is the intrinsic impedance of the medium be- 
tween the two conductors, d is the distance between 
plates, and 0 is the width of the plates. Note that this 
is the exact solution; there are no fringing effects, be- 
cause the total electric field in the o plane is mapped into 
the interior of the period rectangle in the ¢ plane. Thus 


the characteristic impedance of the fin is given by 
Zo = nK/K’, (20) 


where 
sin 09/2 — siny 


Locka rea ie | 
oS sin 69/2 + siny 


(21) 
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Fin Half-Angle, » — degrees 


Fig. 7—The characteristic impedance of the coplanar fin antenna 
~ with equal fin angles. y is the fin half angle, 0 is the angle between 
the axes of the fins. 


A graph of Z» vs w for various values of the parameter 
6) is shown in Fig. 7. 

In the case of arbitrary fin angles, as shown in Fig. 
5(a), a slightly more complicated mapping function is 
needed to obtain the symmetrical configuration in the 
o plane. As shown in Fig. 5(b), the configuration in the 
w=u-+iv plane exhibits no symmetry about the v axis. 
The general bilinear transformation may be employed 
to obtain the desired mapping. This transformation has 
the well-known property of mapping any three points in 
the w plane into any three arbitrarily chosen points in 
the o plane. Let us choose this correspondence between 
points as follows: 


0) > oy, = — 1/k 
W.>02.> —- 1 
W3> 03 = + 1 (22) 
This transformation takes the form: 
We — Wz We— WwW, o +1/k 2 
; ss : p (23) 
OU nL ial tO ao—1l 1— 1/k 


Note that k is as yet undetermined. To fix its value we 
can require w=w, when ¢=0,=1/k. Making this sub- 
stitution in (23) along with the values of w,, we, w3, and 
ws from Fig. 5(b), and solving for k, we find 


hho) 
pr et hee /ance (24) 
Cc Cc 


where 
WA Ww Ws 


A eet ee fy 1 cos % — cos (Wi + 2) 


Ws — Wz We — Wy 2 


sin y, sin ps oe 
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(a) (b) 


Fig. 8—The transformation of equal noncoplanar fins into arcs 
of circles in the w plane. (a) Spherical coordinates, (b) Complex _ 
w plane. 


and 6) = 6,’ +6.’ is the angle subtended by the axes of the 
fins. Thus (23) maps the w plane of Fig. 5(b) into the « 
plane of Fig. 6(a). From this configuration, we can pro- 
ceed directly to the solution by the method outlined for 
the case of two identical fins described earlier, with the 
exception that the modulus k is given by (24) and (25). 
Therefore, 


Zo = nK/K', (26), 


where K and K’ are the complete elliptic integrals of 
the first kind. 


THE NONCOPLANAR FIN ANTENNA 


This section considers an antenna consisting of two 
equal fins, which are not in the same plane, and have a 
common apex such that the fins are normal to the plane 
passing through the center line of each fin, as shown in 
Fig. 8(a). To resolve this problem, it is necessary to de- 
termine how (6) transforms planes through the 0@=7/2, 
o=7/2 axis. For any value of 6 (80/2 is the inclination 
of the plane to 6=0) such that 0<6)<7, there are two 
possible planes passing through the 0=7/2, ¢=7/2 
axis. The equation of these planes in spherical coordi- 
nates is 

80 
Se reat WEEN bat bh (27) 
If the transformation (6) is followed by a transformation 
to cartesian coordinates in the w plane, we find that 
(27) becomes 


(a 


which is the equation of a family of circles centered on 
the u axis passing through the points v= +1. 

It is also necessary to find the locus of the edges of a 
fin of constant half angle y as the 6» angle is varied. It 
can be seen that this locus describes two cones whose 


t ~) + 2 » bo (28) 
(31 0} UP = CSC. = 
2 2 


axes are the lines 


I 


fe T $ 1. d 
=— ===) an 0 
2 2 
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(b) 


Fig. 9—Maps of the fins in the complex £ and ¢ planes. (a) 
A Complex £ plane, (b) Complex ¢ plane. 


; The equation of these cones in spherical coordinates is 


+ sin @ sin ¢ = sin y. (29) 
; By using transformation (6) followed by a transforma- 
5 tion to cartesian coordinates in the w plane, (29) be- 
E comes 


- 


u? = (v + csc py)? = cot? y, (30) 
_ which is the equation of a family of circles centered on 
_the v axis. The representation of a typical set of fins, 
_ along with the above family of circles, is shown in Fig. 
— 8(b). 

In order to map the fins BC and EF of Fig. 8(b) into 
_a more amenable geometry, we apply a bilinear trans- 
_ formation followed by a power transformation. 


t+w 
pene at) 
o = (Ré)*/% (32) 


where 


aa (- + ae 
1 — siny 


‘The bilinear transformation (31) maps the circles cen- 
tered on the uw axis in the w plane in Fig. 8(b) into the 
radial lines in the & plane; it also maps the circles cen- 
tered on the v axis onto concentric circles, as shown in 
Fig. 9(a). Thus the fins BC and EF map into the inclined 
slits as shown. The power transformation (32), con- 
taining the normalizing factor R, yields the colinear 
slits, as in Fig. 9(b). This corresponds to Fig. 6(a), aside 
from a rotation. The coordinates of point £ in Fig. 8(b) 
are determined by the intersection of the two circles G 
and J. The value of k can be found by substituting the 
coordinates of point £ into the transformations (31) and 


(32). Thus, 
( 1 — sin v\7/% 
feet) 
1+ siny 
The solution for the geometry of Fig. 9(b) has been ob- 
tained in a previous problem. It is given by 


0<.00 3 x. (33) 
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Fig. 10—The characteristic impedance of the noncoplanar fin 
antenna with equal fin angles. y is the fin half angle, 00 is the 
angle between the fins. 


Zo = nK/K’ 


where k, the modulus of the complete elliptic function 
K, is given by (33), Fig. 10 is a graph of Zo vs W for 
various values of 4. 

The extension of these results to arbitrary fin orienta- 
tions is straightforward ; however, the results appear in a 
cumbersome form, so no attempt will be made to pre- 
sent them here. 
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Microwave Stepped-Index Luneberg Lenses* 
G. D. M. PEELER} anp H. P. COLEMAN 


Summary—Since materials with a continuous variation of index 
of refraction were unavailable for the construction of spherical Lune- 
berg lenses, an investigation was made of stepped-index lenses in 
which the desired continuous variation of index with radius was 
approximated by a number of constant-index spherical shells. Ex- 
perimental data on two-dimensional, TE, mode, stepped-index 
Luneberg lenses indicated that 10 shells were sufficient for 18-inch 
diameter lenses for use at X band. Expanded polystyrene materials, 
developed with dielectric constants ranging from 1.1 to 2.0 in incre- 
ments of 0.1, were used to construct 10-step, 18-inch diameter, 
spherical Luneberg lenses. Performance data for these lenses are 
given at wavelengths of 3.2, 1.8, and 0.8 centimeters. 


INTRODUCTION 


N RECENT YEARS, the Luneberg lens has come to 
if be considered by many workers in microwave optics 

an ideal wide-angle objective due to its complete 
symmetry. The lack of suitable materials has, however, 
made the construction of spherical lenses difficult. As a 
result, most Luneberg lenses constructed have been 
two-dimensional models! which utilize waveguide or 
geodesic analog techniques to achieve the well-known 
variation of index of refraction m with radius r 


w= /2 — 7. 


These techniques are not adaptable to the construction 
of spherical lenses; the construction of these lenses is 
consequently dependent upon the development of ma- 
terials with a suitable range of refractive indices. This 
paper describes the development of a lens, utilizing ex- 
panded plastic materials. The variation in n is obtained 
by varying the density of the materials. 

It would be desirable, of course, to obtain dielectric 
materials which could produce a continuous variation 
of with position and to construct a lens from these. 
The limitation of available techniques, however, pro- 
hibits this, and consequently it was decided that a 
stepped variation approximating the continuous varia- 
tion of m should be investigated as a means of lens con- 
struction. The first part of the paper discusses the choice 
of stepped positions and index increments, and considers 
certain stepped-n, two-dimensional lenses constructed 
for the purpose of determining the effects of stepping. 
The latter portions of the paper are devoted to the per- 
formance of an 18-inch stepped-spherical lens. 


* Manuscript received by the PGAP, February 8, 1957. 

} Missile Systems Div., Raytheon Mfg. Co., Bedford, Mass. 
Formerly with the U. S. Naval Res. Lab. 

t Microwave Antennas & Components Branch, U. S. Naval 
Res. Lab., Washington, D. C. 

1G. D. M. Peeler, and D. H. Archer, “A two-dimensional micro- 
wave luneberg lens,” IRE TRANS. ON ANTENNAS AND PROPAGATION, 
vol. AP-1, pp. 12-23, July, 1953, and references therein. 


STEPPING TECHNIQUES 


Some conclusions about stepped-index lenses may be 
reached through general considerations. For example, 
if a lens were constructed with a large number of steps 
and with small index increments, it would be expected 
that ray paths through the lens would be only slightly 
perturbed from those in the ideal continuous lens and, 
as a consequence, little deterioration in performance 
would result. On the other hand, for a lens with few 
steps, especially if the distance between steps were 
large compared to a wavelength, ray paths would pass 
through only a few large constant-index regions and 
only slight focusing could be expected. Thus, there 
probably exists a minimum number of steps for a lens 
to provide some specified performance. A choice may 
be made among several distinct methods for determin- 
ing the radii of step positions, and the constant index 


—— 


values between steps so as to approximate the desired - 


continuous index vs radius variation. One method is to 
choose equal index increments, 


Nm — Mm+1 = An = constant, 


and to choose the step positions at the radii of mean 
index, (%m+%m4i1)/2, in the ideal lens. Here m is an 
integer denoting the particular shell. Another method 
is to choose equal dielectric constant steps 


— €n41 = Ae = constant, 


Em 


where €=n? and to choose step positions at the radii of 
mean ¢ in the ideal lens. Still another method is to 
choose increments in either 2 or e so that the ratio of 
successive step values is constant 

Nm € 


mm 
= constant or —— = 
Em+1 


constant 
Mm+t 
and then to choose step positions accordingly. 
Originally it was hoped that phase fronts could be 
calculated for stepped-index lenses, with a range of steps 
for each of the stepping methods, in order to provide a 
basis for determining the minimum necessary number of 
steps and the optimum stepping method for designing 
a Jens with specified performance. However, attempts 
to set up valid phase-front calculations were unsuccess- 
ful. Calculations were performed by introducing the 
approximation that ray paths are straight lines between 
lens steps and that Snell’s law is satisfied at the steps. 
It was evident from these calculations, however, that 
this approximation was too gross to provide any useful 
information. 
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Two-DIMENSIONAL STEPPED- 
INDEX LUNEBERG LENSES 


The construction and measurement of two-dimen- 
sional stepped-index lenses was chosen as the most ad- 
vantageous method of obtaining guides for a stepped- 
index spherical lens design. The equal ¢€ increment 
method of stepping was chosen since ¢€ is usually more 
directly measured than n. In a continuously variable 
Luneberg lens € varies with radius r according to the 
formula e=n?=2—r?, where 0<r<1. Thus, assuming 
equal € increments for a lens with p steps, Ae=1/p, and 
the dielectric constant €,, for the mth region is 


m 
€m = 2——> 


p 


where m=O for the center region. The step positions 
were chosen at radii for which e=2—(m-+1/2)Ae in the 
continuously variable lens, where m=0, 1,---, p—1. 
In this arrangement e deviations from 2—,7? are maxi- 
mum at the steps and equal in magnitude on each side 
of the step. The boundary radii for the mth region are 


oe 
T= . 
p 


Two lenses were constructed, one with p=10 and 
the other with p=20. The polystyrene-filled-waveguide 
technique, previously used for constructing a two-di- 
mensional, continuously variable index lens,? was chosen 
for constructing the stepped-index lenses. The desired 
equivalent dielectric constants were obtained by vary- 
ing the plate spacing stepwise, with constant plate 
spacing between steps. For a waveguide filled with a 
material of dielectric constant e’ propagating the TEto 
mode, the equivalent dielectric constant is 


oN 2 
alae eee ee 
where } is the wavelength and a is the H-plane wave- 


guide dimension. The plate spacing for the mth region 
is obtained by combining (1) and (3) to give 


(1) 


(2) 


(3) 


r 


——————— 4 
2 RET] b) @) 
where m=0, 1,:--, p. (There are p+1 regions.) 

Fig. 1 is a diagram of the construction of the two- 
dimensional lenses designed to operate at \=3.2 cm. 
Polystyrene (e’ =2.503) was used in their manufacture. 
The tolerance on thickness of the dielectric slab and on 
depth of steps in the aluminum coverplates was within 
+0.002 inch. The actual plate-spacing tolerance could 


2 Ibid. 
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not be determined; however, both the polystyrene slab 
and the aluminum cover plates were stress-relieved be- 
fore machining in an attempt to minimize warpage and 
to maintain plate-spacing tolerances near the machining 
tolerances. If a plate-spacing tolerance of +0.004 inch 
is assumed, the error in desired equivalent ¢ is +0.005 
at the lens center and +0.023 at the lens edge. Each lens 
was illuminated with the polystyrene-loaded-waveguide 
horn utilized for the continuously variable lens.* 


Lue Z Z 
POLYSTYRENE 


Fig. 1—Diagram for two-dimensional lens construction. 


The phase front of each lens in the plane of the lens 
and along a tangential aperture, was measured using 
standard techniques.‘ Typical sets of phase front data 
for these two 18-inch-diameter lenses are shown in Fig. 
2. For the 10-step lens, deviations of the phase front 
from a linear phase front are approximately twice the 
deviations for the 20-step lens, although the deviations 
in each case are less than \/10. 

Typical E-plane radiation patterns (in the lens plane) 
for the two lenses are shown in Fig. 3. The half-power 
beamwidths of both patterns are 4.5°; the beamwidth 
factor, defined as beamwidth X Jens diameter/wave- 
length, is 64°. Side lobes are down from peak power 20 
db for the 10-step lens and down 18 db for the 20-step 
lens. These compare favorably with the calculated 
sidelobe level of 17.5 db from a continuously variable 
index lens. Although the phase front data indicate that 
the sidelobe level should be slightly larger for the 10- 
step lens than for the 20-step lens, the phase errors are 
small enough that such comparison has little meaning. 


3 Peeler and Archer, op. cit., Fig. 2. : 
4 S. iver Ed., “Microwave Antenna Theory and Design,” 
Rad. Lab. Ser., Vol. 12, sec. 15:11; 1949. 
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Fig. 2—Typical phase fronts for two-dimensional lenses. 


12 6 O DEGREES 6 lj2 
Fig. 3—Typical E-plane radiation patterns for two-dimensional 
lenses. 


The above data indicate that for this aperture-to- 
wavelength ratio, ten steps are sufficient and that the 
method of stepping is adequate. The data also show 
that if continuously variable index materials are un- 
obtainable, stepped spherical Lunebergs can be con- 
structed. For other aperture-to-wavelength ratios, the 
phase-front data may be scaled to provide some indica- 
tion of lens performance for this method of stepping. 
For 1° beams which require an aperture-to-wavelength 
ratio 4.5 times as large as for these lenses, a 20-step lens 
is indicated if phase errors are to be kept below )/10. 
For beamwidths of one degree or less, the method of 
stepping may be optimized to decrease the number of 
steps necessary for a given maximum phase error. 


SPHERICAL STEPPED-INDEX LENSES 


It was decided to attempt construction of a stepped- 
index spherical Luneberg lens utilizing the technique of 
nesting together molded hemispherical shells of iso- 
tropic, homogeneous materials in which the index varia- 
tion from shell to shell is obtained by a density varia- 
tion. Since this lens should have no bandwidth limita- 
tions other than those imposed by the steps and aper- 
ture-to-wavelength ratio, it was desirable to choose an 
aperture size and a number of steps such that data at 
X band and higher frequencies would indicate limitations 
imposed by steps better than did the two-dimensional 
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lens data. Based on the two-dimensional lens data, 10-— 
equal-e steps and an 18-inch diameter were chosen. 

Since materials with the required dielectric constants 
were not available, it was necessary to develop suitable 
materials. It was not expected that the materials de- 
veloped for the first lens would meet all the requirements 
for use in general weather conditions and over an ex- 
treme range of temperatures; the goal was to obtain 
materials suitable for a laboratory lens with as many 
desirable properties as possible incorporated without 
undue expenditure of time or funds. The specifications 
on materials were as follows: specified dielectric con- 
stants maintained within +0.02; loss tangent 0.001 or 
less; water absorption less than 0.1 lb per square foot of 
exposed skin surface; stable both mechanically and 
electrically over the temperature range —20° to +150° 
F; permit moderate handling without cracking or dent- 
ing. The specified tolerance on dielectric constant is that 
obtained with the two-dimensional lenses. The loss 
tangent specification provides for total lens losses less 
than 0.5 db at \=3.2 cm or less than 2.0 db at 8 mm. 
The water absorption and temperature specifications 
insure against variations in dielectric constant and loss 
tangent under ambient weather conditions. 


Material Development and Lens Manufacture 


It was recognized at the beginning of this project 
that many expanded plastic materials would meet most 
of the specifications, but most of these materials have 
loss tangents greater than 0.001 at densities which pro- 
vide a dielectric constant of 2.0. Foam-in-place mate- 
rials, in general, have density variations from point to 
point too large to satisfy the tolerance on dielectric 
constant. It was decided to use polystyrene as the base 
material since it meets the loss tangent specification 
even in solid form and its mechanical properties are 
reasonably good even at low densities. However, poly- 
styrene materials of adjusted density had to be devel- 
oped which were amenable to molding into spherical 
shells with accurate radii. A contract was consequently 
let to Emerson and Cuming, Inc., Canton, Mass. for 
appropriate material development and lens construc- 
tion’ which was jointly sponsored by Air Force Cam- 
bridge Research Center and Naval Research Lab- 
oratory. 

Several methods for producing polystyrene foams 
were considered and preliminary experiments were run. 
The method which appeared to be most promising was 
to prepare a free-flowing polystyrene powder which 
upon compacting would have a predetermined bulk 
density. When heated above the softening point of poly- 
styrene (about 85°C), the discrete particles of powder 
weld together into a rigid material. By carefully select- 


°G. D. M. Peeler, e¢ al., “Microwave Stepped-Index Luneber 
Lenses,” NRL Rep. No. 4843; October, 1956. sf 
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Fig. 4—Dielectric constant vs density for expanded polystyrene. 
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‘ing particle size and shape, a wide range of densities is 
‘possible. It was found advantageous to perform the 
welding operation in a closed mold, and to incorporate 
‘in the polystyrene a blowing agent which would decom- 
“pose or vaporize during the welding operation and thus 
create positive pressure in the mold. 

_ A polystyrene material, termed expandable beads, 
was made available commercially just prior to the start 
of this program. These are discrete spheroids of poly- 


‘styrene which contain a blowing agent. When heated 


above about 85°C, the spheroids expand and the degree: 


of expansion can be controlled. 

__A series of experiments was performed to develop 
formulations and processing techniques for the desired 
range of dielectric materials. Materials were molded in 
‘small rectangular aluminum molds. When the work had 
progressed to the point that good mechanical quality 
polystyrene-based dielectrics were consistently ob- 
tained, measurements were made to determine the re- 
quired density for the material in each shell. The loss 
tangent of all these materials was less than 0.001. Fig. 
4 shows the experimentally determined variation of ¢ 
with density; for comparison, Fig. 4 includes curves for 
both a linear and a logarithmic variation of ¢ with 
density. Using these experimental data, the formula- 
tions and densities for each shell were determined. The 
hemispherial shells were then cast and trimmed for 


assembly. A set of such shells, comprising one half of a- 


lens, is shown in Fig. 5. The lens was then assembled, 
staggering the seams between adjacent shells and ce- 
menting only the outside layer together. The completed 
lens is shown, mounted for test, in Fig. 6. The completed 
lens is of 17.54-inch diameter, where the outer layer 
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Fig. 5—Exploded view of lens shells. 


Fig. 6—Completed spherical luneberg lens. 


(e=1, 18-inch diameter) may be thought of as being 
comprised of air. 

Microwave Data: Three 18-inch diameter spherical 
Luneberg lenses were constructed. Materials for the 
first lens were formulated on the basis of the linear curve 
in Fig. 4 and therefore the dielectric constants of its 
shells are too high. Materials for the second and third 
lenses were formulated on the basis of the experimental 
curve in Fig. 4. The third lens was delivered directly to 
Air Force Cambridge Research Center. 

H-plane radiation pattern measurements were made 
at three wavelengths: 3.2 cm, 1.8 cm, and 0.8 cm. A 
typical radiation pattern of the second lens at \=3.2 
cm is shown in Fig. 7. Summary data for a number of 
patterns for the first and second lenses are given in 
Table I. Patterns were taken for 18 equally spaced feed 
positions on a circumference of the lens; in some cases, 
patterns were taken for feed positions on three orthog- 
onal planes. Illumination taper refers to the free space, 
primary radiation pattern taper of the flared waveguide 
at +90°, the taper being the same in the E and H planes. 
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Fig. 7—Typical radiation pattern for spherical Luneberg lens. 


TABLE I 
RADIATION PATTERN DATA ON SPHERICAL LUNEBERG LENSES 


Beam- 


: Beam ; Side Lobe 
Lens nation : width No. of 
No » (cm) Taper rien Factor rat Patterns 
(db) “8. (deg) 
1 See 10 4.3 61 sha? 54 
72 A} 74 10 4.0 57 14.1 18 
2 Ale? ‘20 4.4 63 slic! 54 
2 Sin 30 4.8 68 19.4 18 
2 1.8 20 2.6 66 14.9 18 


Beamwidths and sidelobe levels are presented as aver- 
age values in Table I. The beamwidth factor, defined as 
beamwidth Xdiameter/wavelength, was computed in 
each case on the basis of the average beamwidth and an 
18-inch diameter, although the diameter of shell No. 
9 is 17.55 inches. Only the second lens was measured at 
\=0.8 cm; it did not focus for any feed position at this 
wavelength. The design radial feed position was 0.228 
inch from the outer edge of shell No. 9 (the outer edge 
of shell No. 10). For the first lens, the radial feed posi- 
tion should have been slightly inside shell No. 9; in 
order not to destroy the lens, measurements were taken 
with the feed at the outer edge of shell No. 9 with the 
radiation patterns slightly defocused. For the second 
lens, measurements were taken at the feed position for 
best focused patterns, displaced } inch from the outer 
edge of shell No. 9, which corresponds quite well to the 
design focal position. 

The measured gain for the second lens with 20 db il- 
luminated taper at \=3.2 cm is 31.0 db. The gain com- 
puted from 27,000/626:, where 6% and 6y are the E-and 
H-plane beamwidths respectively, assuming beam- 
widths of 4.4°, is 31.5 db. The gain computed from 
47rA/)? is 33.0 db. If the measured gain is equated to 
k4rA/)*, k=0.63, where & is the aperture efficiency. 

Phase fronts from the second lens were measured, 
along a tangential aperture, by standard techniques at 
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Fig. 8—Typical phase front from second spherical lens. \=3.2, 
20 db illumination taper. 


TABLE II 
DIELECTRIC CONSTANTS OF SHELLS IN SECOND SPHERICAL LENS 


Shell No. € Specified 
0 1.85-+0.18 2.0 
—0.20 
1 1.90+0.06 Lg 
—0.03 
2 1.89-+0.02 1.8 
—0.03 
3 1.70+0.03 i7 
—0.1i1 
4 £61 1-0) 08 1.6 
5 1.51+0.02 1.5 
6 1.41+0.06 1.4 
—0.09 
7 1.31-++0.07 igs 
—0.03 
8 1.22+0.03 i leg’) 
—0.02 
9 1.25+0.07 ha | 
—0.14 


A=3.2 cm for a number of feed positions. A typical 


phase front is shown in Fig. 8. The deviations from a 
linear phase front are approximately three times as large 
as those for the two-dimensional lenses shown in Fig. 2. 
Assuming the same phase front deviations at 1.8 and 0.8 
cm, the deviations at these wavelengths are approxi- 
mately 0.4 and 0.9), respectively. These data correlate 
quite well with the radiation pattern deterioration at 
the shorter wavelengths. 

It was desirable to obtain some idea of the point-to- 
point dielectric constant variations in the second lens 
shells. Only one-half of each shell was used for 
sampling. From each hemisphere used, four samples 
were cut from positions on a great circle: two samples 
were taken almost diametrically opposite each other, 
one sample at the center of the hemispherical shell, and 
one at an intermediate position. Average dielectric 
constants and the total deviations are given in Table II, 
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Since the number of samples was small, these deviations 
do not necessarily represent the total deviation in each 
shell. 

It is apparent from the above data that the random- 
ness of the dielectric constant variations throughout the 
lens would permit relaxation of the dielectric constant 
tolerance +0.02 for an 18-inch-diameter lens at \ =3.2 
cm. However, it is also apparent that dielectric constant 
deviations as large as those found in the second lens are 

‘noticeable at 1.8 and 0.8 cm. In the absence of a rigorous 
theoretical study of stepped-index Luneberg lenses for 
comparison, it is difficult to determine what amount of 
performance deterioration is due to the stepping tech- 
nique itself and what amount is due to dielectric con- 
stant variations. 


CONCLUSION 


In the absence of materials to use for constructing 
continuous variation-of-index Luneberg lenses, the 
technique of using a stepped variation of index was in- 
vestigated. Data on two-dimensional, TE» mode, Lune- 
_ berg lenses show that with a reasonable number of steps 
good lens performance can be obtained, although as ex- 
pected the performance is not as good as from continu- 
ous index lenses. Expanded polystyrehe materials were 
developed with dielectric constants from 1.1 to 2.0 in 
increments of 0.1. 10-step spherical Luneberg lenses 
were constructed by assembling molded hemispherical 
shells of these materials. 
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The quality control on the expanded polystyrene ma- 
terials at present appears adequate for small lenses but 
needs improvement for extension to larger lenses. It is 
felt that the technique of using molded hemispherical 
shells lends itself to accurate quality control in large 
quantity production where much of the hand labor 
would be eliminated. 

The expanded polystyrene materials appear reason- 
ably weatherproof and may be satisfactory for limited 
outdoor use. However, a thin, weatherproof cover which 
will absorb ultraviolet radiation is probably desirable 
for prolonged outdoor applications. These materials are 
dimensionally stable although resilient, and seem ca- 
pable of surviving normal mechanical abuse. They do not 
show a tendency to crack and craze as does solid poly- 
styrene. For high-temperature applications, expanded 
materials, such as silicones, teflon, or ceramics, should 
be developed. For applications requiring lightweight 
lenses (an 18-inch diameter Luneberg lens of expanded 
polystyrene weighs 35 pounds) the use of loaded, light- 
weight foams is indicated. 
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Measurements of 1250-MC Scatter Propagation as 
Function of Meteorology* 
D. L. RINGWALT}, W. S. AMENTT, anp F. C. MACDONALD} 


Naval Research Laboratory and the University 

of Florida conducted a 1250-me scatter propaga- 
tion experiment on a 262-nautical-mile overwater path, 
with transmitter at Cape Canaveral, Fla., and receiver 
at Nassau, Bahamas. Airborne refractometer soundings 
were taken at four selected points along the path every 
six hours. Meteorological observers in the aircraft re- 
corded cloud formation, haze-layer heights, etc. Six- 
hourly radiosonde and pilot balloon runs were made at 
Patrick AFB, Grand Bahama (near midpath), Eleu- 
thera, and Mayaguana. The received signal was moni- 
tored and recorded continuously at the Florida Field 
Station and Nassau, and NRL’s Flying Laboratory 
measured field strength in the air. 

The purpose of the program was to ascertain whether 
there is a correlation between low-level refractometer 
soundings (500 to 15,000 feet) and the median signal 
power received on the ground-to-ground link. One hun- 
dred and five full-length refractometer soundings were 
taken during the week. Each profile showed a break of 
about 40 or 50 ” units in a 1000-foot interval centered 
at altitudes from 4000 to 6000 feet, moist oceanic air of 
high refractive index being found below the break. This 


1D) sera the week of December 3-8, 1956, the 


[ * Manuscript received by the PGAP, March 15, 1957. 
} Wave Propagation Branch, Electronics Div., Naval Res. Lab., 
Washington 25, D. C. 


break corresponded with a visible haze layer, a temper- 
ature inversion, and the base of scattered’ cumulus 
clouds. The clouds, and probably the mixing of air at 
the inversion height, produced large random excursions 
of the refractometer outputs. The larger cumulus for- 
mation rose to over 13,000 feet and may be regarded as 
imposing a 10-—40-mile horizontal meso-meteorological 
scale on the path meteorology. As a result, the 105 
refractometer soundings may be regarded as roughly 
independent samples of a single meteorological situation 
having a dominant 40—50 n-unit break. 

Supporting this view, the synoptic weather charts 
throughout the week showed a stationary high-pressure 
area off Cape Hatteras and a stationary, dissolving cold 
front eastward of Cuba. 

Minor aspects of the soundings were: At altitudes 
13,000-14,000 feet, a thin 4 or 5 n-unit layer was gen- 
erally encountered. This layer is believed to represent 
the horizontal leakage of water vapor from the higher 
cumulus clouds. Apart from this layer and away from 
clouds, the air above the inversion was featureless, al- 
though there were occasional wedge-like intrusions, just 
above the inversion, of magnitude of 5-6 n-units. Han 
diurnal trend (pointed out by Cullen Crain) occurred in 
the inversion height at midpath; the inversion height 
was maximum in the afternoon and minimum in the 
early morning. (No corresponding diurnal effect in sig- 
nal was observed.) 
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The 1250-mc signal received at the Florida Field 
Station showed Rayleigh fading, characteristic of a 
scatter-propagated signal, at all times. Hourly medians 

_of the signal were substantially constant throughout the 
week except for the first 18 hours of December 5. Dur- 
ing this period the hourly medians averaged about 10-db 
higher (two medians were 24-db higher). Since no major 
feature in the refractometer soundings correlated with 
this signal enhancement, the wind shear across the in- 
version was investigated. Pibal data, available at 1000- 
foot height intervals, and radiosonde data gave wind 
vectors of considerable scatter at inversion height 
(possibly due to convective mixing there). Hence, wind 
shear across the inversion was estimated from the rela- 
tively stable wind data at heights roughly 1500 feet 
above and below the inversion. Midpath wind shear 
values were high during the high signal period. When 
medians of the hourly signal medians for each of fifteen 
six-hour periods and the corresponding wind shear value 
at midpath are compared by the rank-order method, 
odds three in four, that the correlation of wind shear 
with signal is not accidental, are obtained. However, 
wind shear data from Patrick AFB and Eleuthera, near 
the path’s endpoints, show little correlation with signal 
level and confirmation of this wind shear signal correla- 
tion is needed. 

According to original plans, the Flying Laboratory 
was to fly at 1000 feet from Nassau to Cape Canaveral, 
measuring the dependence of signal on range. The re- 
turn flight to Nassau was to be devoted to refractometer 
soundings with spot signal observations. Those spot 
observations indicated a signal maximum at inversion 
height. This fact (and the failure of attempts to corre- 
late the ground-to-ground signals with the refractometer 
soundings) caused a reassignment of the efforts of the 
Flying Laboratory in favor of more frequent measure- 
ment of signal vs height. 

The inbound runs at 1000 feet showed a Rayleigh 
distributed signal, roughly a three-second interval be- 
tween median crossings, and a range-attenuation rate 
of about 0.17 db per nautical mile. With height, the sig- 
nal increased 5—15 db in the 1000 feet below the inver- 
sion height, diminished 5 to 15 db in the next 1000 feet, 
and then increased with altitude. Near and above the 
inversion, the Rayleigh fading was usually (four out of 
six times) superimposed on a slower fade having a period 
of about one minute with fades of about 10 db. These 
features occurred at ranges beyond 160 nautical miles. 
At closer ranges, a rapid increase of mean signal and 
shallow fading showed that the dominant propagation 
mechanism was diffraction and/or ducting. 
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The attenuation rate of 0.17 db/mile is in agreement 
with much other work. The height gain data (as well as 
the suggested shear correlation) suggest two alternative 
scatter mechanisms. 

1) The consistent 40-50 n-unit break at inversion 
height could have acted as a curved partially reflect- 
ing ceiling, with the signal beneath arriving by multiple 
bounces. 2) This 40-50 n-unit break constituted an 
elevated duct about 1000 feet thick having about 
a 10 m-unit (maximum) n-deficit. This erratically de- 
fined, and hence randomly leaky, duct could have been 
coupled with both receiver and transmitter via scatter- 
ing at either end of the path. 

The slow underlying fade in the signals at inversion 
height suggests the duct mechanism, while the midpath 
shear correlation supports the leaky ceiling picture. 

An additional effect of low-level meteorology was 
noted: A subsidiary haze layer and cloud deck occurred 
on December 7. This layer was at an altitude of 1500- 
2500 feet over Nassau and extended about 70 miles 
along the path, and appeared as an 8-12 n-unit “inver- 
sion” in the corresponding soundings. It served to en- 
hance the ground-to-air signal over Nassau at the ex- 
pense of the ground-to-ground signal. 

Among the conclusions unlikely to be altered by con- 
tinuing analysis of the data are: the prominent 18-hour, 
10-db increase in ground-to-ground signal did not corre- 
late with the refractometer soundings in the 500—10,000- 
foot height region. High signal periods correlate sug- 
gestively with periods of wind shear across the domi- 
nant 40-50 m-unit inversion at 4000-6000 feet. High sig- 
nals with slow underlying fades, observed at inversion 
altitudes, indicate ducting or multiple scattering as the 
dominant mechanism for long-range microwave propa- 
gation in this meteorological regime. The refractometer 
profiles were varied in detail because of the meso- 
meteorological scale imposed by cumulus clouds pro- 
truding through the inversion at scattered points. 

The University of Florida work was performed under 
ONR contract Nonr 580(04). Communications, me- 
teorological observations, and other valuable support 
was provided by the Air Force Missile Test Center, 
Patrick Air Force Base. The Naval Air Development 
Unit, South Weymouth, Mass., gave outstanding sup- 
port in flying a portion of the refractometer schedule. 
Valuable advice and aid in data interpretation were 
given by volunteer participants from Lincoln Labo- 
ratory, Air Force Cambridge Research Center, and the 
Applied Physics Laboratory. Certain Lincoln Labo- 
ratory equipment was vital in obtaining the ground-to- 
air signals. 
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Design Data for Small Annular Slot Antennas* 
W. A. CUMMING} anp M. CORMIER} 


Summary—While several excellent papers have been published 
describing the theory of operation of the annular slot radiator, certain 
design parameters have not been too readily available. This note 
gives susceptance data for practical configurations having electrical 
diameters of from 0.2 to 0.6 wavelength. 


NNULAR slot antennas have a number of air- 
A borne and vehicular applications in the vhf and 
uhf bands. From the theoretical viewpoint this 
type of radiator has been treated as an open-ended co- 
axial line with an infinite conducting flange attached to 
the outer conductor.!~ Various field configurations can 
be set up in the annular aperture thus formed, but since 
most applications require an omnidirectional E-plane 
pattern, dominant-mode excitation is most commonly 
used. Different methods of achieving this excitation 
have been proposed,®* but all basically involved excita- 
tion from a radial transmission line, or a conical trans- 
mission line, or a combination of the two. Excitation by 
a radial transmission line appears to be the most prom- 
ising method when space limitations are severe, and 
when the bandwidth requirement is not excessive. Slot 
susceptance data were therefore obtained with refer- 
ence to this type of feed. 

Following the method of Johnson and Rothe® the di- 
electric-covered annular slot is backed on the shadow 
side of the ground plane by a radial transmission line, 
short circuited at its inner radius. This system is tuned 
to resonate in the dominant mode, and matching to a 
coaxial line source is accomplished through a magnetic 
coupling loop, as shown in Fig. 1. The usual design pro- 


Fig. 1—Annular slot fed by radial transmission line. 


* Manuscript received by the PGAP, August 14, 1957; revised 
manuscript received, January 20, 1958. 
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Ottawa, Can. 
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cedure is to tune the cavity to resonance by matching 
the radial line susceptance to the slot susceptance at the 
inner slot boundary a. The series slot resistance R,’ is 
then matched to the characteristic impedance of the 
coaxial transmission line by suitably adjusting the mu- 
tual impedance Xm, and finally the loop reactance is 
tuned out by means of a series tuning capacitor C. An 
equivalent circuit of this configuration is shown in 
Fig. 2. 

Slot conductance G, can readily be calculated from 
the field configuration in the slot, and is given graphi- 
cally in Fig. 3. The susceptance of the radial transmission 
line B; can also be calculated, and design curves of this 
quantity are presented in Fig. 4. The mutual impedance 
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X,, of the coupling loop is given by 
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Finally, measured values of slot susceptance B,, re- 
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ferred to the inner radius of the slot are given in Figs. 
5, 6, and 7, for slots of various widths and dielectric 
coverings. These data were obtained by placing both 
an input and output coupling loop in the radial trans- 
mission line, and using it as a lossy transmission device. 
At the resonant frequency of the system, the input sus- 
ceptance of the radial transmission line was calculated 
and equated to the slot susceptance. This procedure was 
carried out for a number of slot widths and radial trans- 
mission line dimensions. 
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Kansas, Lawrence, 
Kan., and the Ph.D. 
degree in engineering 
from Yale Univer- 
sity, New Haven, 
Conn., in 1944, 

From November, 
1943 until July, 1957, 
he was a member of 
the faculty of the 
University of Texas, 
Austin, Tex., the last three years as professor 
of electrical engineering. In addition, he 
assisted in the organization of the electrical 
engineering research laboratory of the Uni- 
versity of Texas in 1946 and was assistant 
director of the laboratory from 1948. He 
has worked primarily in the fields of low 
angle tropospheric radio propagation, an- 
tenna systems, and analog computers for 
data reduction. 

In July, 1957, he joined Collins Radio 
Company, Dallas, Tex. as assistant director 
of research and development. 

Dr. Brooks is a member of Sigma Xi, Tau 
Beta Pi, and Eta Kappa Nu. 


F. E. Brooks, JR. 
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E. Gésta R. Carlson was born on January 
12, 1930, in Vadstena, Sweden. He entered 
the Royal Institute of Technology, Stock- 


holm, Sweden in 1948, where he received the 
degree of Civil-ingenjérsexamin in electrical 
engineering in 1954. ; 


Since 1952 he has 
been with the Re- 


National. 
Stockholm, as a re- 
search engineer, where 
he participates in 
wave propagation re- 
search especially con- 
cerning beyond-the- 
horizon phenomena. 

Mr. Carlson is a 
member of Svenska 
Teknologféreningen. 


E. G. R. CARLSON 
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Robert L. Carrel (S’54) was born in Fort 
Wayne, Ind., on February 8, 1933. He re- 
ceived the B.S. degree in electrical engi- 
neering in 1955 from 
Purdue University, 
Lafayette, Ind., and 
the M.S. degree in 
electrical engineering 
from the University 
of Illinois, Urbana, 
in 1957. 

He was employed 
by the Farnsworth 
Electronics Company 
from February to 
September, 1955, as 
an antenna engineer 
in the microwave section. He entered the 
University of Illinois in September, 1955, as 
a research assistant in the Antenna Labora- 
tory, and is currently working toward the 
Ph.D. degree. 

Mr. Carrel is a member of Eta Kappa 
Nu and Tau Beta Pi. 
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H. P. Coleman was born in November, 
1928, at Paeonian Springs, Va. He studied 
at the College of William and Mary, 

Williamsburg, Va., 
and at American 
University, Washing- 
ton, D.C. 

‘ In 1945, he joined 
the U. S. Naval Re- 
search Laboratory, 
Washington, D. C.as 
a summer employee, 
and was engaged in 
vacuum tube devel- 
opment and later in 
accelerator develop- 
ment. Mr. Coleman 

became a member of the Microwave An- 

tennas and Components Branch, Electronics 

Division in 1952. Since then he has worked 

in the field of microwave antennas, prin- 

cipally with scanning systems. 
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K. Folke Eklund was born on July 17, 
1924, in Vasteras, Sweden. After preparatory 
education and military service he entered 
the Chalmers Insti- 
tute of Technology, 
Gothenburg, in 1945 
where, in 1949 he re- 
ceived the Civilin- 
genjor degree in elec- 
trical engineering. 

Since 1949 he has 
been with the Re- 
search Institute of 
National Defense, 
Stockholm, where he 
is a section leader in 
the laboratory for an- 
tennas and propagation. 

He has been concerned mainly with 
_ microwave propagation problems. 

Mr. EklundJis a member of Svenska 
Teknologféreningen. 


h K. F. EkLuND 
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3 John P. German (A’51-SM’55) was born 
in Livingston, Tex., in 1917. He received the 
_ B.S.E.E.degreefrom The University of Texas, 
: Austin, Tex., in 1940, 
and spent the next 
two years working in 
industry. He spent 
three and a half years 
in the U. S. Army, 
serving in the India- 
Burma theater of op- 
erations as a radar 
officer. After World 
War II he returned 
to The University of 
z J. P. GERMAN Texas for graduate 
4 work in electrical en- 
_ gineering. He received the M.S.E.E. degree 
_ in 1949 and the Ph.D. degree in 1955. Dur- 
_ ing this time he also served as a teacher and 
- as a radio engineer for the electrical engi- 
neering research laboratory at The Uni- 
_ versity of Texas. In 1955 he accepted ap- 
_ pointment as associate professor of electrical 
engineering at Purdue University, West 
_ Lafayette, Ind. 

Dr. German is a member of Eta Kappa 
_ Nu, Sigma Pi Sigma, and Sigma Xi. 
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| Bengt A. S. Josephson (M’55) was born 
June 28, 1914, in Stockholm, Sweden. He 

received the Civilingenjérsexamen degree in 
electrical engineering 
in 1942 and the Tek- 
nologie Licentiatexa- 
men degree in 1947 
at the Royal Institute 
of Technology in 
Stockholm, where he 
also worked as lab- 
oratory and research 
assistant from 1940 
to 1943. 

During 1944 and 
1945 he was employed 
; with Telefonaktie- 
- bolaget L. M. Ericsson in Stockholm, and 

was mainly engaged there in development of 


B. JOSEPHSON 
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uhf measuring instruments. He then joined 
the Research Institute of National Defence 
(FOA) in Stockholm, where he is in charge of 
the laboratories for transmission lines, an- 
tennas and propagation. 

In 1946, he was granted a fellowship 
by the American Scandinavian Foundation 
and spent half a year with the RCA Labora- 
tories in Princeton, N. J., and Long Island, 
N. Y. Since 1947 he has been part-time 
teaching uhf techniques at the Royal Insti- 
tute of Technology and since July, 1957 
he has been professor pro tempore in micro- 
wave engineering (during leave from 
FOA). 

Mr. Josephson has been active in Swedish 
and International standardization work on 
rf transmission lines and is chairman of a 
working group on waveguides within the 
International Electrical Commission. 


George D. M. Peeler (M’52) was born in 
Rockwell, N. C. on January 10, 1923. He 
received the A.B. degree cum laude in phys- 
ics in 1943 from 
Catawba __ College, 
Salisbury, N. C. and 
the M.A. degree in 
engineering science 
and applied physics 
in 1948 from Harvard 
University, Cam- 
bridge, Mass. 

He taught physics 
for one year at the 
University of North 
Carolina, Chapel Hill, 
and also attended the 
University of Maryland, College Park. He 
was employed by the Naval Research Labo- 
ratory in Washington, D. C. asa member of 
the Search Radar Branch from 1944 to 1950 
and as a member of the Microwave An- 
tenna and Components Branch from 1950 to 
1956. Since February, 1956, he has been head 
of the Antenna Branch in the Missile Sys- 
tems Division of Raytheon Manufacturing 
Co., Bedford, Mass. 

Mr. Peeler is a member of the American 
Physical Society. 


G. D. M. PEELER 
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Paul I. Richards (A’45) was born in 
Orono, Maine, on February 8, 1923. He re- 
ceived the M.A. degree in 1946 and the 
Ph.D. degree in phys- 
ics in 1947 from Har- 
vard University, 
Cambridge, Mass., 
with a thesis on the 
analytic theory of 
impedance-functions. 

From 1947 to 1952 
he was a physicist at 
Brookhaven Nation- 
al Laboratory, where 
he worked on theoret- 
ical physics and on 
the experimental de- 
velopment of a time-of-flight mass spectrom- 
eter. 
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After spending a year at Transistor Prod- 
ucts, Inc., in 1954 he joined Technical Opera- 
tions, Inc., where he is presently employed in 
theoretical analysis of various problems in 
physics, engineering, and operations research, 
mainly under the sponsorship of various 
government agencies. His recent contribu- 
tions have centered primarily in the trans- 
port theory of nuclear and optical radiation. 

Dr. Richards is a member of the Ameri- 
can Physical Society, Sigma Xi, and Phi 
Beta Kappa. 
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Dipak L. Sengupta (S’55) was born in 
Bengal, India in 1931. He received the B.Sc. 
degree with honors in physics in 1950 and the 
M.Sc. degree in radio- 
physics and_ elec- 
tronics in 1952 both 
from Calcutta Uni- 
versity. He worked 
at that institution as 
a Government of 
India senior research 
scholar until Sep- 
tember, 1954. Dur- 
ing that time he de- 
signed and developed 

D. L. SENGUPTA a special type of 

transmitter and re- 

ceiver suitable for super-refraction studies 
at 300 me. 

At the present time, Mr. Sengupta is in 
the last year of graduate study at the Elec- 
trical Engineering Department of the Uni- 
versity of Toronto. He is also working as a 
part-time assistant on the radioastronomy 
project in the Electrical Department, where 
he is responsible for the development of the 
antenna assembly of the recently installed 
radio telescope. His main interests are in em 
theory, microwave antennas, and microwave 
theory. 
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Hillel Unz (S’53-M’57) was born in 
Darmstadt, Germany, on August 15, 1929. 
He went to Israel in 1931 and served in the 
Israel Defense Forces 
during the war of in- 
dependence in 1948- 
1949. He received the 
B.S. degree in elec- 
tronics, summa cum 
laude, from Israel In- 
stitute of Technol- 
ogy, in 1953. 

From 1953 to 1956, 
he was a research as- 
sistant in the micro- 
wave antenna labora- 
tory (Office of Naval 
Research Project) in the Electrical Engineer- 
ing Department at the University of Cali- 
fornia, Berkeley, Calif., where he received 
the M.S. degree in 1954, and the Ph.D. de- 
gree in 1957. 

Since January, 1957, he has been an 
assistant professor of electrical engineer- 
ing at the University of Kansas, Lawrence, 
Kan. His main interest is in the fields of 
electromagnetic theory, boundary value 
problems, diffraction theory, and antennas. 

Dr..Unz is a. member of Sigma Xi. 


H. Unz 
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PEOPLE... 


like the results of their work to be 
elegant in the mathematical sense of 
“ingenious simplicity,” and we are 
no exception. In fact, we expend a 
great deal of effort in search of this 
kind of simplicity. Applied to air- 
borne antennas, “elegant” obviously 
implies minimum size and weight. 
These characteristics are equally 
important to us. 


With many of our antennas we be- 
lieve that a considerable amount of 
elegance has been achieved. This is 
true of the family of airborne uhf 
monopole type radiators which we 
manufacture. It includes a wide va- 
riety of stubs and annular slots, all 
of which are substantially simpler 
and smaller than can be obtained 
through application of the published 
art. Nevertheless, theoretical con- 
siderations indicate that further sub- 
stantial reductions in size are not 
impossible. 


Efforts to discover elegant means 
of accomplishing this size reduction 
will be a small but important part of 
our activity in the future. We con- 
sider this to be a most interesting 
problem and one to which no exten- 
sion of conventional techniques may 
be expected to yield a solution. 


We are anxious to add a few sen- 
ior engineers to our staff. Our labo- 
ratory is modern and air-conditioned 
and our program is very varied. We 
have about 100 people, and so far 
as we can determine, we have de- 
signed antennas for more different 
aircraft than any other group in the 
country of any size. If you are 
anxious to do original work in this 
field, please contact Arthur Dorne, 
President. 


DORNE & MARGOLIN, unc. 
29 NEW YORK AVENUE 
WESTBURY, N.Y. 


Conta 


MR. DELMER G3 PORTS 
Jansky and Bailey, Inc. 
13839 Wisconsin Ave, N.W. 
Washington 7, D. C. 
Telephone: 
Federal 3-4800 
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AVAILABLE BACK ISSUES OF IRE TRANSACTIONS 
ON ANTENNAS AND PROPAGATION 


Publication Prices 
Group IRE Non- 
Members Members Members* 

AP-1, No. 1, July, 1953 $1.20 $1.80 $3.60 
AP-1, No. 2, October, 1953 $1.20 $1.80 $3.60 
AP-2, No. 1, January, 1954 $1.35 $2.00 $4.05 
AP-2, No. 2, April, 1954 $2.00 $3.00 $6.00 
AP-2, No. 3, July, 1954 $1.50 SV As) $4.50 
AP-3, No. 4, October, 1954 $1.50 $2.25 $4.50 
AP-3, No. 1, January, 1955 $1.60 $2.40 $4.80 
AP-3, No. 2, April, 1955 $1.60 $2.40 $4.80 
AP-4, No. 4, October, 1956 $2.10 $3.15 $6.30 
AP-5, No. 1, January, 1957 $3.20 $4.80 $9.60 
AP-5, No. 2, April, 1957 S145 $2.60 $5020 
AP-5, No. 3, July, 1957 $2.00 $3.00 $6.00 
AP-5, No. 4, October, 1957 $1.70 $2.25 $5.10 
AP-6, No. 1, January, 1958 $3.85 $5.80 $11.55 


* Colleges, Universities, Subscription Agencies, and all Libraries, may 
purchase at TIRE Member rate. 
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The IRE Professional Group on Antennas and Propagation is grateful for = | 
the assistance given by the firms listed below, and invites application for ; . 
Institutional Listing from other firms interested in the field of Antennas and ~ t 
Propagation. 

ANDREW CORPORATION, 363 E. 75th St., Chicago 19, Ill. be 


ANTLAB, INC., 6330 Proprietors Rd., Worthington, Ohio 
Antenna Pattern Range Systems—Recorders & Mounts. 


BLAINE ELECTRONETICS, INC., 14757 Keswick St., Van Nuys, Calif. 


Antennas, Antenna Systems, Transmission Lines, Development and Production. 
Antennas, Paraboloids, Scale Models, Antenna Radiation Pattern Measurement Towers. | 
| 


COLLINS RADIO COMPANY, Cedar Rapids, Iowa | 
Antenna Design and Propagation Research Related for Airborne and Ground Communication Systems. 


DEVELOPMENTAL ENGINEERING CORP., 1001 Conn. Ave. N.W., Washington, D. C. and Leesburg, Va 
Research, Development, Installation of Antennas and Antenna Equipment for Super Power Stations. 


THE GABRIEL LABORATORIES, Div. of the Gabriel Co., 135 Crescent Road, Needham Heights 94, Mass. 
Research and Development of Antenna Equipment for Government and Industry. | 


HUGHES AIRCRAFT COMPANY, Culver City, Calif. 
Research, Development, Mfr.: Radar, Missiles, Antennas, Radomes, Tubes, Solid State Physics, Computers. : 


I-T-E CIRCUIT BREAKER CO., Special Products Div., 601 E. Erie Ave., Philadelphia 34, Pa. 
Design, Development and Manufacture of Antennas, and Related Equipment. 


JANSKY & BAILEY, INC., 1339 Wisconsin Ave. N.W., Washington 7, D.C. 
Radio & Electronic Engineering; Antenna Research & Propagation Measurements; Systems Design & Evaluation. 


MARK PRODUCTS CO., 6412 W. Lincoln Ave., Morton Grove, IIl. 
Multi Element Grid Parabolas, Antennas for Two-Way Communications, R & D. 


MARYLAND ELECTRONIC MANUFACTURING CORPORATION, College Park, Md. 
Antenna and System Development and Production for Civil and Military Requirements. 


THE RAMO-WOOLDRIDGE CORPORATION, Los Angeles 45, Calif. 


TRANSCO PRODUCTS, INC., 12210 Nebraska Ave., Los Angeles 25, Calif. a 
Res., Design, Dev., & Mfr. of Antenna Systems & Components for Missile, Aircraft & Ground Installations. 


WHEELER LABORATORIES, INC., 122 Cutter Mill Road, Great Neck, N. Y. 
Consulting Services, Research and Development, Microwave Antennas and Waveguide Components. 


The charge for an Institutional Listing is $25.00 per issue or $75.00 for four 
consecutive issues. Application may be made to the Technical Secretary, 
The Institute of Radio Engineers, 1 East 79th Street, New York 21, N.Y. 


